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EXECUTIVE  SUMBiARY 


The  general  emphasis  for  1994  would  be  on  increased  software 
development,  testing  of  subelements  and  design  calculations.  For  these 
purposes,  the  constitutive  law  coding  and  development  would  be 
coordinated  by  Nick  Aravas.  and  implemented  in  ABAQUS.  The  initial 
implementation  would  be  the  elastic/plastic  model  for  MMCs  with  interface 
debonding  developed  in  1993  (Leckie).  This  would  be  extended  in  1994  to 
include  creep  and  some  aspects  of  thermomechanical  cycling.  The  code 
would  be  used  for  design  calculations  concerned  with  MMC  rotors, 
actuators  and  vanes  (Leckie).  A  plan  is  being  formulated  to  collaborate  with 
Pratt  and  Whitney  to  acquire  MMC  sub-elements  representative  of  these 
components  during  1994.  E}q>erimental  tests  on  these  subelements  would 
be  capable  of  providing  a  direct  validation  of  the  code  capabilities. 

Constitutive  law  and  fatigue  llfing  software  would  be  created  for  CMCs 
using  continuum  damage  mechanics  (CDM)  approaches  (Leckie. 
McMeeklng).  The  approach  has  been  motivated  by  micromechanics  models 
developed  in  1993  (Hutchinson,  Zok.  E^rans).  These  codes  would  be  used  to 
calculate  stress  redistribution  effects  and  fatigue  life  on  simple  sub¬ 
elements.  such  as  center  notched  and  pin-loaded  plates.  Comparison  with 
e^q^erimental  measurements  needed  to  test  the  fidelity  of  the  models  will  be 
based  on  moir6  interferometiy  and  thermoelastic  emission.  This  effort  is 
coordinated  with  the  NASA  EPM  program  throu^  both  General  Electric  and 
I^tt  and  Whitney.  A  plan  for  acquiring  sub-elements  fiom  DuPont  Lanxide 
is  being  formulated. 

A  new  emphasis  for  1994  would  be  on  the  transverse  properties  of 
CMCs.  The  measurements  and  calculations  performed  in  1993  have 
indicated  a  strategy  for  curved  sections  and  Junctions  that  would  establish  a 
consistent  design  approach.  The  basic  approach  for  resisting  failures  from 
combinations  of  interlaminar  shear  and  transverse  tension  involves  the  use 
of  stitching  and  angle  ply  weaving  patterns  that  inhibit  major  reductions  in 
stiffiiess  when  matrix  cracks  are  induced  transverse  loads  and  bending 
moments.  For  this  purpose,  calculations  would  be  performed  that  combine 
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the  mechanics  of  delamination  cracks  with  models  of  bridging  by  inclined 
fiber  bundles  (Hutchinson.  Ashby.  Evans.  McMeeking).  The  insight  gained 
firom  these  calculations  would  be  used  to  design  and  acquire  sub-elements, 
such  as  C  sections  and  T  Junctions. 

Additional  software  development  will  be  for  creep  and  creep  rupture 
(McMeeking.  The  models  devised  in  1993  and  test  data  relevant  to  MMCs 
will  be  combined  into  a  code  that  predicts  the  creep  and  rupture  of 
unidirectional  MMCs  subject  to  multiaxial  loads.  Some  aspects  of  this  code 
will  also  be  applicable  to  CMCs. 

Two  new  activities  will  be  introduced  in  1994:  thermal  properties  and 
damping.  The  thermal  properties  will  be  studied  on  both  CMCs  and  MMCs 
(Ashby.  Hutchinson).  Measurements  of  thermal  dlfiusivity  will  be  made  by 
the  laser  flash  method  and  related  to  the  properties  of  the  interface  and  the 
density  of  matrix  damage  in  the  material.  Thermal  e}q)ansion  measurements 
will  also  be  performed  with  emphasis  on  determining  hysteresis  effects, 
which  can  be  related  to  the  temperature  dependence  of  the  interfaces 
properties,  through  cell  models.  The  latter  mi^t  evolve  into  a  diagnostic  for 
establishing  relationships  between  the  interface  properties  and 
thermomechanical  fatigue. 

The  processing  activities  in  the  program  will  have  nevdy  established 
goals  in  1994.  The  principal  emphasis  will  be  on  concepts  for  affordable 
manufacturing.  The  issues  selected  for  investigation  will  be  consistent  with 
manufacturing  processes  that  allow  near-net  shape  consolidation  while  still 
yielding  reasonable  combinations  of  longitudinal  and  transverse  properties. 
Performance  models  developed  in  the  program  would  be  used  as  an  initial 
test  of  concept  viability. 

Beyond  these  general  trends,  specific  activities  are  planned  for  1994. 
These  are  elaborated  below.  The  status  of  understanding  and  development 
in  each  of  these  areas  is  summarized  in  Table  1.  Increasing  msignitudes 
between  0  and  1  designate  a  knowledge  range  from  limited  to 
comprehensive. 
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TABLE  IB 

Status  of  Design  Knowledge  for  CMCs 
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2.  CONSmTmVB  LAWS 

Two  approaches  will  be  used  to  create  a  formulation  capable  of 
representing  the  in-plane  properties  of  CMCs.  One  would  be  based  on 
Continuum  Damage  Mechanics  (COM)  (Leckie).  The  other  would  use 
concepts  analogous  to  those  used  in  plasticity  theory  (Hutchinson).  The 
CDM  approach  uses  damage  parameters  that  relate  explicitly  to 
micromechanics  models.  A  potential  function  has  already  been  identified  as 
the  state  variable  which  separate^  represents  the  strain  from  the  elastic 
compliance  change  caused  by  the  matrix  cracks  and  the  inelastic  strains 
associated  with  the  debonding  and  sliding  interfaces.  Derivatives  of  the 
potential  with  regard  to  strain  and  damage  give  the  relationships  between 
variables,  such  as  stress,  interface  sliding  resistance,  matrix  crack  density, 
etc. 

The  first  version  of  the  CDM  model  would  use  the  minimum  number  of 
damage  variables  potentially  capable  of  representing  the  behavior  of 
laminated  or  woven  composites.  Cross  terms  between  the  damage  variables 
would  not  be  considered  at  this  stage.  Moreover,  matrix  cracks  would  be 
introduced  normal  to  the  maximum  principal  tensile  stress,  consistent  with 
the  e]q}erimental  observations. 

The  plasticity  theory  approach  would  seek  a  formulation  based  on 
matrix  cracks  occurring  normal  to  the  maximum  principal  tension.  It  would 
introduce  parameters  that  reflect  the  inelastic  strain  caused  by  interface 
sliding  upon  off-axis  loading  which  would  be  calibrated  from  tests  performed 
in  tension  in  0/90  and  45/45  orientations. 

The  insight  needed  to  characterize  off-axis  loading  effects  will  be  gained 
from  cell  models  (Hutchinson)  in  a  marmer  analogous  to  that  previously 
used  for  axial  loads.  The  principal  objective  will  be  to  understand  trends  in 
matrix  crack  opening  and  interface  debonding/sliding  with  applied  loads. 
The  stress  on  the  fibers  will  be  calculated  with  the  intent  of  predicting 
effects  of  loading  orientation  on  fiber  failure.  The  models  will  be  compared 
with  measurements  made  in  45/45  tension,  using  various  CMCs  (E^mis). 

Calibration  of  the  damage  parameters  for  each  material  would  be  made 
from  hysteresis  loop  measurements  in  accordance  with  procedures 
developed  in  1993.  Experimental  results  obtained  in  0/90  tension.  45/45 
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tension  and  in-plane  shear  will  be  used.  In  future  work,  it  is  hoped  that 
shear  tests  will  not  be  necessary. 

The  validation  of  the  constitutive  laws  will  be  achieved  comparing 
calculations  with  measurements  made  on  sub-elements,  especially  pin- 
loaded  holes  (E>^ans).  The  experimental  results  include  residual  strains 
obtained  Moir6  interferometry  (Fig.  2. 1).  ultimate  loads  for  either  tensile 
or  shear  failure  and  principal  strain  trajectories  delineated  by  matrix 
cracking  patterns.  Acoustic  methods  will  also  be  developed  to  probe  the 
local  values  of  the  elastic  modulus  (Clarke.  Wadley)  which  could  be 
compared  directly  with  the  CDM  predictions. 

3.  FATIGUE  LIPING 
3.1  CMCs 

A  software  program  for  isothermal  low  cycle  fatigue  (LCF)  of  CMCs. 
developed  in  1993  (Fig.  3. 1)  will  be  extended  in  1994.  The  present  program 
asserts  that  fatigue  is  associated  with  cyclic  degradation  of  the  interface 
sliding  resistance,  t,  which  can  be  characterized  by  analyzing  hysteresis 
loops  measured  periodically  during  a  fatigue  test.  With  this  methodology, 
S-N  curves  have  been  predicted  for  both  unidirectional  and  woven  0/90 
composites  tested  in  (Tclic  tension  as  well  as  changes  in  compliance  and 
permanent  strain.  Some  additional  effort  is  required  to  analyze  data  on  0/90 
laminates  in  order  to  validate  the  model  predictions.  The  extensions 
envisaged  for  1994  include  thermomechanical  fatigue  (TMF),  strain 
controlled  LCF  and  off-axis  fatigue  (Zbk.  Evans).  Experiments  are  planned 
which  would  assess  the  effects  of  temperature  cycling  and  of  inclined  fibers 
on  T  degradation,  measured  from  hysteresis  loops.  Various  cell  model 
calculations  (Hutchinson)  will  be  used  to  interpret  the  experiments.  The 
results  will  be  used  to  establish  general  rules  for  interface  degradation  in 
CMCs. 

The  off-axis  e3q)enments  will  also  give  insight  into  the  fiber  failure 
criterion  that  replaces  the  global  load  sharing  (GLS)  results  successfully 
used  for  0/90  loadings.  This  study  will  coordinate  with  the  cell  calculations 
described  above,  and  the  45/45  tensile  experiments. 

Notch  fatigue  studies  wUl  he  Initiated.  These  will  examine  cyclic  stress 
redistribution  and  notch  sensitivity  (E^rans). 
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3.2  lOfCs 

Fatigue  crack  growth  and  notch  strength  studies  in  MMCs  wiU  be 
extended  to  0/90  laminates  (Zok,  Suo).  The  e3q>erlments  concerned  with 
crack  growth  will  be  interpreted  using  crack  bridging  models.  The  utility  of 
such  models  has  been  validated  in  previous  years  through  studies  on 
unidirectional  MMCs.  It  is  envisaged  that  the  fatigue  crack  growth 
characteristics  of  the  unidirectional  and  0/90  configurations  will  be  related 
through  the  volume  fraction  of  fibers  aligned  with  the  loading  direction.  The 
notch  strength  behavior  will  also  be  interpreted  using  crack  bridging 
models.  Such  models  have  been  developed  in  1993  and  found  to  be  useful  in 
rationalizing  the  behavior  of  unidirectional  materials  (Zok.  Suo).  In  all  cases, 
the  mechanical  measurements  will  be  augmented  by  tn-situ  observations  to 
identify  changes  in  damage  mechanisms  with  temperature,  fiber 
architecture,  etc.  Plans  to  study  the  influence  of  panel  thickness  on  fatigue 
and  fracture  resistance  are  also  being  developed,  as  well  as  tests  to 
understand  the  potential  for  crack  growth  in  mixed  mode  loadings  (Hirth. 
Zok). 

Studies  of  the  TMF  response  of  MMCs  loaded  parallel  to  the  fiber  axis 
will  be  initiated  (Zok.  Leckie).  Experiments  will  evaluate  both  in>phase  and 
out-of-phase  loadings.  Models  of  load  shedding  (matrix-fibers)  will  be  used 
to  Interpret  the  hysteresis  loops  and  to  develop  fatigue  life  models  applicable 
to  low  (yde,  high  strain  TMF. 

4.  CREEP  AND  RUPTURE 
4.1  MMCs 

The  considerable  progress  made  in  1993  towards  identifying  and 
imderstanding  the  mechanisms  of  creep  and  rupture  in  unidirectional 
MMCs  containing  non-creeping  fibers  (McMeeking.  Zok)  will  be  used  to 
develop  creep  rupture  software.  The  longitudinal  creep  model  to  be  used 
incorporates  stochastic  fiber  fiiacture  and  interface  sliding  in  a  format 
amenable  to  the  prediction  of  primary  and  tertiary  creep  in  terms  of  matrix 
creep  strength.  Interface  sliding  resistance,  fiber  strength.  Weibull  modulus, 
etc.  The  concepts  would  be  visualized  in  a  rupture  mechanisms  map 

KJS«Zn94 


10 


CFIg.  4.1).  The  transverse  creep  behavior  would  include  interface  debonding, 
which  greatly  accelerates  the  creep,  leading  to  marked  anisotropy.  A 
constitutive  law  for  creep  that  includes  these  effects  will  be  developed 
(Aravas.  McMeeking). 

Additional  e3q}eriments  and  calculations  will  be  conducted  to  assess  the 
effects  of  notches  and  holes  on  creep  rupture  (Zok.  Suo).  E^erience  with 
MMCs  at  ambient  temperature  indicates  that  the  notch  sensitivity  is  largely 
dictated  matrix  properties  (i.e..  strength  and  ductility).  The  reduction  in 
matrix  properties  at  elevated  temperatures  may  lead  to  a  substantial 
elevation  in  notch  sensitivity.  However,  this  behavior  may  be  complicated  by 
the  development  of  alternate  damage  processes,  such  as  shear  bands. 

4.2  CMCs 

Studies  of  the  creep  and  rupture  of  CMCs  will  continue  with  emphasis 
on  materials  containing  creeping  fibers.  A  particular  emphasis  wlU  be  on 
matrix  cracking  that  arises  as  fiber  creep  relaxes  fiber  bridging  tractions 
(McMeeking.  E^vans).  The  experimental  studies  will  be  performed  on  SiC/SiC 
composites.  Hysteresis  loop  measurements  will  be  used  to  monitor  matrix 
damage  during  composite  creep,  using  procedures  devised  in  1993.  Models 
will  be  developed  based  on  time  dependent  fiber  bridging  concepts 
(McMeeking.  Cox). 

It  is  envisioned  that  the  lifetime  of  some  CMCs  will  be  dictated  by  time- 
dependent  rupture  of  the  fibers.  A  lifetime  prediction  tool  for  such  a 
composite  must  incorporate  the  knovdedge  of  fiber  strength  degradation  over 
time.  A  new  activity  will  be  initiated  to  address  this  problem  (Suo.  E^vans). 
The  initial  work  will  involve  a  surv^  of  data  in  the  existing  literature,  and  a 
comparison  with  available  models.  A  new  model  is  being  developed  for  single 
crystal  fibers.  This  model  involves  a  residual  pore  inside  a  fiber  which 
changes  shape,  under  stress,  via  surface  difiiision.  to  become  a  crack.  These 
issues  will  be  viewed  in  the  broad  context  of  fiber  and  composite 
manufacture. 
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5.  TRANSVERSE  PERPORBIANCE  OP  CMCs 

Analyses  and  tests  performed  in  1993  (Ashby.  Hutchinson.  Bao)  have 
highlighted  the  essential  Issues  related  to  components  that  experience 
combinations  of  transverse  tension  and  interlaminar  shear.  In  both 
loadings,  matrix  cracks  form  at  manufacturing  flaws  at  low  stresses,  of 
order  10-100  MPa.  These  cracks  extend  across  the  plies  and  interact 
minimally  with  the  fibers.  Although  the  crack  configurations  differ  for 
transverse  tension  and  interlaminar  shear  loadings,  multiple  cracks  always 
form.  This  multiplicity  of  cracking  causes  a  major  reduction  in  stiffness, 
which  can  cause  unacceptably  large  displacements  and  also  redistribute 
stress  into  other  areas.  The  formation  of  the  matrix  cracks  is  probabilistic  in 
nature  and  governed  the  size  distribution  of  manufacturing  flaws.  Design 
based  on  the  prevention  of  such  transverse  cracks  must  rely  on  weakest 
link  statistics,  usually  with  a  low  Weibull  modulus.  Alternative^,  it  may  be 
assumed  that  cracks  inevitably  form  and,  instead,  reliance  is  placed  on 
controlling  the  diminished  modulus  of  the  material,  after  matrix  cracking 
has  occurred.  This  approach  relies  on  having  3>D  architectures,  with 
transverse  fibers  introduced  locally  either  stitching  or  by  usixig  angle 
plies.  To  explore  this  possibility,  calculations  will  be  performed  (Hutchinson. 
Evans)  to  examine  fiber  architectures  that  lead  to  minimum  stlfihess  loss, 
subject  to  acceptable  in-plane  properties.  Based  on  these  calculations,  sub¬ 
elements  will  be  designed  that  test  out  the  concepts. 


6.  COMPRESSIVE  BEHAVIOR 

The  studies  completed  in  1993  on  the  compressive  failure  of  potymer 
matrix  composites  by  the  growth  of  kink  bands  (Budiansky,  Fleck)  will  be 
extended  to  metal  matrix  composites,  through  a  coordination  with  3M. 
Compressive  failure  of  A1  and  Ti  MMCs  with  small  diameter  fibers  has  been 
observed  by  3M  to  occur  in  accordance  with  the  same  kink  band 
mechanism  known  to  operate  in  PMCs  and  in  C/C  composites.  The  theoiy 
should  thus  extend  to  the  MMCs.  with  the  fiber  misalignment,  the  shear 
yield  strengtli  of  the  matrix  and  its  work  hardening  coefficient  as  the 
principal  variables.  A  comparison  between  the  theory  and  experimental 
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results  would  provide  the  basis  for  specifying  the  compressive  properties  of 
MMCs. 

Compression  failure  of  CMCs  occurs  different  mechanisms  (Ashby). 
The  dominant  failure  modes  are  similar  to  those  that  operate  in  porous 
brittle  solids  such  as  monolithic  ceramics,  concrete  and  rocks.  The  theory  is 
well  established  and  validated  for  these  materials,  ^plications  of  the  theory 
to  various  CMCs  will  be  made  and  applied  to  the  understanding  of  a 
behavior  of  pin-loaded  holes  (Evans.  Ashby). 


7.  THERMAL  PROPERTIES 

A  new  focus  on  the  thermal  properties  of  CMCs  and  MMCs  will  be 
initiated  in  1994.  Calculations  of  the  effects  of  matrix  cracks  in  the  thermal 
expansion  of  CMCs  will  be  made  (Hutchinson).  These  will  be  compared  with 
data  obtained  from  TMF  testing  (Zol^.  The  effects  of  such  cracks  on  the  in¬ 
plane  thermal  conductivity  will  also  be  calculated  (Hutchinson). 
Measurements  will  be  performed  using  the  laser  method  (Ashby). 

Thermal  conductivity  measurements  will  be  initiated  on  H  MMCs 
(Asht^).  These  will  be  used  to  understand  the  effects  of  the  fiber/matrix 
interphases  and  of  matrix  damage  on  the  transverse  and  in-plane  thermal 
conduction. 


8.  MATERIALS  SELECTION 

The  Cambridge  Materials  Selector  software  wiU  be  expanded  in  1994  to 
include  high  temperature  creep  design  with  the  corresponding  data  base 
(Ashby).  This  expanded  version  will  permit  estimates  to  be  made  of 
temperature  limits  for  MMCs  based  on  creep  controlled  TMF  and  on  the 
transverse  creep  of  components  with  unidirectional  reinforcements. 

0.  DESIGN  CALCULATIONS  AND  SUB-ELEMENT  TESTS 

A  larger  fraction  of  the  effort  in  1994  will  be  on  design  and  sub-element 
testing,  particularly  for  MMCs.  Discussions  are  now  in  progress  with  Pratt 
and  WhitnQr,  Textron  and  3M  to  perform  design  calculations  using  the 
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constitutive  equations  developed  at  UCSB  and  to  produce  sub-elements  for 
testing. 

The  design  emphasis  for  MMCs  will  be  on  various  diffusion  bonded 
Joints  with  Ti  matrices  and  monolithic  T1  attachments.  Two  specific 
subelements  are  envisaged.  The  first  involves  unidirectional^  reinforced 
rods  (or  plates),  clad  with  monolithic  metal.  The  purpose  of  the  cladding  is 
to  prevent  eaqposure  of  the  fibers  to  the  environment  and  to  mechanical 
abrasion.  The  design  of  clad  MMC  stn’.ctures  requires  consideration  of 
(i)  the  residual  stresses  resulting  from  thermal  mismatch  between  the 
cladding  and  the  composites  section,  (ii)  the  potential  for  fatigue  cracks  to 
initiate  and  grow  through  the  monolithic  material,  and  (iii)  the  interaction  of 
such  cracks  with  the  composite  section  and  their  influence  on  the  strength 
and  life  of  the  structure.  The  design  and  testing  of  such  subelements  (Zok. 
Leckie)  will  be  augmented  by  calculations  of  crack  growth  and  firacture, 
incorporating  the  effects  of  thermal  and  elastic  mismatch  between  the 
cladding  and  the  composite  (McMeeking).  The  dad  structures  will  also  be 
used  to  initiate  studies  on  the  reinforcement  of  holes  in  composite  sections 
with  monolithic  metal  patches,  as  drawn  in  Fig.  9.1  (Zok,  Suo).  The  second 
subdement  involves  the  attachment  of  a  MMC  actuator  rod  to  a  pin-loaded 
monolithic  section  (Fig.  9.2).  The  critical  design  issues  relate  to  the  strength 
and  fatigue  resistance  of  the  interfaces  between  the  composite  and 
monolithic  matrices.  Design  studies  shall  also  be  completed  on  rotor  rings 
with  special  efforts  made  to  produce  rule-based  design  procedures  which 
would  be  used  by  industry  at  the  conceptual  levd  of  design  to  determine 
sizes  and  the  effldent  disposition  of  material. 

For  CMCs,  the  sub-element  studies  would  be  based  on  the  calculations 
described  above  in  Section  5.  These  would  include  C  sections  and 
T  Junctions  (Fig.  9.3)  Negotiations  for  manufacturing  these  sub-elements 
will  be  initiated  and  tests  performed  at  UCSB. 

10.  AFFORDABLE  ICANUFACTURING 

As  our  understanding  of  composite  mechanics  and  its  interplay  with 
design  and  performance  has  evolved,  it  has  become  increasingly  evident  that 
cost  and  reproducibility,  are  major  constraints.  Even  as  processing 
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developments  make  the  prospect  of  affordable  hi^  temperature  fibers  more 
realistic,  evolving  knowledge  on  the  mechanical  and  thermochemical 
functions  of  interfaces  have  led  to  design  concepts  involving  carefully 
tailored  interphase  layers,  with  unfavorable  impact  on  cost.  Moreover,  if 
affordable  coated  fibers  were  available  today,  fabrication  costs  associated 
with  consolidation  and  pressure  densification  would  often  remain 
prohibitive.  Future  processing  and  manufacturing  activities  are  predicated 
on  these  issues,  especially  the  need  for  new  ideas,  and  the  related 
knowledge  base. 

10.1  MMCs 

Melt  processing  methods  provide  the  more  affordable  options  in 
composite  synthesis  with  the  added  benefit  of  near-net  shape  capability.  For 
continuous  fiber  composites  melt  Infiltration  also  enables  full  density  while 
minimizing  the  consolidation  stresses  that  typically  cause  premature 
reinforcement  failure  in  solid  state  processes.  However,  melt  processing 
requires  a  high  degree  of  thermochemical  ccmpatibility  between  matrix  and 
reinforcement  since  deleterious  diffiisional  interactions  would  be  accelerated 
by  the  liquid  phase.  Conventional  melt  processing  also  exhibits  limited 
ability  to  control  the  volume  fraction  and  spatial  uniformity  of  the 
reinforcements. 

Among  metal  matrices.  Tl  alleys  epitomize  unsuitability  for  direct  melt 
infiltration  owing  to  aggressive  reactivity.  Fiber  clustering  is  also  a  concern, 
even  in  solid  state  processes  based  on  powder  or  foil  matrices.  Composite 
consolidation  vapor  deposition  (PVD)  of  the  matrix  on  the  fibers  provides 
an  avenue  for  improving  homogeneity  of  fiber  spacing.  However,  present 
schemes  require  expensive  pressure  densification  with  its  many  problems.  A 
potential  solution  involves  a  hybrid  manufacturing  route  wherein  part  of  the 
matrix  is  first  applied  to  the  fibers  Ity  PVD.  The  pre-metallized  fibers  are 
then  assembled  into  a  preform  having  the  desired  shape  and  then  infiltrated 
with  the  remaining  matrix  in  liquid  form. 

Direct  infiltration  with  TL  alloys  could  be  feasible  owing  to  the  protection 
of  the  fiber  by  the  PVD  layer,  but  the  high  temperatures  involved  would 
exacerbate  the  diffuslonal  Interactions  at  the  fiber-matrix  interface.  An 
alternate  approach  involves  depositing  the  more  refiractoiy  constituents  of 
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the  matrix  (e.g..  H.  Nb.  Mo.  etc.)  by  FVD  and  then  infiltrating  with  the  lower 
melting  point  constituents  (e.g.  Al).  Based  on  stoichiometric  considerations, 
the  latter  approach  wotild  be  suitable  for  matrices  with  >  25  at.%  Al,  notably 
the  orthorhombic  and  02  alloys.  The  obvious  problem  with  this  approach  is 
the  homogenization  of  the  matrix  after  consolidation,  which  may  require 
lengthy  high  temperature  treatments  in  the  solid  state.  However,  a 
significant  part  of  the  matrix  synthesis  reaction  could  be  effected  in  the 
presence  of  molten  Al.  followed  a  final  heat  treatment  in  the  solid  state. 
While  this  lower  temperature  infiltration  approach  is  evidently  desirable 
firom  a  manufacturing  viewpoint,  it  is  not  clear  that  matrix  homogenization 
can  be  achieved. 

A  program  involving  modeling  and  e>q)erimental  work  will  be  initiated  in 
1994  to  generate  the  knowledge  base  appropriate  to  hybrid  approaches  for 
Tl  matrix  composites  (Levi,  Evans).  Cell  models  (sin^e  fiber  environment) 
would  be  developed  to  study  diffusional  interactions  and 
remelting/solidification  phenomena  as  a  function  of  processing  cycle 
(temperature-time  hlstoiy).  Ebq)eriments  would  be  performed  to  elucidate 
the  relevant  aspects  of  mlcrostructural  evolution  and  provide  the  reaction 
and  interdiffusion  kinetics  needed  to  calibrate  the  models.  Initial 
e3q)eriments  would  be  performed  by  Infiltrating  pure  Tl-wire  preforms  with 
molten  Al  and  subjecting  the  ‘‘composite*  to  dififerent  treatments  in  the 
semi-solid  state.  Subsequent  experiments  would  focus  on  developing  a 
metallization  route  for  Tl-Nb  alloys  on  SiC  fibers  and  on  the  relevant 
interactions  with  infiltrated  Al.  Larger  scale  modeling  issues  would  be 
tackled  in  1995  if  the  proposed  approach  appears  promising. 

Ongoing  activities  on  the  imderstanding  of  microstructure  evolution 
and  its  relationship  to  properties  in  in-situ  TMC  ^sterns  based  on  TiB 
reinforcements  would  be  continued  (Levi).  These  are  by  nature  affordable 
composites  which  ejdiibit  inherent  thermochemical  stability  and  may  be  cast 
into  shapes  using  conventional  Tl  processing  techniques.  A  potential 
application  of  these  materials  would  be  in  joints  with  unidirectionally 
reinforced  composites,  wherein  their  higher  modulus  and  creep  resistance 
combined  with  acceptable  toughness  and  isotropic  properties  could  be 
advantageous.  It  is  also  anticipated  that  these  materials  could  be  used  for 
cladding  in  PVD  or  plasma-sprayed  form,  thereby  reducing  the  potential  for 
fatigue  crack  initiation  in  the  cladding.  Since  TIB  is  thermochemically  stable 
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with  practically  all  Ti  matrices  of  interest  in  fiber  composites,  such 
strengthening  concepts  may  be  readily  implemented. 

lo.a  CMCs 

Measurements  and  observations  in  1993  have  shown  that  strong,  high 
strain  to  failure  CMCs  can  be  fabricated  using  an  inexpensive  method  that 
involves  a)  packing  a  powder  around  fibers  within  a  fiber  preform  using 
pressure  filtration  and  b)  making  the  powder  matrix  strong  by  heat 
treatment  followed  hy  infiltration  with  a  liquid  precursor  that  decomposes  to 
an  inorganic  material.  A  composite  made  this  way.  with  potycxystalline 
alumina  fibers  in  a  silicon  nitride  matrix,  demonstrated  that  the  matrix 
deflects  the  crack.  This  observation  is  significant  since  it  suggest  that  a 
class  of  CMCs  can  be  processed  without  needing  weak  fiber/matrix 
interlaces.  The  potential  of  this  observation  will  be  colored  (Lange.  E^rans), 
processing  a  composite  with  strong,  pofyciystalline  alumina  fibers  in  a 
mulUte  matrix  because  the  thermomechanical  properties  of  muUite  minimize 
thermal  stresses  and  resist  creep.  In  addition,  the  thermal  expansion 
mismatch  is  relative^  small.  Mixed  Al.  Si  metal  alkoxlde  precursors  which 
can  be  gelled  in^situ,  prior  to  decomposition,  will  be  used  to  strengthen  the 
matrix. 

Manufacturing  studies  would  initiate  with  understanding  the  precursor 
infiltration  into  mullite  power  compacts.  The  densification  of  the  matrix 
would  be  determined  as  a  function  of  the  cyclic  infiltration.  Microstructure 
changes  would  be  controlled  to  avoid  flaw  populations  during  densification. 
The  fracture  toughness  and  the  strength  of  the  matrix  woiild  be  determined 
as  a  function  of  the  number  of  precursor  Infiltration  cycles.  Composite 
processing  would  initiate  with  precursor  infiltration  into  alumina  fiber 
preforms  by  pressure  filtration,  with  emphasis  on  the  colloidal  aspects  of 
this  processing  step.  The  goal  would  be  to  determine  the  processing 
conditions  needed  to  produce  a  matrix  that  optimizes  the  ability  to  deflect 
cracks  without  degrading  fiber  strength.  To  optimize  cx>mposite  processing, 
panels  for  testing  tmder  conditions  of  both  strain  and  stress  control  would 
be  manufactured. 
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11.  STRESS  AND  DAMAGE  SENSORS 

The  extensive  exploitation  of  the  optical  fluorescence  method  of 
measuring  stresses  in  sapphire  fiber  and  alumina-containing  ceramic 
composites  begun  in  1993  will  be  continued  in  1994  (Clarke,  WadlQr).  The 
emphasis  is  on  using  the  method  to  understand  basic,  unresolved  issues  in 
stress  redistribution  in  composites  by  the  direct  measurement,  with  high 
spatial  resolution,  of  the  stresses  themselves.  Particular  attention  will  be 
paid  to  determining  the  stress  distribution  associated  with  interfacial 
sliding.  One  of  the  problems  to  be  addressed  relates  to  new  concepts  for 
oxidation  resistant  interfaces  within  MMCs  and  CMCs,  particularly  the 
concomitant  roles  of  fiber  roughness  and  sintering  on  interface  sliding  and 
debonding,  after  exposure  to  high  temperatures  and  cyclic  loadings.  For  this 
purpose,  fibers  with  fugitive,  low  modulus  coatings  will  be  explored  and 
fluorescence  measurements  used  to  understana  stress  evolution  and  its 
connection  vlth  fiber  durability  within  the  composite.  A  second  problem 
relates  to  the  distinction  between  the  line  spring  and  large  scale  sliding 
models  for  fiber  bridging  (Budlansl^.  Hutchinson),  so  as  to  determine  the 
range  of  applicability  of  the  two  models.  The  two  competing  models  predict 
different  distributions  of  stresses  in  the  fibers  within  the  bridging  zone  and 
hence  are  amenable  to  validation  on  the  basis  of  the  measured  stress 
distribution. 

Two  approaches  to  measuring  local  damage  are  under  development  and 
will  be  the  focus  of  the  sensor  activities.  One  is  the  use  of  acoustic  methods 
(Wadley)  to  probe  local  variations  in  the  elastic  modulus  of  CMCs  as  a 
function  of  load.  This  should  provide  a  means  of  mapping  the  distribution  of 
damage  which  can  be  compared  directly  with  the  predictions  of  continuum 
damage  mechanics  models.  The  second  approach  (Clarke)  is  to  detect  the 
third  harmonic  signal  generated  by  the  presence  of  local  damage. 
Preliminary  experimental  results  obtained  in  1993  concerned  w<th  the 
detection  of  crack-like  voids  in  thin  metal  lines,  together  with  computer 
simulation  studies,  have  demonstrated  the  viability  of  the  technique.  This 
work  will  be  extended  in  order  to  detect  damage  accumulation  in  CMCs  and 
MMCs. 
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Abstract:  A  method  to  non- intrusively  measure  stresses  both  at  and 
below  the  surface  of  structural  components  is  developed.  Based  on  the 
piezo-spectroscopic  effect  in  optical  fluorescence,  this  technique 
affords  high-accuracy  stress  measurement  with  spatial  resolution  on  the 
order  of  one  micron.  The  th<'oretical  basis  for  the  piezo-spectroscopic 
effect  is  reviewed,  and  se  •  ral  concepts  relevant  to  its  application  are 
established.  The  applicability  of  fluorescence  spectroscopy  to  stress 
measurement  is  demonstrate’  on  two  sapphire-fiber  reinforced  composites, 
one  with  y-TiAl  and  the  ot.^er  with  a-Ti  as  the  matrix.  The  residual 
stress  distributions  along  the  fibers  are  measured  using  a  through.- 
focusing  procedure,  and  a  mean  volume  stress  is  determined  using  a 
transmission  fluorescence  configuration.  The  effective  temperature 
below  which  residual  stresses  can  no  longer  be  relieved  by  creep  and 
yielding  is  determined  using  an  elastic  model.  In  the  TiAl -matrix 
composite,  this  temperature  is  approximately  400°C  below  the  processing 
temperature,  indicating  extensive  stress  relaxation  during  cooling.  The 
elastic  solution  fails  tp  provide  a  reasonable  estimate  of  the  stress- 
free  temperature  for  the  Ti-matrix  system.  This  is  attributed  to 
extensive  matrix  yielding,  and  a  complete  treatment  will  require 
incorporation  of  inelastic  constitutive  behavior  into  the  model. 


Xeywozde:  stress  measurement,  residual  stress,  piezo-spectroscopy, 
fluorescence,  sapphire,  ruby,  TiAl  matrix,  Ti  matrix 


1 .  Introduction 

The  need  to  study  microscale  properties  of  materials  has  led  to 
the  development  of  various  high-resolution  microscopic  techniques,  such 
as  scanning  electron,  transmiission  electron,  scanning-tunneling,  and 
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atomic  force  microscopy.  However,  advances  in  morphological  and 
crystallographic  evaluation  have  not  been  complemented  with  a 
corresponding  development  in  techniques  for  measuring  micromechanical 
properties.  Although  such  techniques  as  microfocus  X-ray  diffraction 
[IJ  and  high-accuracy  strain-field  mapping  (HASMAP)  [2]  have  made 
possible  the  evaluation  of  surface  strains  with  iit^roved  spatial 
resolution,  these  methods  are  inherently  limited  to  superficial 
measurements.  In  the  present  study,  a  recently-developed  technique 
based  on  the  piezo-spectroscopic  effect  in  optical  fluorescence  is 
en^loyed  to  obtain  stress  information  well  below  the  sanple  surface 
[3,4,5].  This  technique  has  been  shown  to  routinely  combine  micron- 
scale  spatial  resolution  with  a  corresponding  stress  resolution  to 
approximately  10  MPa. 

One  of  the  goals  of  the  spectroscopic  stress  measurement  technique 
is  to  provide  the  means  for  experimentally  determining  spatial  stress 
distributions  in  residually-stressed  or  externally-loaded  components. 
Such  data  is  not  only  crucial  in  validating  finite  element  predictions, 
but  can  be  integrated  into  on-line  quality  control  and  evaluation 
procedures  requiring  the  precise  knowledge  of  stress  levels  within 
structural  components.  In  the  present  study,  we  determine  the  residual 
stress  distribution  along  sapphire  fibers  in  two  composites,  one  with  a 
y-TiAl  matrix  and  the  other  an  a-Ti  matrix.  The  stresses  develop  upon 
cooling  from  the  fabrication  ten^rature  as  a  result  of  the  difference 
in  thermoelastic  properties  between  matrix  and  reinforcement.  These 
stresses  are  only  partially  relieved  by  creep  at  high  temperatures  and 
dislocation  plasticity  at  low  temperatures.  Because  of  these  relaxation 
processes,  particularly  near  the  fiber  ends,  the  stress  distribution 
along  the  fibers  is  not  amenable  to  con^utational  analysis  and  has  to  be 
measured  directly.  To  determine  the  stress  distribution,  we  have 
enployed  the  fluorescence  spectroscopy  technique  and  applied  it  in  two 
configurations:  one,  in  which  successively  deeper  regions  are  probed  by 
focusing  an  optical  lens  down  the  axis  of  the  reinforcing  fibers,  and 
the  second,  in  which  the  whole  fiber  is  probed  simultaneously  in 
transmission.  The  two  procedures  are  detailed  following  a  description 
of  the  physical  basis  of  the  piezo-spectroscopic  technicpie. 


2 . 1  Pi«so-8p«ctsoseoplc  Effect 

The  piezo-spectroscopic  (PS)  effect  is  based  on  the  sensitivity  of 
electronic  transitions  of  substitutional  ions,  or  phosphors,  to  stress- 
induced  changes  in  the  local  crystal  field  of  the  host  lattice.  A 
common  example  is  ruby,  in  which  Cr^*  ions  substitute  for  Al^'^  ions  on 
the  corundum  lattice.  When  ruby  is  excited  with  ultraviolet  or 
sufficiently  energetic  visible-wavelength  radiation,  the  Cr^*  ions 
undergo  an  electronic  transition  from  the  ^A2  ground  state  to  the  broad 
^Ti  and  ^T;  bands,  followed  by  a  rapid  non-radiative  decay  to  one  of  two 
components  of  the  level  (Fig.l).  The  subsequent  decay  from  to  the 
^A;  ground  state  is  accompanied  by  the  radiation  of  photons  having 
frequencies  within  the  visible  spectrum,  a  phenomenon  called  optical 
fluorescence.  In  a  stress-free  single  cr^tal  of  sapphire,  the 
fluorescence  lines  corresponding  to  the  2A— ^^Aj  transition  (R2)  and  the 
E-4^A2  transition  (Ri)  are  separated  by  approximately  29  cm"^  and  occur 
at  frequencies  of  14448  and  14419  cm~^,  respectively  [6] .  However,  as 
the  crystal  field  about  the  Cr^*  ion  is  perturbed,  its  electronic  states 
are  displaced  with  respect  to  one  another,  resulting  in  a  systematic 
frequency  shift  of  the  fluorescence  signal.  Although  the  crystal  field 
can  be  affected  by  external  electric  fields  or  in  the  vicinity  of 


Fig.l  Electronic  transitions  for  a  Cr^'^  ion  in  a  sapphire 

host  lattice,  illustrating  the  origin  of  the  Ri  and  R2 
fluorescence  lines. 


lattice  defects,  the  present  eoqphasis  is  on  the  fluorescence  frequency 
shifts  caused  by  tractions  iiqposed  upon  the  host  lattice  —  the  piezo- 
spectroscopic  effect. 


2.2  i^ppllaetlon  ot  the  98  Effect  to  Streee  Meaeacoaeat 

As  mentioned  above,  the  piezo-spectroscopic  effect  can  be  en^loyed 
to  relate  fluorescence  freqpiency  shifts  to  the  local  stress  tensor. 

This  constitutive  relationship  is  usually  described  through  the  piezo- 
spectroscopic  tensor,  FI : 


AvtCTij)  =  nijOij  (1) 

where  Av  refers  to  the  frequency  shift,  IIij  is  the  ij-th  conponent  of 
the  piezo-spectroscopic  tensor,  and  O^j  is  the  corresponding  coitponent 
the  stress  tensor  (in  the  crystallographic  basis  of  the  host 
lattice) .  The  higher-order  terms  in  the  expansion  of  the  frequency 
shift  have  been  confirmed  to  be  negligible  for  most  purposes  [7,8]  and 
are  therefore  excluded  from  Eqn. (1) . 

In  addition  to  stress  dependence,  the  fluorescence  frequency  has 
varying  degrees  of  sensitivity  to  ten^rature,  applied  electric  field, 
and  dopant  concentration.  To  a  first-order  approxiination,  the  various 
contributions  to  the  fluorescence  line  shift  are  taken  to  be  linearly 
superposable  [9],  such  that  the  net  frequency  shift  can  be  described  as: 

Av  =  Av(Oij)  +  Av(Ei)  +  Av(T)  +  Av(c),  (2) 

where  the  individual  consonants  correspond  to  the  stress,  applied 


•I«ctrlc  tcmpcratura,  and  concantration  dapandanca,  raspactively . 

For  our  Invaatigations,  no  axtamal  alactric  fields  are  applied 
and  it  is  assunad  that  tha  dopants  ara  homoganaously  distributed  on  the 
host  lattlcsr  so  that  Eqn.  (2)  reduces  to  tha  stress-  and  ten^rature- 
dapandant  tarms.  Tha  taiig>aratura  dapandanca  of  frequency  can  be 
described  through  a  scalar  correction  factor: 

AV(T)  =  P(T)T  (3) 

where  P  is  experimentally  found  to  have  a  weak  dependence  on 
tesq>aratura  over  small  excursions,  and  varies  slightly  between  the  Ri 
and  R2  lines  [10-12] . 

Tha  piazo-spactroscopic  tensor  is  a  physically-measurable  material 
parameter  that  reflects  the  symmetry  of  the  host  lattice.  For  the  case 
of  sapphire  as  the  host,  the  site  synmetry  of  Cr^'^  (Pm)  reduces  the 
piezo-spectroscopic  tensor  to  the  form: 

n,  0  o' 

n  =  0  n,  0  (4) 

0  0 

where  11,  and  II,.  refer  to  the  cong>onents  oriented  along  the  a-  and  c- 
axes  of  the  sapphire  unit  cell,  respectively.  By  uniaxially  stressing 
single  crystals  of  ruby.  He  and  Clarice  found  n,  and  11,  to  be  3.45  and 
1.51  cm~^/GPa  for  the  Ri  line  and  2.75  and  2.11  cm*^/GPa  for  the  Rj  line, 
respectively  [7] .  For  the  case  of  a  reinforcing  sapphire  fiber  having 
its  axis  oriented  along  the  crystallographic  c-axis  of  sapphire,  the 
frequency  shift  is  obtained  directly  from  Eqns. (1)  and  (4): 

AvCdij)  *  n,<o„  +  o„>  +  n,o„,  (5) 

where  ,  Om,  and  o„  are  the  radial,  tangential,  and  axial  components 
of  the  stress  tensor,  respectively. 

Although  Eqn. (5)  is  adequate  for  providing  a  semi -quantitative 
measure  of  the  stress  state  within  the  reinforcing  fiber,  it  is 
insufficient  when  a  full  stress  tensor  solution  is  required.  The 
problem  ariseB  due  to  the  lack  of  a  universal  relationship  between  the 
individual  cong>onents  of  the  stress  tensor.  In  most  situations,  such 
relationships  can  only  be  determined  from  numerical  (e.g.,  finite 
element)  calculations.  However,  there  exists  a  limited  number  of 
sin^lified  geometries  and  constitutive  models  for  which  the  relationship 
between  stress  tensor  consonants  can  be  found  analytically.  One  such 
configuration  is  a  cylindrical  fiber  embedded  in  an  infinite  matrix 
having  dissimilar  but  isotropic  thermoelastic  properties.  With  the 
assumptions  that  both  conponents  are  linear-elastic  and  that  no 
displacement  is  allowed  at  the  fiber/matrix  interface  (fixed  boundary 
conditions),  the  conponents  of  the  residual  stress  tensor  can  be 
rigorously  related  using  elasticity  theory: 

®*e  =  ®rr 

^  (6) 

®rr 


where: 


and 


9teJ**AadT 

E;“AodT 


9t 
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and  ara  the  fiber  and  matrix  elastic  moduli  and  Poisson's 
ratios,  respectively,  and  are  allowed  to  have  a  teiq}erature  dependence 
from  room  temperature,  To,  to  the  effective  stress-free  temperature. 

Combining  Eqns. (5)  and  (6),  the  frequency  shift  measured  within 
the  bulk  of  a  fiber  composite  can  be  uniquely  related  to  the  thermal 
residual  stress  state. 


3.1  Threugh-roctts  Stress  Messttrement 

In  the  through-focus  technique,  the  spatially-resolved  stress 
distribution  is  obtained  by  incrementally  focusing  a  laser  beam  along 
the  fiber  axis  and  collecting  the  fluorescence  that  is  excited  about  the 
focal  point  (Fig. 2) .  Both  the  excitation  and  collection  of  the 
fluorescence  signal  are  inplemented  using  a  conventional  optical 
microscope;  as  a  consequence,  the  spatial  resolution  of  the  technique  is 
limited  only  by  the  aberrations  of  the  optical  conponents,  and 
ultimately  by  the  Rayleigh  resolution  limit. 

Although  the  through-focus  measurement  is  nominally 
St rai^t forward,  the  optical  lenses  that  are  used  for  the  excitation  and 
collection  of  the  fluorescence  signal  have  a  finite  depth  of  field. 
Consequently,  the  signal  collected  from  any  depth,  below  the  surface 
is  a  spatial  convolution  of  the  true  signal  distribution,  f(s),  with  a 
system  response  function,  r(z),  which  acts  to  spread  out  the  collected 
information  over  a  finite  volume  in  space  [13] .  Typically,  the 
collection  volume  has  dimensions  of  approximately  30  microns  along  the 
z-axis  (Fig. 2)  and  S  microns  in  the  x-y  plane,  although  the  exact  values 
depend  strongly  upon  the  collection  optics,  including  lenses,  slits,  and 
apertures.  As  a  result  of  the  finite  collection  volume,  the  measured 
frequency  distribution  may  need  to  be  numerically  deconvoluted  in  order 
to  determine  the  true  spatial  frequency  distribution.  Although  the  data 
in  the  following  experiments  were  not  deconvoluted,  such  a  procedure 
becomes  increasingly  important  in  cases  where  large  stress  gradients 
exist  within  the  sanple. 

The  collected  signal,  c(z),  from  position  z  along  the  fiber,  can 
be  described  mathematically  as  a  convolution  of  f(z)  and  r(z) : 

c(z)  »  f(z)  ®  r(z)  *  J*** f(^r(z  -  §)d^,  (7) 

where  the  integral  describes  the  passing  of  a  function  f(z)  through 
r(z).  The  response  function  can  be  experimentally  measured  if  the 
integrated  intensity  under  the  Ri  and  R2  fluorescence  peaks  is  taken  to 
be  the  system  fiuiction,  c(z).  Assuming  that  the  integrated  fluorescence 
intensity  is  only  dependent  on  the  geometry  of  the  sanple  (e.g.,  the 
intensity  is  zero  outside  the  sanple  and  has  a  constant  value.  A,  within 


*  If  the  tenperature  dependence  of  Vt.m  Qf.m  ignored. 


Eqn.  (6)  degenerates  to:  ^ 
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Fig. 2  Through-focus  fluorescence  configuration  illustrating 

that  the  signal  is  collected  from  a  finite  volume.  The 
fiber  axis  is  oriented  parallel  to  the  sapphire  c-axis. 


the  sample),  Eqn. (7)  can  be  solved  explicitly  for  the  response,  r(z): 


r<z)  = 


1  dc(z) 
A  dz 


(8) 


Once  the  response  function  is  knoim,  the  true  function  can  be  recovered 
through  the  application  of  the  convolution  theorem  of  Fourier 
transforms: 


(9) 


Ultimately,  the  experimental  measurement  of  the  stress 
distribution  using  the  through-focus  technique  consists  of  a  two-step 
process:  collection  of  the  frequency  data  at  incremental  depths  below 
the  cosq>onent  surface  and,  if  needed,  numerical  deconvolution  of  the 
frequency  profile  using  a  pre-calibrated  response  function.  However, 
the  technique  becomes  increasingly  cumbersome  to  apply  as  the  component 
size  increases,  motivating  the  development  of  the  transmission 
fluorescence  technique. 


3.2  Vsuaalaaloa  Straas  Maaatt; 


lat 


Although  tha  through-focua  tachniqua  providaa  a  powarful  tool  for 
aicroseala  straas  profiling,  it  suffars  from  savaral  practical 
liaitations.  First,  tha  tachniqua  ralias  on  tha  collaction  of  a  large 
data  sat.  Zn  ganaral,  tha  spatial  resolution  scales  in  direct  proportion 
to  tha  nunbar  of  data  points  collected.  Thus,  whan  tha  region  of 
interest  is  large,  tha  time  required  to  collect  a  sufficient  amount  of 
data  to  provide  tha  prescribed  spatial  resolution  may  be  prohibitive. 

The  second  drawback  derives  from  the  physical  limitations  imposed  by  the 
probing  optics  and  sample  geometry.  When  the  focal  plane  reaches  a 
critical  depth  below  the  surface  of  the  sample  (usually  related  to  some 
characteristic  sangile  dimension,  such  as  grain  size  or  fiber  diameter) , 
the  system  response  function  broadens  dramatically,  resulting  in  a  rapid 
loss  of  resolution.  Finally,  a  thorough  mapping  of  the  stress 
distribution  sMy  require  an  involved  numerical  deconvolution  procedure 
to  correct  for  the  systui  response.  All  of  the  above  limitations  can  be 
circumvented  by  resorting  to  the  transmission  fluorescence  technique, 
wherein  the  whole  sample  volume  is  simultaneously  excited  into 
fluorescence  by  an  external  illumination  source. 


to  collection  optics 
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Illumination 
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Fig. 3  Transmission  fluorescence  configurations  for  a  fiber 
embedded  in  an  opaque  matrix:  (a)  total  internal 
reflection  and  (b)  complete  absorption  at  the  fiber/ 
matrix  interface. 


In  the  transmission  fluorescence  technique,  the  excitation  source 
is  placed  at  one  end  of  the  san^le  and  the  fluorescence  signal  is 
collected  from  the  other  end  using  a  low-magnification  objective  lens 
(Fig. 3).  This  configuration  is  not  only  simple  to  implement,  but  has 
the  additional  advantage  over  the  through-focusing  technique  of  allowing 
components  of  Indefinite  length  to  be  probed. 

In  general,  the  transmitted  frequency,  v,  at  which  the 
fluorescence  line  is  located  can  be  described  using  a  geometrical 
signal -weighting  function,  w: 


I  w  V  dv 

V  *  — - ,  (10) 

wdV 

Jv 

where  v  is  the  frequency  of  the  fluocescence  line  within  an 
infinitesimal  fiber  volume  element,  dV.  In  analyzing  stresses  within 
embedded  fibers  using  the  transmission  fluorescence  geometry,  two 
extreme  conditions  are  possible:  total  internal  reflection  and  conplete 
absorption  at  the  fiber/matrix  interface  (Fig. 3).  Neglecting  the 
effects  of  absorption  within  the  bullc  of  the  fiber,  it  can  be  shown  that 
for  the  case  of  complete  reflection  at  the  fiber/matrix  interface,  the 
weighting  function  is  independent  of  position  and  the  transmitted  signal 
is  simply  the  volume-averaged  frequency  over  the  illuminated  region  of 
the  sanple: 


V  =  ~ j^vdV.  (reflection)  (11) 

In  contrast  to  the  above  situation,  for  the  case  of  complete 
absorption  at  the  fiber  interface,  the  weighting  function  is  sensitive 
to  the  numerical  aperture  of  the  collection  optics.  The  full  treatment 
of  this  problem  is  too  extensive  to  provide  physical  insight;  however, 
with  a  few  sioplifying  assumptions,  a  useful  approximation  for  the 
transmitted  frequency  can  be  made.  It  is  assumed  that,  for  the  case  of 
a  cylindrical  fiber  embedded  in  a  continuous  matrix,  there  is  no  radial 
or  angular  dependence  of  frequency  (or  equivalently,  of  stress) .  This 
assusqption  is  expected  to  hold  well  within  the  bullc  of  the  sample 
(Eqn.(6)),  where  the  fiber  can  be  treated  as  an  equivalent  inclusion 
using  the  Eshelby  formulation  [14].  (The  premise  that  the  fiber  stress 
state  is  constant  is  Icnown  to  fall  near  the  fiber  ends,  where  surface 
relaxation  distorts  the  stress  field.)  Additionally,  the  sample 
chiclcness,  t,  is  assumed  to  be  simultaneously  greater  than  the  fiber 
diameter  and  less  than  a  characteristic  length,  to,  defined  by  the 
system  geometry: 


to  s  R/Sz  ,  (12) 

where  R  is  the  filler  radius  and  Oj  is  the  angle  the  refracted  beam  ma)ces 
with  the  fiber  axis.  Oj  can,  in  turn,  be  related  to  the  numerical 
aperture,  NA,  of  the  collection  lens  and  the  refractive  index,  n;it>erf 
the  fiber:  62*sin~^  (NA/nf^oar)  *  Under  the  former  assumptions,  the 
transmitted  frequency  can  be  calculated: 
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V  =  -Aii 
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X(*)vd2 


-r 


X{z)dz 
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(13) 


It  is  interesting  to  note  that  even  for  the  case  of  cosoplete  absorption 
of  the  signal  at  the  fiber/matrix  interface,  the  weighting  function  has 
only  a  wea]c  polynomial  dependence  on  depth.  As  a  result,  for  very  short 
fibers  (t*to),  Eqns. (11)  and  (13)  predict  very  similar  frequencies  in 
transmission,  ma)cing  the  determination  of  the  exact  extent  of  reflection 
at  the  interface  unnecessary.  However,  for  long  fibers,  the  assunptions 
behind  Eqn. (13)  fail,  and  it  can  be  easily  shown  that  reflection  at  the 
interface  is  crucial  to  obtaining  a  signal  in  transmission. 


4.  Bxpeslaental  Procedure 

The  characteristics  of  the  two  composites  used  in  the  evaluation 
of  fiber  stress  distributions  are  summarized  in  Table  1.  The  samples 
were  machined  to  orient  the  free  surfaces  normal  to  the  fiber,  axes,  and 
the  exposed  ends  were  metallographically  polished  to  a  one-micron  finish 
to  remove  machining  stresses  and  attain  an  optically  smooth  surface. 

The  details  of  the  through-focus  and  transmission  fluorescence 
experiments  are  described  in  Sections  3.1  and  3.2,  respectively.  The 
general  data  collection  and  analysis  procedures  for  both  configurations 
were  similar.  The  fluorescence  spectra  were  collected  using  a  triple- 
grating  Raman  spectrometer.'^’  The  line  frequencies,  widths,  and 
intensities  of  the  fluorescence  signals  were  obtained  by  fitting  the 
spectra  to  a  doiible  pseudo-Voigtian  function,  one  for  each  of  the  R- 
lines,  using  a  commercial  curve-fitting  software  pac)cage.*^  With  each 
fluorescence  spectrum,  a  characteristic  neon  line  was  recorded  as  a 
spectral  reference  to  correct  for  possible  mechanical  drift  in  the 
spectrometer.  The  saiif>le  temperature  was  monitored  to  correct  for  the 
tenqperature  dependence  of  the  fluorescence  shift  (Eqn. (3)). 

The  two  conposite  sanples  exhibited  a  relatively  wide  range  of 
thermoelastic  and  mechanical  properties.  The  titanium  aluminide  was  not 
expected  to  undergo  extensive  plastic  deformation  far  below  the 
processing  temperature,  whereas  the  titanium  matrix  was  considered 
susceptible  to  plastic  yielding  under  the  thermal  residual  stresses.  In 


TABLE  1 — Conposite  systems  used  for  residual  stress  measurements 


Sample  ribex 

Coatina 

Matrix 

Consolidation 

IMC  sapphire 

NMC  rubv 

M0/AI2O3 

none 

Y-TiAl 

O-Ti 

HIP*  at  1070'’C 

UHP*  at  900®C 

*HIP  “  hot  isostatic 

press 

*UHP  - 

uniaxial  hot  press 

Raman  Microprobe  Model  T64000,  Instruments  SA. 
Spectra  Calc,  Galactic  Industries  Corporation. 


addition,  to  ensure  that  the  transmission  configuration  could  be  applied 
to  arbitrary  specimen  sizes,  both  a  thin  and  a  thiclc  sample  were 
examined;  the  titanium-matrix  cosf>osite  was  approximately  15  millimeters 
thick  (corresponding  to  a  fiber  aspect  ratio  of  70:1),  compared  to  0.7 
millimeters  for  the  aluminide  sample  (aspect  ratio  of  5.5:1). 


5.  Keeulte  and  Discussion 

Figures  4a  and  b  show  the  frequency  distributions  along  the  fiber 
axes  in  the  TiAl-  and  Ti-matrix  composites,  respectively.  Marked 
differences  in  the  residual  stress  distributions  are  immediately 
apparent  between  the  two  conposites.  For  the  case  of  the  interroetallic- 
matrix  conposite,  the  stress  state  (related  to  frequency  through 
Eqns. (5)  and  (6))  builds  up  rapidly  from  the  surface-relaxed  to  the  bulk 
state.  The  stress  reaches  a  plateau  value  within  approximately  one 
fiber  diameter,  indicating  conplete  load  transfer  via  the  interface  and 
full  plane-strain  constraint. 

In  contrast  to  the  TiAl  matrix,  the  more  compliant  titanium  matrix 
exhibits  much  longer  load  transfer  lengths,  as  evidenced  by  the 
steadily-increasing  stress  at  depths  up  to  three  times  the  fiber 
diameter  (Fig. 4b).  In  fact,  the  depth  corresponding  to  conplete 
constraint  was  beyond  the  maximum  possible  probing  depth  (approximately 
600  microns,  or  three  fiber  diameters) .  To  measure  the  transfer  length, 
it  was  necessary  to  rely  on  transmission  fluorescence  data.  Once  the 
bulk  stress  state  was  measured  in  transmission,  the  load  transfer  length 
was  estimated  from  the  slope  of  the  frequency  curve  of  Fig. 4b  to  )3e  nine 
times  the  fiber  diameter. 

The  frequency  corresponding  to  the  signal  obtained  in  transmission 
is  superimposed  onto  the  vertical  axis  in  Figs. 4.  In  order  to  verify 
the  validity  of  the  transmitted  frequency,  signals  were  collected  from 


z/R 

z/R 

(a) 

(b) 

Fig. 4  Ri  frequency  shift  distribution  along  the  reinforcing  sapphire 
fiber  in  (a)  TiAl -matrix  and  (b)  Ti-matrix  composites.  The 
transmitted  frequency  is  indicated  by  the  shaded  box  in  the 
upper  left  c'rner.  The  total  fiber  lengths  for  the  TiAl  and 
Ti  coirposites  were  0.7  and  15  mm,  respectively. 


TABLE  2 — Thazinoelastic  properties  of  TiAl-  and  Ti-  matrix  con^osites 


Seaple 

■#(©»•) 

B«(OPa) 

vt 

v« 

Aorxio"*  "c"*) 

Ti/AljOj 

380 

110 

0.26 

0.32 

2.35 

TiAl/AliOj 

380 

160 

0.26 

0.24 

3.55 

^Values  for  B,  v,  and  Aa  are  room  temperature  quantities.  Aaea„  -Q;. 


both  ends  of  the  specimens  to  assure  invariance.  For  the  thick  Ti/Ruby 
sample,  the  signa.^  was  collected  from  several  focal  planes.  Invariance 
of  these  signals  guarantees  that  extensive  reflection  at  the  interface 
obtains,  in  accordance  with  Eqn. (11) .  In  all  cases,  the  transmitted 
signals  were  consistent  to  within  experimental  accuracy.  As  expected, 
the  frequency  shift  determined  from  transmission  falls  slightly  below 
the  plateau  value  for  the  TiAl  conposite  (Fig. 4a).  For  the  Ti -matrix 
composite,  the  transmitted  frequency  is  almost  twice  the  maximum  value 
measured  using  the  through-focus  technique  (Fig. 4b) .  Extrapolation  of 
the  spatially-resolved  data  to  the  transmitted  value  reveals  that  the 
bulk  stress  state  is  not  reached  until  9  fiber  diameters  below  the 
surface,  as  mentioned  above. 

The  residual  stress  conponents  within  the  bulk  of  each  fiber  were 
calculated  from  the  plateau  (i.e.,  fully-constrained)  frequency  shift 
using  Eqns. (5)  and  (6)  in  conjunction  with  the  appropriate  thezmoelastic 
constants  (Table  2)  and  their  known  temperature  dependencies.  The 
results  are  presented  in  Table  3.  Using  these  values  in  the  elasticity 
solution  for  the  thermal  residual  stress, the  effective  elastic  cooling 
range,  AT,;;,  was  calculated  for  the  two  composites  (Table  3) .  This 
range  corresponds  to  the  effective  freezing,  or  stress-free, 
temperature,  below  which  creep  and  dislocation  plasticity  mechanisms 
become  ineffective  (i.e.,  material  is  essentially  linear-elastic)  snd 
above  which  the  residual  stresses  are  fully  relaxed.  For  the  TiAl- 
matrix  composite,  the  effective  stress-free  temperature  is  found  to  be 
approximately  670°C  (400°C  below  the  peak  processing  temperature), 
indicating  that  extensive  relaxation  had  occurred  via  high-temperature 
creep,  as  well  as  possible  microcracking  and  interfacial  slippage  at 
lower  temperatures . 

When  the  residual  stresses  were  used  to  cal>  ulate  T,{f  for  the 
Ti/r\iby  composite,  the  resulting  temperature  was  found  to  exceed  the 
melting  temperature  of  titanium.  This  apparent  discrepancy  suggests 
that  the  degree  of  yielding  in  this  composite  is  so  extensive  that  the 


TABLE  3 — Residual  stresses  and  effective  freezing  temperatures 
of  TiAl-  and  Ti-  matrix  composites  based  on  linear- 
elastic  model 


Sample 

aii(MPa) 

arr(MPa) 

A'i;,,(°C) 

Ti/AljOa 

-1325 

-285 

N/A 

TiAl/Al,0-, 

-955 

-290 

650 

Solutions,  derived  using  boundary  conditions  set  forth  in  section  2.2, 


are:  a. 


=  Ej“AadT  and  0„  =  (91  +  l)E5"AadT.  AT^ff  h  Tg,,  - 


2800 


Fig. 5  Graphical  solution  of  the  piezo-spectroscopic  equation  for 
Ti/ruby.  The  stress  states  corresponding  to  the  elastic 
model,  as  well  as  to  the  extreme  case  of  con^lete  radial 
relaxation  are  indicated.  Av-3.95  cm*^,  11,-3.45  cm*^/GPa, 
Ilc-l.Sl  cm'i/GPa,  and  =  <»„• 


elastic  model  used  to  relate  the  axial  and  radial  stress  con^onents 
within  the  fiber  (Eqn. (6) )  is  not  applicable.  As  described  in  Section 
2.2,  the  piezo-spectroscopic  relation  given  by  Eqn. (5)  does  not  uniquely 
define  the  full  stress  tensor.  For  the  axisynsnetrical  case  of  a 
cylindrical  fiber  embedded  in  a  matrix,  an  additional  relation  between 
the  axial  and  radial  conponents  is  necessary  to  determine  the  stress 
tensor  uniquely.  A  graphical  representation  of  the  piezo-spectroscopic 
solution  for  the  Ti/ruby  conf>osite  defines  a  one-dimensional  surface  in 
stress  space  (Fig. 5).  The  elastic  model  defines  one  possible  solution, 
but  it  is  clear  that  other  solutions  may  be  indicated  by  amending  the 
constitutive  model  to  include  matrix  plasticity  and  cracking  [15] .  In 
all  cases,  however,  the  piezo-spectroscopic  data  provides  the  necessary 
physical  constraint  upon  the  solution. 


6 .  Conelttslons 

The  applicability  of  fluorescence  spectroscopy  to  stress 
measurements  in  both  the  through-focus  and  transmission  configurations 
has  been  experimentally  established  in  two  sapphire-fiber  conposites. 

For  the  case  of  the  metal-matrix  conposite,  the  load  transfer  length  was 
found  to  be  several  times  greater  than  that  for  the  intermetallic-matrix 
composite.  This  trend  is  ascribed  to  the  lower  stiffness  and  yield 
strength  of  the  titanium  con^ared  to  titanium  aluminide.  Experiments  on 
the  800-micron  thick  TiAl  con^osite  confirmed  the  validity  of  the 
transmission  fluorescence  method  through  the  strong  correlation  between 


the  axial  stress  profile  and  the  stress  measured  in  transmission.  The 
applicability  and  utility  of  the  method  were  demonstrated  on  a  long 
metal-matrix  composite  (t/R-70),  for  which  transmission  fluorescence  was 
found  to  be  the  only  means  by  which  to  estimate  the  bulk  stress  state 
and  load  transfer  length. 

Having  established  the  fluorescence  method  a*  a  viable  technique 
for  microroechanical  measurements,  current  emphasis  is  being  directed 
toward  a  strong  correlation  between  frequency  shift  measurements  and  the 
full  stress  tensor,  as  applied  to  components  of  arbitrary  geometry  and 
loading  conditions.  This  involves  the  relation  of  both  frequency  shifts 
and  line  shapes  to  the  local  stress  tensor  using  finite  element  and 
analytical  constitutive  calculations. 


Mot*:  Since  preparation  of  this  manuscript.  He  and  Clarke  [7]  have 
confirmed  that  the  piezo-spectroscopic  coefficients  of  ruby  are  distinct 
for  the  Ri  and  R2  lines.  As  a  result,  an  independent  relation  is 
obtained  from  the  frequency  shift  of  each  of  the  two  fluorescence  lines 
(Eqn.(5)).  This  should  allow  the  principal  stress  conponents  in  an 
axisymmetrical  system  to  be  uniquely  determined  without  the  need  for  a 
constitutive  law. 
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Abatfact — ^The  retidiial  stresses  in  c-axis  sapphire  fibers  in  a  y-TiAl  matrix  and  in  a  polycrystalline  ALO, 
matrix  as  a  function  of  distance  below  a  surface  are  determined.  They  are  obtained  from  the  shift  in 
frequency  of  the  characteristic  R^  fiuorescence  line  of  chromium  in  sapphire  obtained  by  focusing  an 
optical  probe  at  different  depths  in  a  sapphire  fiber  intersecting  the  surface  of  the  composite.  The  method 
is  described  together  with  its  calibration.  Both  the  axial  and  radial  components  of  the  residual  stress  in 
the  fiber  are  observed  to  rrary  over  a  depth  of  approximately  the  fiber  diameter  and  ate  then  almost 
independent  of  depth. 


1.  INTRODUCTION 

When  a  composite  material  is  cooled  down  from 
its  fabrication  temperature,  thermal  mismatch  re* 
sidual  stresses  develop,  which  can  have  a  profound 
influence  on  the  mechanical  properties  of  the  com¬ 
posite.  For  example,  in  many  ceramic  fiber  reinforced 
metal  matrix  or  glass  matrix  composites  the  thermal 
expansion  coefficient  of  the  matrix  a.  is  con- 
skferably  larger  than  that  of  the  fiber  %  and  so  large 
compressive  stresses  can  develop  in  the  fiber.  The 
axial  cmnpressive  stress  creates  a  large  stress  differ¬ 
ence  between  debonded  and  bonded  regions  of 
the  interface  when  the  material  is  under  a^ilied 
tension,  facilitating  interfadal  crack  growth  [1,2], 
a  pre-requisite  for  fiber  pull-out.  On  the  other  hand, 
the  large  radial  compressive  stress  produces  large 
friction  forces  in  the  debonded  region,  which  retards 
fiber  sliding  and  consequently  shortens  fiber  pull¬ 
out  length  [2, 3).  It  is  therefore  possible  to  optimize 
the  mechanical  properties  of  such  composites  by 
tailoring  the  residual  stresses  through  appropriate 
processing  conditions  and  the  use  of  interfadal 
coatings  [4]. 

Unfortunately,  accurate  information  about  the  re¬ 
sidual  stresses  in  embedded  fibers  is  usually  unavail- 
aUe  since  thdr  measurement  is  particularly  difficult. 
As  a  result,  estimations  based  on  the  thermal  and 
elastic  properties  of  the  individual  material  com- 
ponoits  are  usually  used  in  various  analyses.  This 
approach  requires  accurate  knowledge  of  the  effective 
temperature  change  AT,  which  is  usually  an  unmea¬ 
surable  quantity,  and  is  difficult  to  estimate  without 
detailed  knowMge  of  the  stress  relaxation  mechan¬ 
isms,  especially  near  the  interface.  Furthermore,  such 
estimates  can  not  take  into  account  the  variations 
from  fiber  to  fiber  due  to  irregularities  in  coatings 
and  interfadal  roughness:  details  which  are  import¬ 
ant  in  understanding  individual  fiber  pull-out/push- 


out  tests.  Direct  measurement  of  the  residual  stresses 
in  fibers  is,  in  prindple,  possible  by  a  variety  of 
diffraction  techniques.  For  instance.  X-ray  diffraction 
using  ultra-fine  beams  can  provide  strain  in¬ 
formation  of  individual  fibers.  However,  it  suffers 
from  having  a  relatively  small  penetration  depth 
(— 50^m)  (5,6]  and  can  therefore  only  measure 
stresses  near  the  surface,  which  could  be  significantly 
different  from  those  deeper  inside  the  composite. 
Neutron  scattering,  whilst  capable  of  probing  fiber 
stresses  throu^out  the  composite  (7,8],  has  in- 
sufiident  spatial  resolution  to  study  individual 
fibers  or  to  ascertain  stress  distributions  along  the 
fibers. 

An  alternative  method  of  measuring  stresses  in 
fibers  is  by  using  piezo-spectroscopic  effects.  The 
basis  of  such  methods  is  that  the  characteristic  lu¬ 
minescence  lines  (Raman  or  fluorescence)  of  the 
material  diiff  with  stress.  By  utilizing  an  optical 
microscope  to  define  a  probe,  it  is  relatively  easy  to 
select  regions  of  interest  with  spatial  resolution  of 
only  a  few  microns  (9]  dther  on  the  surface  of  a 
composite,  or  if  the  matrix  is  tran^Mirent,  within  the 
composte.  For  example,  in  the  case  of  polydi- 
acetylene  fibers  embedded  in  an  optically  transparent 
epoxy  (10],  Raman  spectra  were  collected  from  differ¬ 
ent  positions  along  a  fiber  to  obtain  the  stress 
distribution  along  it.  Since  optical  spectra  can  be 
collected  and  analyzed  rvith  relative  ease  compared  to 
most  diffraction  methods,  such  measurements  can  be 
performed  routinely. 

In  this  paper,  we  present  a  new  method  of 
measuring  residual  stresses  in  embedded  transparent 
fibers,  one  that  utilizes  optical  fluorescence,  and 
illustrate  its  application  to  determining  the  residual 
stresses  in  sapphire  fibers  in  y-TiAl  and  AljO, 
matrices.  The  theoretical  basis  of  the  fluorescence 
method  3vas  described  in  detail  in  an  earlier  work 
(11].  Briefly,  the  frequency  shift  of  a  characteristic 
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Fig.  1.  Schematic  diagram  of  the  optical  anangement  used  to  collect  fluoreioence  from  a  region  in  an 
embedded,  transparent  fiber  below  iu  surface.  The  stress  induced  shift  in  frequency  of  a  characteristic 
luminescence  line  is  used  to  determine  the  stress  in  the  probed  region. 


fluorescence  line  is  stress  dependent  and  is  given 
by 

Av  «  n^jtr  J  (I) 

where  11,^  ate  the  piezo>spectroscopic  coefficients 
relating  frequency  shift  to  stress,  and  aH  is  the  stress 
tensor  represented  in  the  crystallographic  frame  of 
reference. 

The  method  is  to  use  an  optical  microscope  to 
focus  a  light  beam  on  the  end  of  an  exposed  sapidiire 
fiber  (Fig.  I),  measure  the  frequency  shift  of  a 
characteristic  luminescence  liiK  and  then  repeat  the 
measurement  at  successively  deeper  points  of  focus 
into  the  fiber.  Using  the  lenses  of  the  microscope  to 
form  an  optical  probe  enables  the  luminescence  to  be 
collected,  for  any  particular  point  of  focus,  from  a 
relatively  small  depth  of  field.  This  faciliutes  hi^ 
spatial  resolution  (12, 13]  along  the  fiber  but  at  the 
same  time  means  that  the  observed  frequency  shifts 
have  to  be  deconvoluted  using  a  previously  estab¬ 
lished  calibration  of  the  microscope  convolution 
function  (the  depth  of  field  function).  This  is  de¬ 
scribed  in  the  following  section.  Also  described  is  a 
test  of  the  accuracy  of  the  stress  profiling  method  in 
which  the  stress  through  the  thkkiiess  of  a  ruby  beam 
under  four-point  loading  is  measured  and  compared 
with  the  known  elastic  solution. 

To  illustrate  the  application  of  the  technique  we 
have  measured  the  residual  stresses  in  Saphikon 
sapphire  fibers.  They  have  sufficient  levels  of  chro¬ 
mium  impurity  to  produce  sufficiently  intense  fluor¬ 
escence  lines  R,  and  R,  of  ruby  when  excited  by  argon 

tThe  general  procedures  for  optical  fluorescence  measure¬ 
ment  and  analysis  were  de^bed  in  detail  in  (I  I]  and 

therefore  not  included  in  this  work. 


ion  laser.  Because  the  fibers  are  single  crystal  with  their 
c-axis  along  the  fiber  axis,  the  residua)  stresses  inside 
the  fiben  are  expected  to  have  cylindrical  symmetry. 
As  a  result,  there  are  only  two  independent  stress 
components,  the  radial  cooqwnent  tr,  and  the  axial 
component  e,.  This  simplifim  equation  (1)  to  be 

Av  *  2n,o,  +  (2) 

where  n.  and  11,  are  the  piezo-spectroscopy 
coefficients  for  the  a  and  e  direinions.  They  have 
been  measured  to  have  values  of  2.70  and 
2.15cm~'GPa~',  respectively  [14].  A  positive  shift 
implies  tension  and  a  negative  one  implies  com- 
jxession.  In  Section  3,  measurements  of  line  shifts 
along  fibers  in  both  y-HAl  and  AI2O,  matrices  as  a 
function  of  depth  are  made,  ffinoe  the  frequency  shifts 
are  dependent  on  both  the  axial  and  radial  stress 
components,  their  sqraration  requires  the  solution  of 
the  appropriate  dastk  equations.  This  is  described, 
using  the  results  of  a  finite  element  computation  in 
Section  4.  From  this,  the  stresses  as  a  function  of 
depth  are  determined. 

r  THE  THKOUGH-roCUS  DEPTH 
PR(»1L1NC  METHODt 

2.1.  Depth  of  field  fmctUm 

The  depth  resolution  of  the  c^tica)  mkroprobe 
used  for  the  fluorescence  measurements  is  determined 
by  the  numerical  aperture  of  the  objective  lens  and 
tte  size  of  the  colfection  aperture  u^  (12, 13].  As 
illustrated  in  Fig.  1,  when  the  exciution  laser  is 
focused  at  a  distance  z  below  the  top  surface  of  a 
transparent  materia)  with  uniform  concentration  of 
fluorescing  species,  only  the  signal  generated  in  a 
small  volume  above  and  below  the  focal  plane  is 


QING  MA  ud  CLARKE:  RESIDUAL  STRESSES  IN  SAPPHIRE  COMPOSITES 


1819 


COM  Of 


Fig.  2.  niustntion  of  the  effect  of  the  finite  depth  of  field 
of  the  optical  system  used.  When  the  excitation  laser  is 
focused  at  a  distance  z  below  the  top  surface,  only  the  signal 
generated  in  a  small  volume  above  and  below  the  focal  plane 
is  effectively  collected. 


collected.  The  total  signal  collected  is  then  the  inte¬ 
gral  of  the  fluorescence  from  different  depths 

g(M)dM  (3) 

where  ;(»)  is  the  depth  of  field  function  of  the  lens 
used  to  describe  the  relative  collection  efficiency  as  a 
function  of  u.  the  distance  from  the  focal  plane. 

By  moving  the  focal  plane  from  a  position  above 
the  top  surface  to  a  position  below  the  surface  of  a 
thick  sapphire  disk  (/  -*  oo)  containing  adequate  chro¬ 
mium  concentration,  and  recording  the  fluorescence 
fine  intensity  as  a  function  of  depth  z,  the  depth 
resolution  function  g  is  obtained  from  equation  (3), 
namely 


g(z)  = 


m) 

dr 


(4) 


The  depth  resolution  function  for  our  optical  mi¬ 
croprobe  using  a  SO/O.SS  objective  lens  and  a  SO  pm 
collection  aperture  was  measured  using  the  above 
method  and  is  shown  in  Fig.  3.  As  can  be  seen,  the 
collection  efficiency  is  maximum  at  the  focal  plane, 
and  decrease  rapidly  with  distance  on  either  side  of 
the  focal  i^ne. 


2.2.  Measurement  of  the  frequency  shift  depth  profile 

We  now  suppose  that  there  exists  a  stress  field  that 
only  varies  with  depth  z.  From  equation  (2),  the  shift 
of  a  fluorescence  line  associated  with  the  increment  in 
signal  from  a  narrow  slice  of  material  at  depth  z  -I-  u  is 

Av(z  -I-  u)  a  2n,0,(z  -I-  M>-|-  n,ff.(z  ■¥  «).  (5) 

The  measured  shift  with  focal  point  positioned  at 
depth  z  is  then  a  weighted  average  of  the  signals  from 
all  depths 

_  I  Av(z -I- «)g(u) dir 

Av(z)aiLri  -  - - .  (6) 

g(M)d« 


Fig.  3.  The  depth  of  field  function  of  the  microprobe  used. 
It  describes  the  relative  collection  efficiency  as  a  function  of 
distance  from  the  focal  plane.  Details  ate  given  in  the  text. 


Provided  the  depth  of  field  function  g(z)  is  known, 
the  real  depth  profile  Av(z)  can  be  deconvoluted  from 
equation  (6). 

To  establish  that  the  foregoing  methodology  is 
correct,  the  stress  field  through  the  thickness  of  a 
dilute  ruby  beam  in  4-point  loading  was  measured 
and  compared  with  the  elastic  solution  for  the  stress. 
A  numbw  of  fluorescence  spectra  were  taken  as  the 
laser  beam  was  focused  at  different  depths  from  the 
top,  tension  surface  to  the  bottom,  compression 
surface  of  the  beam.  The  measured  shifts  are  plotted 
in  Fig.  4  as  the  solid  dots.  The  deconvoluted  data  is 
shown  by  the  solid  line.  The  deconvolution  was 
obtained  by  assuming  that  the  real  shifi  Av(z)  is  a 
linear  fuiKtion  of  depth  a  +bz  since  it  is  known  that 
the  stress  distribution  varies  linearly  with  the  depth  z 
for  a  rectangular  beam  under  4-point  hading.  This 
function  is  then  convoluted  with  the  depth  of  field 
function  g(z).  By  fitting  the  resulting  convoluted 
function  to  the  measured  shifts,  the  coefficiente  a  and 
h  are  obtained.  The  initially  assumed  function,  with 
the  coefllcients  obtained  by  fitting,  is  then  regarded  as 
the  deconvoluted  shift  Av(z)  directly  related  to  the 
stress  distributions  by  equation  (2).  The  straight  solid 
line  in  Fig.  4  therefore  represents  the  real  sUfts.  As 


Fig.  4.  The  fluorescence  line  shift  in  the  region  between  the 
tension  surface  (:  -  0)  and  the  compression  surface 
{;  •  513  fint)  of  a  ruby  beam  under  4-point  bending.  The 
solid  line  representing  the  true  shift  is  obtained  by  deconvo- 
luting  the  measured  shift  with  the  depth  of  field  function. 
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expected,  and  required,  the  zero  shift  position  co¬ 
incides  with  the  neutral  plane  of  the  beam.  The  stress 
distribution  calculated  by  using  equation  (2)  is  also  in 
good  agreement  with  that  obtained  from  beam  theory 
and  loading  conditions.  As  a  further  check,  the 
convoluted  shift,  obtained  by  convoluting  the  derived 
values  of  the  real  shift  with  the  instrument  function 
f  (;)  is  shown  by  the  dashed  curve.  It  has  excellent  fit 
with  the  measured  data  except  for  the  slight  deviation 
near  the  compression  surface  which,  we  believe,  is  due 
to  the  fact  that  the  focus  deteriorates  after  light  has 
passed  all  the  way  through  the  ruby  beam.  These 
observations  thus  confirm  that  the  through-focus 
method  can  indeed  accurately  measure  the  depth 
profile  of  the  stresses. 

3.  MEASURED  FREQUENCY  SHIFTS  FROM 
EMBEDDED  SAPPHIRE  FIBERS 

The  depth  profiling  method  has  been  applied  to 
two  different  sapphire  fiber  composites:  the  first,  a 
y-TiAl  matrix  composite  with  carbon  black/alumina 
double  fiber  coatings  on  the  sapphire  fibers  and,  the 
second,  a  polycrysulline  AI^O)  matrix  with  a  thin 
molybdenum  coating  on  the  fibers.  In  each  case,  the 
fiber  diameter  is  about  120  fim.  Thin  slices  of  com¬ 
posite  with  parallel  surfaces  perpendicular  to  the 
sapphire  fibers  were  cut  to  about  4S0  ^m  and  600  ftm 
thick  for  the  y-TiAl  matrix  and  the  AI2O)  matrix 
composites  respectively.  The  top  surface  of  the  spec¬ 
imens  were  then  polished  to  obtain  an  optical  finish 
at  the  fiber  ends. 

The  measured  and  the  deconvolutsd  frequency 
shiA  as  a  function  of  depth  into  a  san>hire  fiber  in 
V-TiAl  matrix,  as  well  as  the  convoluted  function  are 
plotted  in  Fig.  S.  Since  y -TiAl  has  considerably  larger 
thermal  expansion  coefficient  than  saf^hire  (Table  1), 
the  residual  stress  developed  by  thermal  mismatch  in 
the  fiber  should  be  compressive  in  both  radial  and 
axial  directions.  This  is  consistent  with  the  large 
negative  shifts  observed  when  the  probe  is  focused 
deep  in  the  interior  of  the  specimen.  Near  the  surface, 
the  axial  stress  approaches  zero  as  required  by  force 
balance.  The  deconvoluted  shift  beromes  positive 
near  the  surface  indicating  that  the  radial  stress 


Tabic  I 


Propenv 

Sapphire 

y-TiAl 

AUO, 

£(GPa) 

434 

173 

250 

V 

o.:7 

033 

0.27 

1  (10  *'  C) 

8.J«i).9.0<r( 

11-13 

8.3  <  a  <  9.0 

becomes  tensile  near  the  surface.  In  the  case  of  the 
sapphire  fibers  in  polycrystalline  AI.Os.  the  residual 
stress  is  mainly  caused  by  the  anisotropy  of  the 
thermal  expansion  coefficient  of  sapphire.  As  shown 
in  Fig.  6,  the  relatively  large  negative  shifts  near  the 
surface  indicate  that  the  radial  stress  near  the  surface 
is  compressive,  since  the  axial  stress  must  approach 
zero  due  to  force  balance  requirement,  and  thus  does 
not  contribute  to  the  shift.  The  magnitude  of  the 
negative  shift  becomes  progressively  smaller  with 
increasing  depth  implying  that  the  negative  shift 
associated  with  radial  compressive  stress  is  partially 
compensated  by  the  positive  shift  resulting  from  the 
build  up  of  the  axial  tensile  stress.  These  observations 
ate  consistent  with  the  fact  that  the  thermal  expan¬ 
sion  coefficient  of  polycrystalline  AKOj  is  larger  than 
that  of  the  sapphire  in  the  a -direction,  but  smaller 
than  that  in  the  r-direction.  A  notable  feature  of  the 
depth  profile  of  the  shift,  in  both  systents,  is  the  near 
exponential  dependence.  This  is  indicative  of  the 
fibers  being  fully  bonded  to  the  matrices  and  the 
fiber/matrix  systems  being  elastic. 

4.  ANALYSIS 

The  residual  stresses  in  the  sapphire  fibers  are 
related  to  the  observed  frequency  shifts  presented  in 
the  previous  section  by  equation  (2).  In  order  to 
obtain  actual  values  of  the  two  stress  components,  an 
additional  relationship  between  the  two  components 
is  evidently  required.  Two  difficulties  are  associated 
with  establishing  such  a  relationship.  First,  the  re¬ 
sidual  stresses  originate  from  thermal  expansion  mis¬ 
match  but  the  temperature  on  cooling  at  which  stress 
relaxation  no  longer  is  effective,  and  stresses  begin  to 
build  up,  is  not  known.  Secondly,  existing  equations 
relating  the  residual  stresses  in  a  fiber  embedded  in  an 
infinite  matrix  [2]  to  the  thermal  mismatch  strain 


Fig.  5.  The  line  shift  depth  profiles  (r  0  at  the  surface)  in 
a  c-axis  sapphire  fiber  embedded  in  the  y-TiAl  matrix. 


Fig.  6.  The  line  shift  depth  profiles  (r  =  0  at  the  surface)  in 
a  r-axis  sapidiire  fiber  embedded  in  a  polycrystalline  AIjO, 
matrix. 
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Aa  AT  aic  unfortunately  not  valid  for  the  thin  speci- 
mena  uted  in  this  work.  Furthermore,  rigorous  ana¬ 
lytical  solutions  for  the  residual  stress  field  are  not 
available  and  the  an>roximate  models  used  in  analyz¬ 
ing  stresses  in  fibm,  e.g.  those  in  analyzing  fiber 
pull-out/push-out  tesu  [13-17]  breakdown  near  the 
fiber  end.  We  have  therefore  had  to  resort  to  finite 
element  calculations  of  the  residual  stress  field  for  the 
sapphire/y-TiAl  and  sapphire/ AIjO]  systems  and  use 
an  iterative,  self  consistent  procedure  to  calculate 
the  residual  stresses  and  the  temperature  at  which  the 
stresses  become  fro«n  in.  This  is  described  in  the 
following  paragraphs.  First,  however,  it  should  be 
noted  that  the  stresses  can  be  estimated  directly  from 
equation  (2).  For  the  y-TiAl  matrix  system,  the 
residual  stresses  in  the  fiber  are  expected  to  be 
compressive.  If  it  is  assumed  that  the  magnitude  of 
the  radial  component  is  similar  to  the  axial  com¬ 
ponent  in  the  interior  of  the  composite,  then  using  the 
shiA  from  deep  in  the  fiber  we  have  from  equation  (2) 


ar 


Av^" 

2n,-i-n, 


300  MPa. 


For  the  AI2O]  matrix  system,  due  to  the  anisotropy 
of  the  tberml  expansion  coeffident  of  sapphire,  the 
fiber  should  be  under  axial  tension  and  radial  com¬ 
pression.  If  we  assume  that  the  radial  compression 
stress  is  about  the  same  in  the  interior  aitd  near  the 
surface,  we  again  can  estimate  that 
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200  MPa. 


Quantitative  stress  distributions  can  be  obtained  by 
comparing  in  detail  the  measured  shiAs  with  those 
calculated  from  the  elastic  field  in  the  fiber/matrix 
system.  Since  the  fiber  coatings  are  relatively  thin  in 
both  composites,  we  expect  their  effects  on  tte  elastic 
field  to  be  negligible  compared  to  those  of  the  fiber 
and  matrix.  Therefore,  in  the  calculations,  the  elastic 
system  is  simplified  as  consisting  solely  of  an  isotropic 
matrix  of  modulus  and  Poisson's  ratio  v.  and  a 
fiber  of  modulus  E,  and  Poisson's  ratio  v,  (Table  1). 
The  qjecimen  was  represented  by  a  thin  disk  of  the 
same  thickness  as  the  corresponding  specimen  with 
one  cylindrical  fiber  sitting  at  the  synunetry  axis.  The 
diameter  of  the  disk  was  chosen  to  be  about  ten  times 
that  of  the  fiber  in  order  to  be  consistent  with  the  low 
fiber  volume  fraction  of  both  composite  systems.  The 
thermal  expansion  coefficients  of  the  sapphire  fiber  in 
both  the  radial  direction  a,  and  axial  direction  a,  are 
also  listed  in  Table  1.  The  mean  thermal  expansion 
coefficient  for  y-TiAl  is  known  to  be  1 1-13  x  I0~*/°C 
depending  on  the  stress  frozen-in  temperature.  How¬ 
ever,  that  of  the  polycrystalline  AljOa  is  expected  to 
vary  from  sample  to  sample  deprading  on  the  re¬ 
sidual  matrix  porosity.  As  a  result  of  the  complicated 


stress  relaxation  processes  that  occur  during  cooling, 
particularly  near  the  interface,  the  effective  thermal 
expansion  coefficient  a.  for  both  matrices  and  the 
effective  temperature  change  AT  were  adjusted  in  the 
compuution  to  o^in  the  best  agreement  with  the 
experimenul  frequency  shift  data.  The  following 
procedure  was  adopted.  First,  reasonable  estimates 
for  the  matrix  thermal  expansion  coefficient  and 
temperature  change  a™  and  AT"  were  used  in  the 
finite  element  calculation.  The  resulting  stresses 
e"(z)  and  were  then  used  in  equation  (2)  to 
calculate  the  shiA  as  a  function  of  r,  AvS(z).  This 
function  and  the  measured  shifts  are  then  normalized 
by  their  values  in  the  interior  of  the  material 

Av 

(7) 


By  comparing  the  calculated  and  measured  normal¬ 
ized  shifts  for  a  series  of  values  of  a",  the  value  that 
gives  the  best  agreement  is  considered  as  the  true  at.. 
Then,  by  requiring  that  the  calculated  shiA  be  equal 
to  the  measured  shiA  in  the  interior,  the  effective 
temperature  difference  is  proportionally  scaled 

Av" 

AT-AT-^.  (8) 
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Finally,  the  stress  distributions  are  obtained 


'  Ar- 
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(9) 


Since  the  thermal  expansion  coefficient  of  y-TiAl  b 
considerably  larger  than  ot«  and  a,  of  sapphire,  the 
diffnence  between  a,  atxl  a,  was  unimportant.  As  a 
result,  it  was  observed  that  the  fitting  of  the  normal¬ 
ized  shiA  was  not  sensitive  to  the  actual  value  of  at.. 
Figure  7  illustrates  the  fitting  between  the  measured 
and  the  calculated  normal  shiAs  for  at.  - 12.  The 
effective  temperature  drop  AT  calculated  from 
equation  (8)  is  '~730°C,  which  agrees  with  at.-T 
curve  for  y-TiAl  [18].  This  temperature  change  is 
about  300°C  smalls  than  the  difference  between  the 


Fig.  7.  A  comparison  of  the  calculated  and  the  measured 
nonnaUzed  shifts  as  a  function  of  depth  for  the  sapphire/y- 
TiAl  system. 
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Fig.  8.  The  radial  and  axial  stresses  as  a  function  of  depth 
in  the  sapphire/y-TiAI  system  determined  from  the  Suor- 
esoenoe  data. 


Fig.  9.  A  comparison  of  the  calculated  and  the  measured 
normalized  shifts  at  a  function  of  depth  for  the  tap- 
phite/AljO,  system. 


processing  temperature  and  the  room  temperature 
indicating  significant  stress  relaxation  occurring 
during  cooling.  The  stress  distributions  along  the 
fiber  calculated  from  equation  (9)  are  plotted  in 
Fig.  8.  The  axial  stress  approaches  «ro  as  expected 
and  the  radial  component  is  positive  at  the  surface  as 
was  indicated  by  the  frequency  shift  measurement. 
The  diflerences  between  the  calculated  and  the 
measured  normalized  shifts  as  shown  in  Fig.  7  are  due 
to  the  simplifications  used  in  the  model,  e.g.  neglect¬ 
ing  the  eflects  of  coatings  and  assuming  perfect 
interfadal  bonding.  However,  the  finite  element  cal¬ 
culations  shows  that  the  stress  distributions  in  the 
inlerkM'  of  the  specimen  are  not  sensitive  to  small 
changes  in  the  value  of  the  elastic  modulus  and  the 
thermal  expansion  coefikaent  of  the  matrix.  This  is 
because  the  values  of  the  stresses  are  bounded  by  the 
measured  shifts  through  equation  (2)  and  scaled  by 
equatkms  (8)  and  (9).  These  observations  suggest  that 
the  stress  distributions  deduced  from  the  fluorescence 
datt  and  ittostiated  in  Fig.  8  are  relatively  accurate 
especially  in  the  interior  of  the  specimen.  The  possi¬ 
bility  of  some  interfadal  debonding  will  be  discussed 
in  more  detail  in  the  next  section. 

As  expected,  the  normalized  shift  calculated  for  the 
AljO,  matrix  system  depends  much  more  sensitively 
on  the  value  of  a.  used  in  the  computation.  As  shown 
in  Fig.  9,  a. «  8.73  gives  the  bmt  fit  between  the 
calculated  and  the  measured  normalized  shifts.  The 
effective  temperature  difletence  obtained  using 
equation  (8)  is  about  1300°C,  which  is  about  200*’C 
lower  than  the  processing  temperature  again  as  indi¬ 
cation  of  stress  relaxation  occurring  during  cooling. 
The  stress  distributions  in  the  sapphire  fiber  are 
plotted  in  Fig.  10.  The  observation  that  the  calculated 
stress  distributions  in  the  interior  are  not  sensitive  to 
changes  of  also  applies  in  this  system.  The  values 
of  strem  in  the  interior  differ  less  than  5%  for  the 
three  ajt  shown  in  Fig.  9.  Since  the  shape  of  the  riiift 
profile  depends  sensitively  to  the  change  of  a.,  the 
fitting  in  Fig.  9  also  provides  means  of  estimating  the 
local  thermal  expansion  coefikmt  of  the  Al]03 
matrix,  which  could  be  different  from  the  average 


value  of  the  whole  sample  because  the  density  of  the 
matrix  is  slightly  lower  near  the  fibers.  The  resulting 
value  of  0.  can  then  be  used  in  other  analyses,  such 
as  for  push-out  tesu. 

S.  MSCUSSKW 

The  analysis  in  the  last  section  is  based  on  the 
assumption  that  the  fiber/matrix  systems  are  dastic 
without  any  debonding  occurring  at  the  interface. 
This  requires  that  the  interfadal  shear  strength  be 
larger  than  the  maximum  interface  shear  stress.  To 
examine  this  assumption,  the  interfadal  shear  stress 
is  also  calculated  in  the  finite  element  computation. 
The  resuhs  correqxmding  to  the  best  fit  conditions 
determined  in  the  last  section  are  plotted  for  the 
sapphire/y-TiAl  and  the  sapphire/AljOj  systems  in 
Fig.  1 1.  In  both  systems,  the  magnitude  of  the  shear 
stress  has  a  maxiinimri  at  the  surface  and  decreases  to 
zero  with  increasing  depth.  This  is  expected  from  the 
stress  distributions  in  Figs  8  and  10  and  the  force 
balance  requirement 


ftfdg. 

Tdr 


(10) 


Whilst  the  calculated  interfadal  shear  stress  for  the 
sapphite/Al20}  system  is  smaller  than  the  interfadal 


Fig.  10.  The  radial  and  axial  stresses  as  a  function  of  depth 
in  the  sappliire/Al20,  system  determined  from  the'  fluor¬ 
escence  data. 
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Oiplh  Mow  Pfeor  Top  SmImo  x  Itun) 

Fig.  11.  The  interfadal  thear  ttrenes  ia  the  sapphire/r-TiAl 
ud  the  sapplure/AliOi  system*  colculoted  from  the  finite 
elemrat  oompt^tion  ud  the  fiuoresoenoe  data. 


not  significantly  affected  by  possible  debonding  near 
the  surface. 

fi.  CONCLUSIONS 

By  utilizing  the  optical  fluorescence  from  chro¬ 
mium  impurities  in  sapphire  fibers,  we  have  deter¬ 
mined  the  residual  stress  distributions  in  embedded 
fibers  in  both  a  y-TiAl  matrix  and  a  polycrystalline 
AI2O]  matrix  composite.  The  axial  and  radial  com¬ 
ponents  of  the  residual  stress,  produced  by  differen¬ 
tial  thermal  contraction  on  cooling  after  fabrication, 
in  the  fibers  are  determined  as  a  function  of  depth 
below  the  composite  surface.  For  both  composites, 
the  two  stress  components  vary  over  a  distance  of 
about  the  fiber  diameter  from  the  surface  and  then 
are  almost  independent  of  distance  further  into  the 
composite. 
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NO0OI4-92-J-18O8  (DRC).  The  authors  are  grateful  to  Dr  T. 
J.  Mackin  and  J.  I^vis,  UCSB,  for  providing  the  composite 
samples  used  in  this  work. 
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Abstract 


The  high  temperature  processing  of  fibrous  metal  matrix  composites  has 
been  jnedicted  to  result  in  the  development  of  large  residual  stresses  because  of  a 
difference  between  the  coefficients  of  thermal  expansion  (CTE)  of  the  fibers  and 
matrix.  Efforts  to  measure  the  stresses  and  verify  this  explanation  for  their 
development  have  been  hampered  by  the  lack  of  a  convenient  measurement 
methodology.  We  have  used  a  fibo*  optically  coupled  method  to  measure  both 
the  R1  and  R2  optical  luminescence  line  shifts  of  single  crystal  a-'Al203  fibers 

embedded  in  a  CP  titanium  matrix  composite.  The  line  shifts  depend  linearly  upon 
weighted  sums  of  the  fiber's  {vinciple  stresses,  and  we  find  the  average  measured 
shifts  for  all  the  Hbers  in  a  composite  sample  to  be  roughly  in  agreement  with 
concentric  cylinder  residual  stress  model  predictions.  However,  significant  fiber  to 
fiber  variability  is  also  observed.  Using  both  the  R-line  measurements  and  the 
constraining  relation  that  within  the  fib^  the  radial  and  hoop  stresses  are  equal, 
we  have  been  able  to  independently  estimate  the  fiber’s  axial  and  radial  stress 
components.  The  axial  stress  is  found  to  be  mote  than  twice  that  predicted  by  a 
CTE  difference  model,  whilst  the  radial  (and  hoop)  direction  stresses  are  smaller 
than  expected.  We  propose  that  the  stress  in  these  composites  is  affected  by  two 
inesently  unmodelled  effects:  an  elevation  of  composites  residual  stress  due  to  the 
CTE  difference  between  the  tooling  used  for  (hot  isostatic)  processing  of  the 
composite,  and  a  relaxation  of  the  transverse  stress  components  for  some  fibers  by 
radial  matrix  cracking.  The  result  has  important  consequences  for  the  manufacture 
of  metal  matrix  composites  and  for  efforts  to  predict  their  mechanical  properties. 
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1.0  Introduction 

Metal  matrix  composites  (MMCs)  reinforced  with  continuous  ceramic 
fibers  have  been  the  subject  of  much  recent  research  and  development  due  to 
their  favorable  combination  of  mechanical  properties  (e.g.,  high  specific  stiffness, 
ambient  and  high  temperature  strength,  etc.)  [1].  These  traits,  together  with 
potentially  superior  creep  resistance,  make  this  class  of  materials  promising 
candidates  for  high  temperature  structural  applications  such  as  aircraft  engines  or 
the  leading  edges  of  proposed  hypersonic  vehicles  [2,3].  Several  issues  must, 
however,  be  resolved  before  MMCs  can  be  put  into  high  temperature  service. 
One,  of  interest  here,  relates  to  the  usually  large  residual  stress  that  forms  on 
cooling  from  the  high  temperatures  used  for  composite  processing  and/or  during 
the  thermal  cycling  that  accompanies  their  subsequent  use  [3]. 

It  is  normally  assumed  that  thermal  stresses  arise  in  MMCs  because  of  a 
coefficient  of  thermal  expansion  (CTE)  mismatch  between  a  composite's  matrix 
and  fiber.  Because  the  fibers  exhibit  a  lower  CTE  than  the  matrix,  they  impede 
matrix  contraction  near  the  fiber  surface  on  composite  cooling,  resulting  in  the 
development  of  stress  both  in  the  fiber  and  the  surrounding  matrix.  For  some 
fiber-matrix  systems,  the  CTE  difference  can  be  a  factor  of  two  or  more,  and  the 
matrix  stresses  are  sufficiently  large  that  they  exceed  the  matrix  yield  strength,  or, 
in  less  ductile  materials,  its  fracture  strength  resulting  in  matrix  plasticity  and/or 
radial  and  axial  matrix  cracking  on  cooling  [1,2,3]. 

These  internal  stresses  and  their  associated  damage  can  be  responsible  for 
{xremature  failure  in  service,  and  are  a  significant  source  of  quality  loss  during 
manufacturing.  Even  in  the  event  that  the  thermal  stresses  alone  are  not  large 
enough  to  cause  damage  on  first  cooling,  the  thermal  component  adds  to  that 
applied  mechanically  during  thermal-n^hanical  cycling  and  results  in  a  short 
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thermomechanical  fatigue  life  [3].  Understanding,  measuring,  and  learning  how  to 
minimize  these  thomal  stresses  is  essential  as  we  seek  reliable  ways  of 
manufacturing  MMCs  and  extending  their  high  temperature  capability. 

Many  models  have  been  proposed  for  calculating  the  residual  stresses  of 
MMCs  due  to  the  CTE  mismatch  [2,7-9].  They  range  from  the  simplest:  an 
elastic  analysis  of  two  concentric  cylincters  [7],  up  to  the  finite-element  method 
using  a  matrix  constitutive  behavior  that  includes  both  elastic  and  the  inelastic 
deformation  due  to  matrix  plasticity  and  creep  [8].  Recendy,  Pindera  and  Freed 
have  developed  an  elastic-plastic  multiple  concentric  cylinders  model  capable  of 
calculating  both  fiber  and  matrix  stresses  [9].  It  has  been  implen^nted  on  a 
personal  computer  allowing  the  rapid  convenient  calculation  of  stresses. 

Figure  1  schematically  shows  a  Pindera  and  Freed  model  predicted  stress 
field  in  and  around  a  single  crystal  aluminum  oxide  fiber  embedded  in  a  titanium 
matrix  after  cooling  from  a  stress-free  (processing)  temperature  of  90(fC.  (The 
assumptions  and  physical  properties  used  to  make  this  calculation  are  presented 
later.)  We  note  that  the  model  predicts  (1)  the  stress  components  in  the  fiber  are 
all  compressive,  (2)  the  fiber's  hoop  (0^99)  and  radial  (a^rr)  stresses  are  equal,  (3) 

the  fiber's  axial  stress  is  significantly  larger  than  its  other  components,  (4)  the 
matrix  hoop  (0F*^9)  stress  is  equal  in  magnitude  but  of  opposite  sign  to  that  of 

the  fiber,  and  (5)  none  of  the  fiber's  stress  components  vary  with  position  in  the 
fiber,  whereas  both  the  matrix  radial  and  hoop  stresses  decay  with  distance  from 
the  fib^-matiix  interface. 

It  has  been  difficult  to  experimentally  verify  these  models  and  to  confirm 
that  the  fiber/matrix  CTE  difference  is  the  explanation  for  residual  stress 
development  X-ray  diffraction  can  be  used  to  measure  residual  strains  in  the 
matrix  (by  detecting  the  change  in  lattice  parameter),  but  x-rays  normally 
penetrate  only  a  few  microns  in  metals,  so  it  is  necessary  to  remove  successive 
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layers  of  material  from  the  surface  (measuring  the  strains  at  each  newly  created 
surface)  in  order  to  measure  the  strain  as  a  function  of  position  in  the  composite 
[1],  However,  removing  layers  is  both  destructive  and  changes  the  strains  at  the 
newly  created  free  surface.  Substituting  neutrons  for  x-rays  solves  the 
penetration  and  sample  destruction  problems,  but  requires  a  volume  large  enough 
(10-30  mm^)  to  provide  sufficient  diffracted  neutron  intensity  to  reliably  detect 
the  small  diffraction  peak  shifts  due  to  strain  [11].  Since  the  strain  measured  by 
this  method  is  an  average  over  the  probe  volume,  it  does  not  provide  information 
about  the  matrix  stresses  as  a  function  of  position  near  fibers.  Energy  dispersive 
x-ray  diffractometry  [12],  which  uses  high  energy  synchrotron  x-radiation  instead 
of  neutrons,  has  the  needed  penetration  capability,  reasonable  spatial  resolution 
(~50  pm)  and  precision  (±50  MPa)  provided  the  matrix  grain  size  is  small  (<10 
pm),  and  has  been  used  to  measure  matrix  stress  in  SCS-6  reinforced  Ti3Al+Nb 
alloys.  However,  neither  Neutron  or  Synchrotron  diffraction  techniques  are 
convenient  solutions  to  the  measurement  problem  because  most  researchers  do 
not  have  convenient  access  to  these  facilities,  nor  the  resources  to  mount,  what 
quickly  becomes,  a  very  significant  measurement  campaign.  Thus,  there  still 
exists  a  need  for  a  convenient,  accurate  method  for  measuring  residual  stresses 
deep  within  metal  matrix  composites. 

Ma  and  Qaike  [13]  and  Hough,  et  al  [14]  have  sought  to  exploit  the 
natural  luminescence  of  a-Al2C)3  fibers  in  order  to  do  this.  a-Al203  of  the 

purity  used  for  composite  fibers  (both  single  and  polycrystalline  forms)  contain 
sufficient  chromium  that  a  detectable  optical  luminescence  can  be  observed 
following  illumination  with  green  wavelength  light.  In  composite  grade  single 
crystal  fibers.  Hough,  et  al  [14]  have  shown  that  reasonably  efficient  light  piping 
occurs,  and  thus  luminescence  can  potentially  be  excited  and  detected  from  deep 
within  the  composite  using  fiber  optic  methods  of  coupling.  The  luminescence 
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peak  wavelength  (which  can  be  precisely  measured  with  a  modem  spectrometer) 
depends  linearly  upon  a  weighted  combination  of  the  stress  components  in  the 
fiber,  and  since  the  fiber  stress  is  in  mechanical  equilibrium  with  that  of  the  matrix, 
a  measurement  of  the  luminescence  line  shift  affords  the  opportunity  to  infer  the 
matrix  residual  stress.  However,  care  must  be  taken  with  luminescence 
measurements  and  their  interpretation  because  factors  other  than  stress  (e.g. 
temperature)  can  affect  the  measured  line  shift. 

Here  we  review  the  underlying  mechanisms  responsible  for  the 
luminescence  of  Cr(in)  doped  (to  intensify  the  luminescent  response)  a-Al203, 

and  identify  the  best  values  for  the  (piezospectroscopic)  constants  that  relate  the 
luminescence  line  shift  to  stress,  and  other  factors  that  are  responsible  for 
(anomalous)  line  shifts  due  to  temperature  and  chromium  ion  concentration 
variations.  The  design  of  a  fiber  optically  coupled  spectroscopic  system  with 
sufficient  peak  position  resolution  for  residual  stress  measurement  is  next 
described.  We  then  explain  our  composite  fabrication  method  for  producing 
samples  with  embedded  CrOU)  doped  single  crystal  AI2O3  fibers.  The  measured 

peak  shifts  for  the  fibers  in  the  "as  processed"  condition  are  compared  with  those 
predicted  from  fiber  stresses  calculated  using  the  Pindera-Freed  thermal  stress 
model,  and  the  validity  of  the  model's  predictions  for  realistic  MMCs  discussed. 

It  will  be  shown  that  the  containment  used  to  consolidate  composites  introduces 
a  large  additional  residual  stress  in  conqx)site  components  that  is  not  included  in 
the  current  models,  and  this  could  adversely  affect  composite  quality. 
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2.0  Principles  of  Luminescence 

2.1  Luminescence  in  the  Doped  AI2O3  System 

Luminescence  arises  from  the  radiative  (emission  of  a  photon)  transition 
between  two  energy  states  of  an  atomic  or  molecular  system.  The  conditions  for 
luminescent  emission  are  rather  stringent  and  only  a  relatively  small  percentage  of 
molecules  luminesce.  The  two  types  of  luminescence  (fluorescence  and 
phosphorescence)  are  characterized  by  the  multiplicities  of  the  ground  and  the 
emitting  states.  Fluorescence  is  charactraized  by  initial  and  final  states  of  the 
same  multiplicity,  and  the  transition  is  thus  not  forbidden  by  spin  selection  rules. 

It  is  a  "highly  allowed"  and  therefore  a  short  lived  process  (typically  requiring 
only  nanoseconds  to  depopulate  the  excited  state).  Phosphorescence,  on  the 
other  hand,  involves  a  spin  flip  or  change  of  multiplicity  between  the  excited  and 
ground  states.  Formally,  such  transitions  are  forbidden,  but  in  practice,  various 
mechanisms  relieve  the  forbiddenness  of  the  transition  and  make  it  somewhat 
allowed.  This  forbiddenness  causes  phosphorescence  to  have  a  long  lifetime  (i.e., 
a  long  time  is  required  to  depopulate  the  excited  state)  ranging  from  hundreds  of 
microseconds  to  even  seconds.  The  free  CrOH)  |qji  thiiee  outer  d  electrons 
designated  d^)  has  a  quartet  ground  state  and  quartet  and  doublet  excited  states. 
For  this  case  fluorescence  accompanies  the  quartet-quartet  process  and 
phosphorescence  is  seen  from  the  doublet-quartet  transition. 

The  incorporation  of  CrC®)  in  the  near  octahedral  sites  of  the  a-Al203 

(trigonal  symmetry)  lattice  is  responsible  for  the  luminescence  of  aluminum  oxide 
fibers.  The  octahedral  crystal  field  of  the  unstrained  lattice  splits  the  5-fold 
degenerate  d-orbitals  of  the  free  ion  ,  into  a  (lower  energy)  3-fold  degenerate  t2g 
orbital  and  an  upper  eg  orbital.  The  notation,  from  group  theory,  is  for  an 
octahedral  point  group,  where  the  "g"  subscript  denotes  an  orbital  that  is  ' 
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symmetric  with  respect  to  inversion  (gerade);  all  d-electron  derived  excited  states 
will  be  g  in  a  centrosymmetric  point  group.  The  d^-clcctrons  of  Cr^I^)  occupy 
the  three  t2g  orbitals  in  the  unstrained  lattice. 

The  description  of  the  optical  spectroscopy  of  ruby  builds  on  the  orbital 
and  state  diagram  of  an  ion  in  such  an  octahedral  site.  Figure  2  [IS].  This  d^ 
configuration  gives  rise  to  a  ^A2g  ground  state.  Excited  states  arise  from  spin 
flips  within  the  t2g  orbital  and  promotions  with  or  without  spin  flips  of  a  t2g  to 
an  Cg  orbital.  The  excited  t2g  configuration  yields  doublet  exci^  states 
including  a  ^Eg  state,  which  is  the  lowest  excited  state  in  the  ion;  ruby 
phosphorescence  arises  from  the  ^Eg— ^^A2g  transition.  Since  this  transition  is 
formally  spin  forbidden,  the  ruby  phosphorescence  is  also  forbidden  and  is 
therefore  long-lived  (several  milliseconds).  Similarly,  direct  absorption  to  this 
state  (^Eg«-^A2g)  is  forbidden,  and  the  transition  is  too  weak  to  be  optically 
pumped  effectively.  Two  quartet  excited  states  (^ig  and  ^2g)  arising  from  the 
t^gCg  configuration  give  rise  to  allowed  absorptions  (^lg<-^A2g  and 

^2g<”^A2g)  (shown  as  shaded  regions  in  Figure  2).  These  allowed  transitions 
occur  in  the  blue  and  green  parts  of  the  spectrum  and  give  ruby  its  red  color. 
They  absorb  intensely  enough  to  be  pumped  efficiently  by  a  variety  of  lasers 
including  ionized  Ar,  HeCd,  and  green  HeNe.  Once  excited,  the  ^ig  and  ^2g 
states  decay  nonradiatively  with  unit  efficiency  to  the  lower  ^Eg  state,  which 
then  emits  the  luminescence  of  interest  with  near  100%  quantum  efficiency. 

In  reality,  the  picture  is  a  little  more  complicated  than  presented  above. 

The  Cr(^  ion  is  located  in  an  octahedral  lattice  site  with  a  slight  trigonal 
distortion.  The  notation  for  its  energy  levels  must  then  switch  to  that  used  for 
spin  orbit  states  in  a  noncentrosymmetric  trigonal  crystal  field.  This  distortion, 
together  with  spin-orbit  coupling,  splits  the  ^Eg  state  into  an  E  and  a  2  A  state 
separatedin  wavenumber  (energy)  by  about  30  cm’ The  E  and  2  A  states 
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give  rise  to  a  doublet  in  the  phosphorescence  spectrum,  designated  R1  and  R2. 

In  the  absence  of  an  externally  imposed  lattice  strain,  these  states  are  14400  cm*^ 
and  14430  cm~^  above  the  ground  state  [15].  Decay,  via  radiative  transitions, 
thus  results  in  luminescence  with  two  spectral  lines  (the  R1  and  R2  lines, 
respectively)  with  wavelengths  of  approximately  693  nm  and  694  nm  (Figure  2). 
2.2  Luminescence  and  Stress 

Schawlow  [16],  during  development  of  the  first  ruby  lasers,  was  the  Hrst  to 
observe  that  a  compressive  stress  altered  both  the  R  peak  positions,  and 
subsequently  was  able  to  predict  this  phenomenon  using  a  quantum  mechanics 
approach  [17, 18].  At  about  the  same  tiire,  Forman  et  al  at  the  Institute  of 
Standards  and  Technology  (formerly  the  National  Bureau  of  Standards)  also 
observed  a  systematic  change  in  the  emission  wavelength  with  pressure  from  a 
small  ruby  chip  placed  in  a  diamond  anvil  cell  (DAQ  used  for  studies  of  high 
pressure  phase  equilibria  [19. 20].  The  NBS  group  examined  the  responses  of 
many  other  luminescing  crystals  and  detmnined  that  the  very  intense  peak 
emission,  low  background  intensity,  and  moderate  line  width  of  the  ruby  signal 
made  it  the  best  material  for  measuring  stress  [20].  They  calibrated  the 
hydrostatic  stress- wavelength  relationship  of  the  R1  line  up  to  20  GPa  and 
reported  it  to  be  linear  [21].  Bell  et  al  have  extended  the  calibration  to  SSO  GPa 
in  order  to  measure  the  pressure  in  DACs  of  improved  design  [22-24].  They 
reported  that  the  calibration  only  exhibited  measurable  nonlinearities  for 
pressures  above  29  GPa  [24].  (Recent  unpublished  studies  by  Clarice,  et  al 
indicate  a  small  nonlinearity  may  be  present  in  the  R1  peak  shift  for  stresses  less 
than  3GPa.) 
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The  fibers  in  a  composite  are  not  under  purely  hydrostatic  loading.  Grabner  [25] 
has  developed  a  relationship  between  the  tensor  components  of  stress  in  a 
multiaxially  loaded  crystal  to  the  scalar  line  shift: 

Av  =  HjO^  (1) 

where  Av  is  the  line  shift  due  to  the  application  of  stress  (in  cm'^),  <r^  is  the 
ai^lied  stress  defined  in  the  crystal  coordinates,  and  ny  are  the  components  of 
the  crystal's  piezospectroscopic  coefficient  tensor.  This  can  be  rewritten  in  terms 

of  principle  applied  stresses  using  the  standard  transformation  convention: 

c 

-  aimOinn^jn  (2) 

where  aifn  ajn  are  coordinate  transformation  matrices  and  amn  is  the  stress 
tensor  in  the  coordinates  of  the  principle  applied  stresses.  Applying  (1)  and  (2)  to 
a  residually  stressed  ruby  fiber  like  that  in  Hg.  1  with  its  c-axis  aligned  with  the 
O2Z  (axial)  principle  stress  yields: 

Av  *  rin  (<Jrr+<y0e)  +  (3) 

Schawlow,  et  al  [16, 17]  reported  that  FIi  1  and  1133  for  the  R1  line  have 

values  of  3.0  cm'^GFa'^  and  1.8  cm’^GPa*^,  req)ectively.  Different  values  oi  the 
riy  for  the  R1  peak  shift  have  since  been  reported  and  these,  together  with  the 

piezospectroscopic  constants  for  the  R2  line,  are  sununarized  in  Table  1.  The 
values  presented  in  Table  1  were  first  nwai  jred  by  Schawlow  in  1961  [17],  and 
shortly  thereafter  (1962)  by  Kaplyanskii  and  Przhevuskii  [26].  Feher  and  Sturge 
[27]  lat^  reported  average  (i.e.  [R1  +  R2]  2)  values  that  were  lower  than  the 
separate  Rl,  R2  values  of  other  studies.  Recently  Lipkin  and  Clarke  have 
reported  new  measurements  of  piezospectroscopic  constants  [28,  29]  that  are 
similar  to  those  of  [26] . 


Table  1:  Values  of  Fin,  1133  and  hydrostatic  loading  constants  reported  by 

various  studies. 
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Piezospectroscopic  Coefficient 

nil 

(cm^/GPa) 

033 

(cm*/GPa) 

2nii+n33 

(cm^/GPa) 

1  Source 

R1 

R2 

R1 

R2 

R1 

R2 

3.0 

2.8 

1.8 

2.3 

7.8 

7.9 

Kaplyanskii/n^vuskii  26j 

3.2 

2.8 

1.4 

1.9 

7.8 

7.5 

Upkm  &  Clarke,  et  ai  [2^  291 

3^5 

18 

1.5 

2.1 

8.4 

7.d 

Feher  and  Sturge  [27] 

2. 

IS 

7.. 

55 

It  is  interesting  to  note  that  if  line  shifts  for  the  R1  and  R2  peaks  are  accurately 
measured,  and  the  Ilij  values  for  the  two  R  lines  are  sufficiently  different  (and 

accurately  known)  one  could  construct  two  equation  (3)’s  and  use  the 
constraining  relation  Orr  =  090  in  the  fibo-  to  independently  calculate  (Txz  ^ 

(=  <100). 

Using  a  Raman  microprobe  measureimnt  system  [30],  Ma  and  Qarke  have 
recently  reported  a  stress  variation  close  to  the  polished  surface  of  a  sapphire 
fiber  embedded  both  in  poly-crystalline  alumina  and  y-HAl  [13].  They  were  able 
to  confirm  Liang  and  Hutchinson's  recent  prediction  [31]  that  the  residual  stress 
in  the  fiber  (for  fibers  oriented  perpendicular  to  the  free  surface)  is  significantly 
perturbed  near  the  sample's  free  surface  because  of  the  traction-fiee  conditions  at 
the  newly  created  free  surface.  They  ob^rved  that  the  axial  component  of  the 
fiber  stress  decreases  from  the  bulk  value  at  about  a  fiber  diameter  below  the 
surface,  to  ^  ;ro  at  the  free  surface.  They  also  found  that  the  radial  stress  varied 
strongly  along  the  fiber  length,  and  even  became  tensile  within  about  25  pm  (a 
nfth  of  the  fiber  diameter)  of  the  fiber  surface  end.  Thus,  if  measurements  of 
interior  stresses  are  to  be  obtained  with  a  flber  optic  coupling  approach,  it  is 
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impcnrtant  that  the  inteirogated  fiber  have  reasonably  good  light  piping 
c^bilities  and  that  care  is  taken  to  efficiently  couple  (align)  with  the 
inteirogated  fiber  so  that  guided  optical  wave  modes  are  excited,  and  the 
resulting  luminescence  excited  originiues  from  deep  within  the  composite. 

23  Temperature  EfTects 

The  R  lines  of  ruby  are  also  shifted  by  temperature.  Recent  reports  of  the 
temperature  effect  have  been  made  in  the  context  of  the  DAC.  Wunder  and 
Schoen  [32]  report  that  the  wavenumber  shift  of  the  R1  peak  due  to  ^lied 
stress  has  a  linear  temperature  dependence  of  -0.149  cm' Vc  near  room 
temperature  It  is  usually  assumed  that  stress  and  temperature  shifts  are  additive, 
thus  equation  (3)  becomes: 

Av  =  nii  (oii-KJ22)  +  n33O33-01^9(Tstr-Tref)  (4) 

where  T^f  is  a  reference  temperature  (e.g.,  25‘’Q  and  Tsu  is  the  temperature  of 
the  stressed  sample  when  its  line  shift  was  measured- 

More  recent  work  by  Munro  et  al  [33],  has  thoroughly  addressed  the  issue 
of  accurate  R1  and  R2  peak  position  determination  for  hydrostatic  loading  when 
the  measurement  temperature  also  varies.  They  propose  that  the  peak  positions 
are  given  by: 

Vk  =  otici  +0^2?+  a|c3(T-Tref)  (5) 

where  V]^  is  the  peak  position  (in  cm'^)  k=lX  represents  the  R1  and  R2  peaks, 
ojq,  with  i-1-3,  are  constants  for  the  two  R  peaks,  and  P  and  T  are  the 
hydrostatic  pressure  and  sample  temperature,  respectively.  Table  2  lists  the  o^q 
constants  proposed  by  Munro  et  al. 


Table  2:  Constants  ajd  proposed  by  Munro’s  et  al  for  hydrostatic 

loading  [33]. 
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- VSS. - 

otjci  (cm*^) 

ait2  (cm'^/GPa) 

aic3  (cm’^/K) 

R1 

14^2.5±0.4 

-7.59ibW 

-0.130±0.(X)6 

- R2 - 

14432.  liK).4 

-7.6i5±o.dM 

-0.122±6.006 

Hence*  the  change  in  peak  position  between  a  strained  (n^)  and  unstrained 
sanqile  when  the  measurements  are  made  at  different  temperatures  (provided 
Tief  is  the  same  for  the  stressed  and  unstressed  measurements)  is: 

Av=Vjj-  =  ak2(Pstr ■  Po) ®k3fCrstr ‘ ^ref)^ unstr " Tref)] 

=  <*k2(Pstr  ”  Po)  +  ®k3CTstr '  Timstr)  (6) 

where  pstr  is  the  pressure  in  the  stressed  sample,  Po  is  th©  pressure  in  the 
"unstressed"  sample,  Tstr  is  the  temperature  of  the  stressed  sample,  and  Tnnttr  is 
the  temperature  of  the  unstressed  piece,  and  the  0(^(pstr '  Po)  would  be 
equivaloit  to  equation  (3).  Note  equation  (6)  indicates  that  for  Tstr  -  'i'unstr 
•0.5**C  thoe  would  be  a  tenqioatute  adjustment  to  the  line  shift  of  only  -0.06  cm* 
which  is  on  the  order  of  other  experimental  errors. 

2.4  Chnnnium  Concentration  Effects 

Thm  is  also  a  shift  in  R  peak  position  with  chromium  ion  concentration 
when  the  temperature  and  ai^lied  stress  are  both  held  constant  [34],  It  has  been 
found  diat  the  line  shift  data  fit  an  equation  of  the  form: 

Av  =  99[a]  (7) 

where  AvT  is  the  line  shift  (in  cm'^)  and  [Cr]  is  die  chromium  concentration  (in 
weight  percent).  Ma  and  Clarke  [13]  have  examined  the  fundamental  origin  of 
this  relation  and  conclude  that  it  arises  from  the  strain  field  created  by 
substituting  a  Cr(iii)  ion  for  an  Al  atom  in  the  a-Al203  lattice.  B^ause  the 
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dilational  strain  sensed  by  luminescence  is  only  a  third  of  that  measured  by  x-ray 
diffnKtion,  they  argue  that  it  is  the  "image”  component  of  the  strain  field 
that  is  responsible  for  the  concentration  dependent  frequency  shift. 

The  practical  consequence  of  this  is  that  for  a  system  in  which  there  is  a 
difference  in  stress,  temperature,  and  concentration  between  measured  and 
reference  states,  the  observed  shift  in  R  line  frequency  is: 

Av  =  oq(2(Pstr "  Po)  ■*"  ^kSlH'str  “  "^unstr)  ^^[Crlstr "  [Crlimstr)  (8) 
where  [Cr]sQ>  is  the  chromium  concentration  of  tlw  stressed  sample  and  [CrlmisQ. 
is  the  chromium  concentration  of  the  unstressed  sample. 

To  summarize,  the  predicted  fiber  stresses  in  MMCs  due  to  cooling  after 
consolidation  processing  are  in  the  500  MPa  range.  From  equation  (3),  a  stress  of 
this  magnitude  results  in  an  apfHXiximate  peak  shift  of  ~3.8  cm*^  (i.e.  ~1.8  A 
change  in  wavelength).  For  the  shift  to  be  detected  spectroscopically  a  high 
resolution  (<0.2  A)  spectrometer  must  be  used,  and  temperature  corrections 
should  be  made  if  there  is  more  than  a  0  difference  in  measurement 
temperature.  It  is  also  important  that  the  composition  of  the  stressed  and 
(unstressed)  rrference  fibers  be  the  same  (or  independently  measured)  if  potential 
errors  from  equation  (7)  are  to  be  avoided.  Hnally,  an  experimental  approach  is 
required  that  interrogates  fibers  deeply  within  the  composite  so  that  the  end 
effects  reported  by  Ma  and  Qarke  [13]  and  Liang  and  Hutchinson  [31]do  not 
dominate  the  measurement. 


3.0  Experimental 

3.1  Sample  Preparation 

Two  composite  samples  were  fabricated  by  the  foil-fiber-foil  method  from 
five  sheets  of  0.25  mm  thick  CP  titanium  foil  (supplied  by  Alfa,  Johnson  Matthey 
Products)  and  three  Ci^^-doped  single-crystal  AI2O3  fibers.  The  Hbers  for  one 
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composite  sample  (COMPS)  were  nominally  doped  with  0.03wt%  and 
wm  sui^lied  by  Saphikon,  Inc.  They  were  grown  by  the  edge  defined  film  fed 
growth  (EFG)  technique  [35].  The  attenu^on  coefficient  of  this  fiber  was 
measured  by  the  cutback  technique.  Using  the  632.8  nm  line  of  a  HeNe  laser,  the 
attenuation  coefficient  of  a  bare  Hber  was  measured  to  be  about  0.52  cm*^. 

Hbers  for  a  second  composite  sample  (COMPS)  were  nominally  doped  with  0.05 
wt%  CiOn).  These  were  manufactured  using  the  laser  pedestal  growth  (LPG) 
method  by  Lasergenics  Corporation  [36].  The  linear  attenuation  coefficient  of 
this  fiber  was  found  to  be  0.04  cm*^ . 

Figure  3  shows  the  arrangement  used  for  the  composite  layups.  Square 
foil-fiber-foil  layups  were  placed  between  two  0.75  mm  thick,  10.16  cm  diameter 
307  stainless  steel  disks.  The  fibers  were  colinearly  positioned  approximately  5 
mm  iqiart  They  were  then  shipped  to  Kin-Therm.  Inc.  where  they  were 
evacuated  in  a  millitonr  range  vacuum  and  sealed  by  electron  beam  welding.  The 
resulting  can  was  then  placed  in  a  hot  isostatic  press  (HIP)  ( Asea  Brown  Boveii 
Autoclave  Systems,  Inc.,  Minihipper),  heated  to  900**C,  pressurized  at  100  MPa 
for  one  hour,  cooled  at  43*’C/min.  down  to  250'’C,  and  air  cooled  to  room 
temperature  (about  25*'Q.  The  can  ends  were  then  cut  away  to  expose  the 
transverse  section  of  the  fibers.  Some  loss  of  fiber  alignment  was  noted  to  have 
accompanied  the  consolidation  process. 

The  samples  were  mounted  in  epoxy  and  the  exposed  Hber  ends  polished 
to  a  0.05  pm  alumina  finish,  and  examined  by  scanning  electron  microscopy. 

Near  some  of  the  fibers  radial  matrix  cracks  had  formed  on  cooling,  Fi^uic  4.  The 
number  of  these  radial  cracks  varied  from  fiber  to  Hber  (and  from  one  end  of  the 
fiber  to  the  odier)  but  were  present  to  a  greater  or  lesser  extent  on  at  least  one 
end  of  each  of  the  fibers  in  the  composite.  A  ring  of  reaction  product  was  also 


16 


seen  at  the  fiber-matrix  interface.  It  has  been  reported  to  be  Ti3Al  [37].  The 
reaction  product  thickness  varies  with  consolidation  conditions;  for  those  used 
here  it  was  3  pm  thick.  A  debond  crack  had  also  formed  between  the  fiber  and 
the  reaction  product.  This  presumably  formed  (during  sectioning)  in  response  to 
the  tensile  radial  stress  that  has  been  predicted  to  develop  within  ~25  pm  of  the 
freshly  exposed  free  surface  [13, 31]. 

3J2  Spectroscopic  Measurements 

The  spectroscopic  measurements  reported  here  were  conducted  using  a 
fiber  optically  coupled  set-up  shown  in  Hgure  S.  An  Omnichrome  model  532- 
MAP  air-cooled  argon  laser,  set  to  lase  at  514.5  nm,  was  used  to  excite 
fluorescence  in  the  a-Al203  fibers.  The  laser  beam  was  filtered  through  a  CUSO4 

solution  filter  to  eliminate  laser  plasma  lines,  particularly  the  one  at  696 J43  nm. 
The  laser  beam  was  then  focused  onto  one  arm  of  a  bifurcated  fiber  optic  bundle 
(supplied  by  Electro  Hberoptics  Corp.)  by  a  5  cm  focal  length  convex  lens.  The 
(1.6  mm  diameter)  common  end  of  the  fiber  optic  bundle  was  placed  as  close  as 
possible  to  the  end  of  the  ruby  fiber  to  be  interrogated  in  order  to  obtain  the 
maximum  luminescence  signal 

In  practice,  it  was  found  that  the  fib^s  axis  was  sometimes  not 
perpendicular  to  the  polished  composite  surface,  and  so  was  not  always  aligned 
with  the  fibCT  optic  bundle.  This  prevented  the  efficient  coupling  of  the  exciting 
light  into  the  ruby  fiber  and  reduced  the  luminescence  signal  detected.  By 
rotating  the  composite,  it  was  usually  possible  to  align  the  ruby  fiber's  axis  with 
the  axis  of  the  fiber  bundle  which  allowed  for  a  much  improved  signal.  Further 
signal  improvements  were  achieved  by  focusing  the  laser  beam  onto  only  a  single 
fiber  in  the  fiber  optic  bundle  and  positioning  this  over  the  ruby  fiber. 
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The  second  arm  of  the  bifurcated  fiber  optic  bundle  transmitted  the 
luminescence  signal  to  a  0.64  m  single  grating  monochromator  (Instruments  S  A, 
model  HR-640).  Interfacing  the  fiber  bundle  to  the  monochromator  was  achieved 
by  transmitting  the  luminescence  signal  from  the  fiber  bundle  through  20  mm  and 
60  mm  focal  length  lenses  and  then  focusing  upon  the  monochromator  slit.  The 
slit  used  here  was  20  mm  in  height  and  about  1 1  |im  in  width.  The 
monochromator  was  equipped  with  a  1 10  x  1 10  mm  flat  grating  with  1200 
grooves/mm.  The  resolution  at  the  R  line  luminescence  wavelength  using  this 
grating/slit  combination  was  about  0.17  A. 

A  Spectraview-2D  CCD  array  was  attached  to  the  monochromator  as  a 

detector.  The  liquid  nitrogen  cooled  array  had  a  spectral  response  from  4(X)  to 

1000  nm.  The  signal  from  this  array  was  then  transmitted  to  a  personal  computer 

equipped  with  PRISM,  a  collection  of  computer  programs  and  routines  that  have 

been  developed  by  Instruments,  SA,  for  the  purpose  of  acquiring  and  analyzing 

spectroscopic  data  collected  with  their  equipment  The  PRISM  program  was 

used  to  position  the  monochromator  to  the  desired  region  of  the  spectrum,  to 

average  and  smooth  peaks,  and  to  make  a  preliminary  determination  of  peak 

position  wavelength.  For  more  precise  peak  position  determinations,  a  Simplex 

program  was  used  to  fit  a  double  Lorentzian  intensity  curve  to  the  Rl,  R2  line 

spectrum.  The  standard  deviation  of  this  measurement  was  0.1  cm*^ 

(corresponding  to  about  a  50  MPa  stress  error). 

33  Spectroscopic  Experiments 

Spectra  were  obtained  from  each  of  the  three  fibers  in  the  two  composite 
samples.  Two  unstressed  reference  samples  from  the  same  original  fiber  used  in 
each  of  the  composites  were  also  measured  and  used  to  find  the  peak  shifts 
associated  with  compositing.  At  least  three  spectra  were  collected  from  each  end 
of  each  composited  fiber.  The  Simplex  method  was  used  to  determine  the 
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position  of  the  R  peaks,  and  the  position  values  reported  were  then  averaged  for 

each  fiber  end.  The  average  for  each  fib^  end  was  then  subtracted  from  the 
average  for  the  appropriate  unstressed  fiber  to  provide  Av^  The  temperature  was 

recorded  vidth  each  spectrum  taken,  and  the  ai^nopriate  temperature  correction 
made  using  equation  (6). 

4.0  Results 

Hgure  6  shows  two  representative  luminescence  spectra  that  include  both 
R  peaks.  One  corresponds  to  the  signal  from  a  composited  fiber  and  the  other  to 
that  of  an  uncomposited  fiber.  Shifts  due  to  compositing  can  be  seen  for  both 
peaks.  In  order  to  measure  the  R  line  shifts  associated  with  the  compositing 
process,  luminescence  spectra  were  obtained  from  both  ends  of  each  fiber  in  the 
COMP3  and  COMPS  samples  and  from  uncomposited  fibers  from  the  same  fiber 
batches  used  to  produce  the  composites.  The  peak  positions  and  measurement 
temperatures  were  measured  three  times  for  all  fibers  in  both  saiiq>les,  together 
with  the  reference  fibos. 

The  average  R1  and  R2  temperature  corrected  peak  shifts  (Av  i  and  Av2) 

are  given  for  both  ends  of  each  COMPS  sample  fiber  in  Table  4  and  for  the 
COMPS  sample  in  Table  S.  Tstj-  -  Tynstj-  ranged  from  -0.37®C  to  +0.1®C  so  that 

the  temperature  correction  varied  from  •O.OS  to  -fO.Ol  cm~^.  Compositing  can  be 
seen  to  have  resulted  in  a  red  shift  associated  with  the  development  of  a 
predominantly  compressive  stress  in  the  fiber. 
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Table  4;  R  line  peak  shifts  from  the  composite  containing  fibers  doped 
with  0.03wt%  CK®):  (i)  end  A,  and  (ii)  end  B. 

(i)EndA  (ii)EndB 


Fiber 

1 

2 

3 

^3:25 

1 

Fiber 

1 

2 

3 

43:23 

406  1 

Table  S:  R  line  peak  shifts  from  the  composite  containing  fibers  doped  with 
0.05wt%  Ci(ni):  (i)  end  A,  and  (ii)  endB. 

(i)EndA  (ii)EndB 


Fiber 

1 

2 

3 

43:64 

44:86  1 

Fiber 

MM 

^¥iT:ernru 

mam 

3 

Average 

43:87 

43:28  1 

The  errors  in  Tables  4  and  5  include  contributions  associated  with  the 
temperature  correction.  These  arise  from  uncertainty  in  the  temperature 
conection  coefficient  (±0.006  cm'^/K,  shown  in  Table  2),  and  uncertainty  in  the 
two  measurement  temperatures  (±0.S*’Q.  The  tables  have  root  mean  errors 
calculated  using  the  equation: 

Mm  =  {(Av  s  )2  +  (Av +  (Tstt  -  T„TOtt)^Aoik3)^  (10) 

+  (OM) VTstr)^  +  (ATunar)^] ) 

where  A(A^)  is  the  etror  in  the  peak  shift  (in  cm'^),  AO)^  is  the  error  in  the 
temperature  coefficient,  and  ATstr  and  ATynstr  ^  estimated  uncertainties  in 
the  temperature  measurements  [38]. 
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5.0  Discussion 

Large  R1  and  R2  peak  shifts  have  been  observed  to  accompany  the 
incorporation  of  a-Al203  fibers  in  a  CP  titanium  matrix.  They  can  be  compared 

to  predicted  peak  shifts  by  combining  the  expected  state  of  fiber  stress  (from  a 
residual  stress  model)  and  equation  (3)  using  previously  reported  values  for  the 
piezospectroscopic  constants.  We  have  chosen  to  use  the  residual  stress  model 
of  Pindera  and  ^’reed  [9]  to  calculate  the  residual  stress  in  the  fibers  because  it  can 
accorrunodate  fibers  with  either  isotropic,  transversely  isotropic,  or  orthotropic 
elastic  and  CTE  constants  (c-oriented  sapphire  fibers  can  be  reasonably 
approximated  as  transversely  isotropic),  and  because  it  allows  the  matrix  to  be 
either  elastic  or  elastic-plastic.  Because  the  residual  stresses  calcula!^  d  by  the 
Pindera  and  Freed  model  are  based  on  the  local/global  stiffness  matrix  method  of 
Buflo*  [39],  it  facilitates  the  analysis  of  problems  with  many  concentric  cylindrical 
layers.  Thus,  the  effects  of  an  interfacial  reaction  product  (Ti3  Al)  layer  in  the 

samples  used  here  could  also  be  incorporated  into  the  stress  predictions. 

To  calculate  the  stress,  we  assumed  the  fiber  to  be  an  intinite  cylinder  with 
transversely  isotropic,  temperature  dependent  elastic  and  CTE  components. 
Values  for  these  have  been  obtained  from  the  literature  (Table  6).  The  matrix  has 
been  modelled  as  a  concentric  cylinder  with  either  isotropic  elastic,  or  elastic 
plastic  with  a  linear-hardening  constitutive  behavior.  The  expansion  coefficients 
(a).  Young's  moduli  (E),  yield  stresses  (Gy),  and  linear  strain  hardening  parameters 

(SHP)  of  the  matrix  cylinder  were  all  allowed  to  be  dependent  upon  temperature. 
For  one  case,  a  3  pm  thick  matrix  shell  was  included  and  given  the  thermo¬ 
physical  and  mechanical  properties  of  'n3Al  (the  reaction  product).  The  model 
assumed  continuity  of  displacements  at  the  layer  interfaces  (no  debonding)  and  a 
uniform  temperature  throughout  the  body  at  any  moment  during  cooling.  Both 


21 


are  reasonable  for  the  interior  of  the  composite  studied  here  (but  not  for  the  free 
surface  of  the  sectioned  composite  because  of  radial  debonding).  It  computes  the 
axial,  radial,  and  hoop  (=:radial)  fiber  stress  components,  and  those  of  the  matrix  as 
a  function  of  radial  position. 

Table  6:  Physical  properties  used  in  model:  (a)  for  sapphire,  (b)  for  titanium,  and 

(c)forTi3AL 


(a)  Sapphire* 
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(c)Ti3A1** 


T(«C) 

— E7GP35 — 

OyCMPa) 

^4 

r47 

332" 

143 

508 

132 

444 

300 

120 

163 

700 

105 

163 

500 

57 

116 

v»0.3  SHP=3GPa  ot=I0.4x 

Sj^hiie  moduli  calculated  from  [40].  'Htanium  elastic  moduli  and  yield  stress 
taken  ^m  [41].  a  taken  from  (i)  for  titanium  from  [42}  (recommend^  values 
used),  and  (ii)  for  sapphire,  [43].  a'  and  a”  for  sapphire  t^n  from  [44]  (for  a’: 
values  obtained  by  differentiating  the  equation  recommended  for  calculating 
strain  as  a  function  of  temperature;  for  a":  curve  127  used  for  room  temperature 
a-axis  value,  curve  58  used  for  other  a-axis  values;  curve  126  used  for  room 
temperature  c-axis  value,  curve  59  used  for  other  c-axis  values). 

*’^astic  modulus  and  Poisson's  ratio  for  'n3Al  are  from  [45].  Oy  for  'n3Al  at 
T^00°C  is  from  [46].  Remaining  data  are  for  Ti-24A1-1  INb  (at%),  which  has  low 
temperature  properties  similar  to  those  of  II3AI,  taken  from  [47]. 

The  stresses  predicted  by  the  model  are  only  as  valid  as  the  thermophysical 
constants  (Table  6)  that  go  into  the  calculations.  The  elastic  constants,  E33  and 
El  1,  in  Table  6  are  axial  and  transverse  Young's  moduli  for  c-oriented  sapphire 
fibers  whilst  V13  and  V12  are  the  axial  and  transverse  Poisson's  ratios, 
respectively.  They  were  calculated  from  compliance  constant  data  using  a 
transverse  isotropy  approximation  [41]. 

The  residual  stresses  for  several  different  cases,  (see  Fig.  7)  have  been 
analyzed  using  various  assumptions  (presence  or  absence  of  a  reaction  product, 
elastic  or  elastic/plastic  matrix  behavior,  and  variation  of  fiber  volume  fraction). 
Three  different  reported  values  for  the  thermal  expansion  coefficients  of  the 
sapphire  fiber  have  been  examined  since  this  property  dominates  the  predicted 
shifts,  and  there  exists  significant  variability  in  the  reported  values. 
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The  ratio  of  fiber  to  sample  area  is  0.2%  for  the  samples  studied  here,  and 
so  we  have  begun  our  stress  analysis  (Case  I  in  Fig.  7)  by  taking  a  fiber  volume 
fraction  of  0.2%.  The  stress-free  temperature  of  QOO^C  (the  consolidation 
condition)  was  used  and  the  sample  was  allowed  to  cool  to  2S°C.  Transversely 
isotropic  elastic  and  CTE  properties  (the  values  denoted  by  a  in  Table  6)  were 
used  for  the  fiber.  The  CTE  values  were  those  recommended  for  the  Saphikon 
fiber  by  its  producer.  The  matrix  was  assumed  elastic  with  no  reaction  product 
present.  It  was  also  assumed  that  stresses  began  to  develop  immediately  upon 
cooling  (i.e.,  there  was  no  range  of  temperatures  for  which  creep  was  able  to  relax 
the  stress).  The  resulting  stresses  and  their  distribution  are  those  shown  earliei  in 

figure  1.  The  stress  components  in  the  fiber  were  independent  of  position  and 
calculated  to  be  <i,T=<^ee=-194  MPa  and  <j22=-869  MPa.  Substituting  these 

(position  independent)  fiber  stresses  into  equation  (3)  and  using  the  values  of 
Feher  and  Sturge  for  the  piezospectroscopic  constants  [25],  the  predicted  R  line 
shift: 

Av  =(2.7cm'^/GPaX2)(-0.194GPa)+(2.15  cm' VGPaX-0.869  GPa) 

^2.92  cm'l 

This  predicted  line  shift  is  in  the  range  of  the  measured  values  (sec  Tables  4  and  5) 
but  somewhat  smaller  than  the  average  values  -  a  conclusion  that  is  independent 
of  the  choice  of  Ilij  values  in  Table  1. 

The  residual  stress  model  predictions  are  sensitive  to  changes  in  the  fiber 
fraction  and  the  thermophysical  parameters  of  the  model.  Since  there  is  some 
uncertainty  regarding  the  best  choices  for  each  of  these,  the  model  has  been  used 
to  investigate  the  stability  of  the  predicted  shifts  for  the  range  of  possible  values 
for  these  parameters. 

The  first  change  (Case  II)  was  made  to  the  fiber  volume  fraction(Vf).  For 
the  Av  values  calculated  above,  the  fiber  volume  fraction  was  -0.2%  since  this 
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was  the  ratio  of  fiber  to  sample  cross  sectional  areas.  However,  only  three  fibers 
were  (vesent,  and  it  could  be  argued  that  the  central  fibers  may  have  "seen"  a 
different  volume  fraction  to  the  outer  ones.  To  obtain  a  measure  of  the  model's 
sensitivity  to  this,  the  fiber  fraction  was  arbitrarily  increased  to  1.5%.  The  stress 
components  were  (Fig.  7)  similar  to  the  first  case  run,  and  the  peak  shift  was  -2.82 
cm*^  (within  the  error  of  the  experiments,  the  same  as  Case  I).  Thus,  for  the  very 
low  volume  fractions  considered  here,  variability  in  fiber  fraction  can  be 
discounted  as  a  contributor  to  the  disciepancies  between  model-predicted  and 
measured  line  shifts  and  to  the  variability  in  the  line  shifts  between  individual 
fibers. 

The  issue  next  addressed  (Case  m)  was  that  of  matrix  plasticity,  which  can 
relax  the  residual  stresses.  The  residual  stress  calculations  of  Case  I  were 

reperformed,  but  this  time  matrix  plasticity  was  allowed.  The  stress  components 
calculated  were  If  compared 

with  Case  I,  it  is  clear  that  no  significant  plastic  relaxation  occurred,  and  this  can 
also  be  eliminated  as  a  source  of  the  diffmences  in  predicted  and  measured  line 
shift. 

The  next  consideration  (Case  IV)  was  the  observation  that  a  Ti3Al  reaction 
zone  had  formed  around  the  fiber  during  processing  (prior  to  cooling),  as 
observed  in  Figure  4.  The  effect  of  this  can  be  incorporated  by  introducing  a 
concentric  ring  with  the  properties  of  I13AI  (Table  6c)  into  the  model.  When  tl>is 

was  performed,  the  model  again  predicted  similar  fiber  stresses  to  Case  I  and  a  line 
shift  of  -2.93cm' and  we  conclude  that  reaction  products  do  not  significantly 
affect  the  stress. 

The  Hnal  consideration  (Case  V)  was  that  of  uncertainty  in  literature  values 
for  the  CTE.  There  is  a  significant  range  of  CTE  values  reported  in  the  literature, 
and  so  we  have  calculated  the  peak  shift  using  two  other  CTE  values.  For  one 
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case,  the  a!  values  of  Table  6  were  used  [44].  Otherwise,  the  physical  properties 
are  the  same  as  in  Case  I.  These  GTE  values  are  significantly  larger  than  those 
used  in  the  cases  previously  presented,  resulting  in  a  reduced  stress  and  a  peak 
shift  of  only  -1.41cm'^.  This  is  approximately  half  the  value  calculated  for  the 
previous  cases  and  signiHcantly  less  than  any  of  the  measurements.  No  plasticity 
was  predicted  to  have  occurred  in  this  case.  The  lowest  reported  CTE's  were  also 
used  (the  a"  values  are  given  in  Table  6)  to  calculate  the  stress.  The  resulting 
stress  components  were  predictably  much  larger  (Hg.  7)  and  result  in  a  peak  shift 
of  -3.S9  cm~^ .  A  small  amount  of  matrix  plasticity  was  also  predicted  to  have 
occurred,  but  enforcing  an  elastic  matrix  response  in  the  model  did  not 
signiftcantly  change  the  predicted  stress. 

From  these  residual  stress  calculations  we  conclude  that  the  predicted 
peak  shift  lies  in  the  range  of  >1.41  and  >3.59  cm*^,  depending  on  the  values 
taken  for  the  expansion  coefficient  of  sapphire.  Tables  4  and  5  show  that  these 
predicted  values  of  Av  overlap  the  lowor  end  of  the  values  measured  in  COMP3 

and  COMPS.  The  troubling  issues  are  the  significant,  reproducibly  measured 
variability  in  Av  values  between  the  fibers  within  both  composites  and,  for  some 

fibers,  measured  shifts  that  w(;re  significantly  larger  than  any  of  the  predictions. 
The  line  shift  values  observed  here  have  recently  been  independently  veriried  by 
Lipkin  and  Qarke  [28]  and  cannot  be  attributed  to  a  systematic  measurement 
error.  The  variability  is  also  significantly  larger  than  the  random  experimental 
error.  The  calculations  above  show  that  it  cannot  be  accounted  for  by  different 
amounts  of  matrix  plasticity,  fiber  volume  fraction  variations  or  interfacial  reaction 
product  layer  effects. 

The  most  obvious  reason  for  measured  Av  that  are  smaller  (and  more 
variable)  than  predicted  is  the  radial  matrix  cracking  seen  in  the  micrographs  of 
Figure  4.  Highly  variable  amounts  of  crack;  s  observed  for  all  the  fibers  in 
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CX>MP3  and  COMPS.  Since  a  radial  crack  provides  a  mechanism  to  elongate  the 
matrix  in  the  hoop  direction,  this  can  relax  both  the  Abers'  hoop  (and  the 
mechanically  coupled  radial)  stress,  and  thus  reduce  the  peak  shift.  The 
variability  in  cracking  (which  may  extend  along  the  Aber)  would  then  account  for 
the  variability  in  line  shift  observed. 

A  second  possibility  is  that,  for  some  Abers,  the  luminescence  signal  was 
obtained  from  the  end  of  the  Aber  (whm  Aee  surface  effects  result  in  lower 
stresses  than  an  inAnite  Aber  residual  stress  model  would  predict)  due  to 
inefAcient  optical  coupling/transmission  in  the  Aber.  Ma  and  Clarke  have  shown 
that  signiAcant  relaxadon  of  the  stresses  does  occur,  and  within  a  Aber  diameter 
or  so  of  the  Aber  end,  the  radial  stress  can  even  be  tensile  [13, 28].  Lipkin  and 
Qarke  [28]  confirm  that  similar  relaxations  occurred  in  these  samples.  Indeed, 
metallography  of  the  samples  indicated  that  radial  debonding  of  the  Aber>matrix 
intoface  had  occurred  at  the  free  surface,  consistent  with  this  hypothesis. 
Furthermore,  we  note  that  the  COMPS  sample  with  Abers  of  poorer  transmissivity 
gave  slightly  smaller  peak  shift.  Thus,  it  is  possible  that  for  some  (misaligned  or 
highly  attenuating)  Abers,  signals  were  more  heavily  weighted  by  end  effects 
which  contributed  >  mailer  than  expected  peak  shift. 

It  is  not  so  straightforward  to  explain  why  the  measured  Av  values  were 

sometimes  greater  than  the  predicted  values.  One  possibility  could  have  been  a 
variation  in  the  concentration  of  Cr(^)  in  the  Abers.  Ma  and  Clarke  have 
reported  that  changes  in  the  Cr(^)  concentration  creates  different  amounts  of 
strain  in  an  AI2O3  matrix  which  can  be  observed  as  a  shift  in  the  ruby 

luminescence  peaks  [13].  Thus,  enhanced  concentrations  of  Cr(I^)  in  a 
composited  Aber  compared  to  an  unstressed  reference  could  give  an  "anomalous" 
increase  in  the  apparent  compressive  stress.  EDAX  analysis  on  the  Aber  ends  of 
the  COMP3  sample  indicated  this  not  to  be  the  case  here.  The  second  possibility 
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is  that  stress  existed  in  the  pre-composited  fiber.  Though  small  stresses  can 
develop  in  fibers  grown  by  the  EFG  or  LPG  processes,  the  small  temperature 
difieiences  [48]  make  it  unlikely  that  these  are  sufficient  to  account  for  the 
anomalous  shift  observed  here,  especially  since  all  the  fiber  segments  came  from  a 
single  original  piece  of  fiber. 

Further  insight  into  this  anomalous  shift  can  be  gained  by  recalling  that 
equation  (4)  for  the  two  measured  (R1  and  R2)  line  shifts  is  really  two  linear 
equations  in  three  unknowns  (Cn-,  <706  and  and  a  constraining  relation  (On- 
=  t^06)  exists  between  two  of  the  unknowns.  Thus,  in  principle,  one  can 

determine  both  of  the  unknown  components  of  stress.  Using  the  Lipkin  and 
Clarice  [28]  piezospectroscopic  constants  of  Table  1,  and  the  average  line  shifts  of 
Tables  4  and  5.  we  find  the  stress  components  shown  in  Table  7. 

Table  7:  Stress  components  recovered  from  R1  and  R2  line  shift  data. 

(a)  0.03wt%  Cr  doped  fibers 

(i)EndA  (ii)EndB 


Rber 

GzziQPdi) 

aiT=O00(GPa) 

1 

1.14 

2 

IF! 

3 

TIE 

Fiber 

OyyfGPa) 

<yrr=<y00(GPa) 

1 

1.19 

IHi 

1:48 

(b)  0.05wt%  Cr  doped  fibers 
(i)EndA  (u)EndB 


Fiber 

Oy^fGPa) 

OiT=<y00(GPa) 

1 

2 

127 

3 

I5S 

Fiber 

OyyfGPa) 

<yrr=<y00(GPa) 

'■UH 

r:95 

125 

Even  though  the  stress  values  obtained  are  obviously  sensitive  both  to  error  in 
line  shift  measurement  and  in  the  piezospectroscopic  constants,  it  is  clear  from 
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Table  7  that  1)  significantly  lower  than  expected  radial  (and  hoop)  stresses  are 
present  in  these  fibers  and  2)  the  axial  stresses  are  much  greater  than  predicted  by 
a  concentric  cylin^r's  analysis  of  the  fiber-matrix  CTE  mismatch. 

The  lower  than  expected  radial  stress  values  can  be  attributed  to  the  radial 
matrix  cracks,  but  the  high  axial  stress  cannot.  One  explanation  for  the 
anomalous  high  axial  stress  result  is  the  existence  of  a  second  source  of 
compressive  snkmn  the  composite.  Such  a  stress  could  arise  from  a  CTE 
mismatch  with  the  stainless  can  used  for  consolidation  (Hg.  3).  Since  the 
coefficient  of  thermal  expansion  of  the  circular  stainless  steel  canister  used  for 
sample  preparation  (17.3  x  10"^  ®C"^)  was  more  than  twice  that  of  titanium  (8.7  x 
the  entire  composite  would  be  subjected  to  an  additional  compressive 
stress  on  cooling.  At  high  temperatures  during  the  cooling,  the  matrix  would 
creep  or  plastically  deform  in  response  to  this,  and  this  would  transfer  a  significant 
fraction  of  the  can  load  onto  the  fibers  (particularly  in  the  axial  direction). 

To  can  gain  additional  insight  into  this,  we  have  incorporated  into  the 
concentric  cylinder's  analysis,  an  additional  ring  with  the  properties  of  the 
stainless  steel  tooling.  Using  a  volume  fraction  of  50%  for  the  can  and  relaxing 
the  pressure  used  for  consolidation  monotonically,  the  concentric  cylinder 
analysis  predicts  an  axial  (compressive)  fiber  stress  of  -2730  MPa  and  a  radial 
stress  of  -280  MPa  after  cooling.  Some  of  this  would  be  relaxed  when  the  can 
material  is  removed,  but  the  inelastic  matrix  deformations  incurred  during  cooling 
ensure  that  an  additional  compressive  fiber  stress  would  be  present,  and  a  larger 
than  otherwise  predicted  peak  shift  for  both  R  lines  occurs. 

We  propose  that  the  stress  in  the  fibers  of  HIP  consolidated  composites 
arises  from  a  combination  of  fiber/matrix  and  compositeAooling  CTE  mismatch  in 
conjunction  with  a  relaxation  of  the  transverse  stress  components  by  radial  matrix 
cracking.  It  is  interesting  to  note  that  the  large  measured  compressive  stress  that 
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remains  in  the  fibers  studied  here  would  have  a  significant  effect  upon  the 
apparent  insitu  strength  of  the  fiber  and  thus  the  composite's  strength  in  tension. 
Compressive  strength  would  also  be  adversely  impacted  due  to  fiber  buckling  at 
lower  applied  loads.  The  additional  stress  would  also  affect  the  extent  of  fiber- 
matrix  debonding,  and  thus  pullout  behavior,  during  crack  growth.  These  tooling 
associated  stresses  are  ignored  in  the  modelling  of  mechanical  properties,  and  this 
needs  to  be  addressed,  since  mechanical  properties  are  likely  to  depend  upon  the 
material  used  for  HIP  tooling.  In  the  extreme  case  where  tooling  materials  are 
used  whose  expansion  coefficients  are  much  less  than  that  of  the  composite,  the 
fiber  could  even  suffer  a  tensile  axial  residual  stress  with  the  risk  of  fiber 
fragmentation  during  cooling  and  a  resulting  composite  of  very  poor  quality. 

6.0  Summary 

The  underlying  principles  of  a  luminescence  method  for  measuring  stress 
have  been  reviewed  and  applied  to  the  problem  of  residual  stress  determination  in 
metal  matrix  composites.  It  has  been  shown  that  the  change  in  residual  stress 
associated  with  the  compositing  of  Cr  doped  a-Al203  ribers  in  a  titanium  matrix 

can  be  detected  via  a  fiber  optically  coupled  luminescence  measurement.  The 
luminescence  peak  shifts  are  smaller  in  magnitude  than  those  predicted  from 
thermal  stress  models,  and  a  reproducible  variability  in  the  line  shifts  for  di^erent 
fibers  in  the  composite  has  been  observed.  Using  both  R1  and  R2  line  data,  the 
axial  and  transverse  stress  contributions  to  the  line  shifts  have  been  separated,  we 
find  the  axial  stresses  to  be  more  than  two  times  that  of  the  predictions  whilst 
those  in  the  radial  (and  hoop)  directions  are  signiricantly  smaller.  We  attribute  the 
large  axial  stress  to  a  CTE  difference  between  the  HIP  containment  and  the 
composite  which  is  not  included  in  most  residual  stress  analyses.  The  low 
measured  radial/hoop  fiber  stress  components  are  proposed  to  result  from  radial 
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matrix  cracking.  These  two  "unmodeiled"  phenomena  have  offset  one  another  in 
the  R1  or  R2  luminescence  measurement  and  resulted  in  the  fortuitous  rough 
agreement  between  concentric  model  predictions  and  individual  peak  shift 
measurements.  The  stress  developed  because  of  the  HIP  tooling  could 
significantly  affect  the  quality  of  metal  matrix  composites  and  introduces  an 
additional  dependence  of  mechanical  properties  upon  the  detailed  manner  of 
composite  manufacture. 
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ELASTIC-PLASTIC  MODEL  PREDICTED  STRESSES 
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Figure  1,  Schematic  of  Pindera-Freed  model  predictions 
components  in  and  near  a  composited  fiber. 
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Figure  2.  Diagram  showing  the  electronic  structure  of  ruby  for  those 

transitions  in  the  visible  region  and  the  phenomena  that  cause  the 
splittings.  (Redrawn  from  [15]). 
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Figure  3,  Schematic  representation  of  spatial  arrangements  of  Ti  foil  and 
ruby  fibers  prior  to  consolidation.  (Square  stainless  steel  sheets 
are  shown  for  simplicity:  circular  disks  of  stainless  steel  were  used 
in  practice.) 


Figfure  4.  Scanning  electron  micrograph  of  the  fibers  in  COMP3:  (a)  the 
middle  fiber,  and  (b)  closeup  of  the  middle  fiber. 
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Figure  5.  Schematic  diagram  of  the  spectroscopic  system  used  to  fiber 
optically  measure  the  residual  stress  in  MMC's. 
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Figure  6.  Graph  of  peak  intensity  versus  wavelength  showing  a  typical  peak 
shift  associated  with  compositing. 
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Figure  7.  Model  predicted  R1  line  shifts  for  various  sapphire/CP  titanium 
model  systems 
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OPTICAL  FLUORESCENCE  FROM  CHROMIUM  IONS 
IN  SAPPHIRE;  A  PROBE  OF  THE  IMAGE  STRESS 

QING  MA  aad  D.  R.  CLARKE 

Materwli  Deputment,  University  of  California.  Santa  Barbara.  CA  93106.  U.S.A. 

(Rectwtd  7  October  1992) 

AftMnci— The  concentration  dependent  frequency  shiA  of  the  fluorescence  from  chromium  ions  in 
sapphire  is  shown  to  be  a  direct  probe  of  the  image  stress  produced  by  the  chromium  substituting  for 
the  smaller  aluminium  ions  in  sapphire.  Our  analysis  is  based  on  a  comparison  of  the  frequency  shift 
produced  when  the  volume  of  a  sapphire  crystal  is  changed  by  an  external  pressure  to  that  produced  by 
a  concentration  of  substitutional  solute.  The  volumetric  strain  sensed  by  the  fluorescing  chromium  ion 
when  the  lattice  is  expanded  by  substitutional  solute  is  approximately  one-third  of  that  determined  from 
X-ray  lattioe  parameter  measurements,  in  accord  with  the  Eshelby  analysis  of  the  strain  held  due  to  a 


concentration  of  elastic  defects  in  a  flnite  body. 


1.  INTRODUCTION 

The  question  as  to  the  lattice  distortions  produced  in 
a  crystal  by  a  substitutional  solute  is  of  long-standing 
interest,  and  has  been  the  subject  of  considerable 
intellectual  inquiry.  Although  it  has  been  recognized 
throughout  that  the  only  proper  description  is  one 
that  takes  account  of  the  atomic  nature  of  the  solute 
and  the  appropriate  inter-atomic  potentials,  early 
investigators  did  not  have  available  to  them  this 
information.  In  the  absence  of  such  information,  it 
was  necessary  to  address  the  question  using  contin¬ 
uum  descriptions  of  the  solute  and  its  effect  on  the 
lattioe.  One  approach,  that  might  be  termed  the 
elastic  defect  model,  was  to  consider  a  solute  as  a 
center  of  dilatation  supplemented  by  point  singular¬ 
ities  of  higher  order.  This  enabled  tte  methods  of 
continuum  elasticity  pertaining  to  centers  of  pressure 
to  be  employed  in  finding  solutions.  For  'Example, 
Bilby  modeM  an  interstitial  atom  as  a  misfitting 
elastic  sphere  forced  into  a  smaller  hole  in  a  material, 
and  derived  solutions  for  the  interaction  of  a  solute 
with  a  dislocation  [1].  Later,  Eshelby,  in  one  of  the 
classic  papers  in  materials,  provided  a  rigorous  foun¬ 
dation  for  the  elastic  defect  model  by  deriving  the 
forces  acting  on  an  elastic  singularity  in  a  solid  [2].  In 
subsequent  seminal  papers  [3-5],  E^lby  set  forth  a 
comprehensive  framework  for  the  computation  of 
elastic  displacements,  strains  and  stresses  resulting 
from  the  presence  of  point  defects,  inclusions,  precipi¬ 
tates  and  transforming  regions.  By  envisaging  a  series 
of  hypothetical  cutting,  transforming  and  welding 
operations,  and  the  forces  and  displacements  in¬ 
volved,  Eshelby  formulated  analytical  expressions  for 
the  stresses  and  strains  created  by  a  region  transform¬ 
ing  but  constrained  by  a  surrounding  matrix. 

Two  concepts  introduced  by  Eshelby  are  of  par¬ 
ticular  importance.  One  is  that  inside  any  finite  body 


containing  defects  there  exists  an  additional  displace¬ 
ment  field,  and  by  extension  strain  and  stress  fkid, 
that  arises  from  the  boundary  condition  that  at  the 
surface  of  the  body  the  normal  stress  must  be  zero. 
By  analogy  with  electromagnetic  fields,  Eshelby 
termed  this  additional  elastic  field  the  image  field.  The 
image  field  is  independent  of  position  inside  the  body 
and  independent  of  the  size  of  the  body.  The  second 
key  concept  was  that  of  the  “stress-free”  strain,  the 
shape  change  of  the  transforming  region  that  would 
have  occurred  in  the  absence  of  any  constraint.  (This 
same  strain  has  been  denoted  as  the  “eigenstrain"  by 
Mura  (6).  It  also  appears  as  the  “natural”  strain  in 
recent  literature  on  the  virtual  crystal  approximation 
(VGA)  approach  to  elearonic  structure  calculations 
[7].)  Using  the  concept  of  a  stress-free  transformation 
strain,  Eshelby  was  able  to  derive  general  expressions 
for  the  strains  produced  in  the  matrix  as  a  result  of 
the  traiuformation  as  a  function  of  the  shape  of  the 
ellipsoidal  inclusion.  The  matrix  strain  is  shown  to 
consist  of  two  components,  the  strain  in  the  matrix  if 
it  were  of  infinite  extent  and  the  image  strain. 

The  existence  of  two  components  to  the  strain 
field  in  a  body  containing  defeas  naturally  raises  the 
question  as  to  which  strains  are  actually  being  deter¬ 
mined  when  measuring  the  lattice  parameter  '  / 
diffraction  methods.  Although  there  was.  for  a 
period,  considerable  confusion  concerning  this  ques¬ 
tion,  it  now  appears  resolved  that  Bragg  diffraction 
methods,  using  X-rays  and  neutrons,  provide  a  direct 
measure  of  the  spatially  averaged  total  strain  in  the 
body.  (Diffuse  scattering  provides  additional  infor- 
nution  on  the  distribution  of  the  lattice  distortions 
but  does  not  include  information  about  average 
lattice  strain.)  Since  the  image  strain  was  shown  by 
Eshelby  to  be  a  constant  fraaion  of  the  total  strain, 
the  magnitude  of  the  image  strain  has  been  deduced 
from  tte  concentration  dependence  of  the  lattice 
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panuneter.  However,  there  have  not  been,  as  far  as 
the  authors  are  aware,  any  direct  measurements  of 
the  image  strain. 

In  the  course  of  rather  different  studies  (8-10), 
directed  towards  the  application  of  optical  fluor¬ 
escence  methods  to  the  measurement  of  stresses  in 
composite  materials,  we  have,  we  believe,  found  that 
the  fluoresence  from  substitutional  solutes  of  chro¬ 
mium  directly  probes  the  image  stress  in  single  crys¬ 
tals  of  sapphire  (ALO,).  In  this  short  paper,  we 
describe  how  we  reached  this  conclusion.  It  is  first 
necessary,  however,  to  re-examine  the  data  in  the 
literature  on  the  chromium  solid  solution  in  alumi¬ 
num  oxide,  together  with  our  own  experimental  dau. 

2.  CHROMIUM  SUBSTITUTION  IN  SAPPHIRE 

Solid  solutions  of  chromium  in  aluminum  oxide 
provide  a  rich  data  base  to  re-examine  the  question 
of  the  short-range  and  long-range  distortion  induced 
by  the  substitution  of  a  mis-fitting  solute.  In  the 
literature  there  exist  data  on  the  chromium  concen¬ 
tration  dependence  of  density  [11],  lattice  parameter 
determined  by  X-ray  diffraction  (1 1],  and  frequency 
of  characteristic  optical  fluorescence  lines  (12).  In 
addition,  both  lattice  parameter  and  fluorescence  line 
shifts  as  a  function  of  hydrostatic  pressure  have  been 
reported  for  selected  chromium  ion  concentrations 
(13-15).  The  system  AljOj-CrjOs  also  has  the  advan¬ 
tage  that  there  is  complete  solid  solution  between  the 
two  oxide  phases  [16]. 

The  substitutional  nature  of  the  chromium  ion  has 
been  confirmed  by  density  and  X-ray  measurements 
(below).  This  is  consistent  with  both  the  existence  of 
the  continuous  solid  solution  with  CrjO],  the  simi¬ 
larity  in  ionk  radius  between  Cr^*  (O.OMnm)  and 
Al^'*'  (0.0S7  run),  and  the  fact  that  the  corundum 
structure  is  a  rather  close-packed  structure  and  does 
not  contain  any  interstices  as  large  as  the  chromium 
ion.  Structure  analysis  as  a  function  of  hydrostatic 
pressure  further  shows  that  the  corundum  structure 
of  sapphire  compresses  uniformly,  with  aU  the  dis¬ 
tances  decreasing  equally  and  all  the  interatomk 
angles  remaining  unchanged  (17). 

Chromium  in  sapphire  is  also  a  chromophote, 
giving  rise  to  the  characteristk  red  fluroresence  of 
ruby.  The  fluorescence  originates  from  electronk 
transitions  in  the  chromium  ion  produced  by  appro¬ 
priate  excitation.  Electrons  excited  from  their  ground 
state,  *jl2,  can  lose  their  excess  energy  by  a  variety  of 
radiative  (and  non-radiative)  processes.  The  most 
pronounced  transitions  being  the  ff-fluorescence  lines 
at  1.790  and  l.794eV  in  the  red.  (A  recent  paper 
tabulates  all  the  radiative  transitions  in  the  visible 
part  of  the  spectrum  [18].)  The  diromium  ion  can 
substitute  for  any  of  the  aluminum  ions  in  the 
corundum  crystal  structure  of  sapphire  with  the  result 
that  they  are  octahedrally  coordinated  to  neighboring 
oxygen  and  aluminum  ions.  The  octahedron  is  not 
regular,  but  elongated  along  the  c-axis  of  the  lattice. 


resulting  in  splitting  of  the  R -fluorescence  lines  to 
form  the  doublet  of  R,  and  R,  lines.  X-ray  refine¬ 
ments  [19, 20],  in  highly  doped  (4  at.%  Cr)  sapphire 
show  the  chromium  ion  to  be  slightly  displaced  along 
the  r-axis  from  the  position  normally  oa:upied  by  the 
'uminum  ion  so  that  the  chromium  site  symmetry  is 
.irther  reduced,  but  nevertheless  retains  a  three  fold 
rotational  symmetry  about  the  c-axis.  The  effect  of  an 
applied  stress  is  to  distort  the  crystal  fkid  surround¬ 
ing  the  chromium  ion,  change  its  potential  energy  and 
thereby  alter  the  energies  of  the  radiative  transitions. 

2./.  Optical  fluorescence  data 

The  frequencies  of  the  R,  and  R^,  as  well  as  the 
extended  neighbor  fluorescence  N  lines,  have  all  been 
reported  to  shift  with  pressure  and  with  chromium 
concentration.  The  pressure  dependence  is  well 
known  [13, 14]  and  continues  to  be  subject  to  further 
refinement  broause  of  iu  central  importance  as  a 
pressure  gauge  for  high  pressure  studies.  The  re¬ 
ported  pressure  dependence  of  the  frequency  shift  Av 
may  be  expressed  as 

Av  =  -(2n„-t-n„)R  (1) 

where  11^  is  the  piezo-spectroscopic  coefficient.  The 
original  work  of  Forman  et  al.  [13]  is  reproduced  in 
Fig.  1(a)  for  the  R,  line  with  the  sum  211,,  -f  flj, 
having  a  value  of  7.7  cm" '  GPa" '.  The  results  for  Rj 
line  is  very  similar  and  therefore  not  included  here. 
Recent  work,  examining  the  fluoresence  from  beams 
loaded  in  four  point  bending,  has  shown  that  the 
equation  is  equally  valid  for  tension  [8, 9]. 

The  chromium  concentration  dependeixe  of  the 
fluorescence  energy  has  been  the  subject  of  just  one 
previous  study  [12].  Our  own  results  have  proven  to 
be  essentially  identical  to  those  of  Kaplyanskii  et  al. 
[12],  whkh  are  reproduced  in  Fig.  1(b).  The  data  may 
be  expressed  by  the  linear  dependence 

Av  =  99c,  (2) 

where  Av  is  measured  in  wavenumber  and  c,  is  the 
chromium  concentration  (in  wt%). 

2.2.  X-ray  lattice  parameter  rrteasuremenu 

X-ray  lattice  parameter  measurements  have  been 
made  as  a  function  of  pressure  on  rubies  with  a 
number  of  different  chromium  concentrations 
[IS,  17, 18).  These  studies  indicate  an  isotropic  com¬ 
pressibility  with  the  interatomk  distances  decreasing 
equally  with  pressure  and  the  interatomk  angles 
remaining  unchanged  up  to  at  least  4.6  GPa  [IT].  The 
results  of  d' Amour  et  al.  [IS]  is  shown  in  Fig.  1(c), 
whkh  can  be  fitted  to  the  following  linear  relations 

a  «  0.47664(1  -  0.001 173R) 
and 


c«  1.3007(1  -0.001332R) 


(3) 
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where  the  lattice  parameters  a  and  c  are  measured  in 
nanometers,  and  P  is  measured  in  GPa.  They  also 
found  that  the  density  increased  with  increasing 
pressure  according  to  the  Murnaghan  equation  of 
sute. 

The  X*ray  lattice  parameter  measurements  were 
also  made  as  a  function  of  chromium  concentration 
[II]  [Fig.  1(d)],  over  a  range  of  chromium  concen¬ 
tration  up  to  2  wt%  by  Jan  el  al.  They  indicate  that 
the  lattice  parameters  and  derived  density  are  given 
by  the  lineiu’  relationship 

0-0.47591(1 0.0527c.) 

cm  1.29894(1  +  0.0452c.) 


and 

P;r  =  3.9864(1  -I- 0.340c.)  (4) 

where  the  density  px  is  measured  in  g/cm\ 

2.3.  Density  data 

Jan  et  al.  have  also  made  precision  density 
measurements  [11]  using  the  picnometry  technique. 
Their  results  indicate  a  linear  dependence  of  the 
density  on  the  chromium  concentration  given  by  the 
relationship 

p  =3.9860(1  -t- 0.341c.).  (5) 

The  linear  dependence,  as  clearly  shown  in  Fig.  2.  is 
consistent  with  chromium  being  a  substitutional  solute 


1814 


QING  MA  and  CLARKE;  OPTICAL  FLUORESCENCE  IN  SAPPHIRE 


Fig.  2.  Linear  dependence  of  Uie  density  of  alumina  oxide 
on  the  chromium  concentration  (1 1]. 


[21].  The  physical  density  determined  above  is,  within 
experimental  error,  the  same  as  the  X-ray  density 
measurements  by  Jan  et  al.  in  equation  (4).  This 
demonstrates  the  equivalence  of  the  solute  induced 
macroscopic  strain  and  the  mean  strain  measured  by 
X-ray. 

2.4.  Frequency  shift  as  a  function  of  mean  strain 

From  the  results  presented  in  Fig.  1,  the  frequency 
shift  can  be  related  to  the  mean  volumetric  strain 
AF/K  by  using  the  following  relationship 

AF  „An  Ac 

—  -2 —  H - . 

F  a  c 

Therefore,  from  equations  (1)  and  (3),  and  equations 
(2)  and  (4),  we  obtain  the  following  expressions  for 
the  fluorescence  shifts 

Av'-2100^ 

AvO-660^  (6) 

depending  on  whether  the  strain  is  produced  by  an 
external  pressure  or  by  solute  doping.  The  super¬ 
scripts  P  and  D  denote  pressure  and  solute  doping 
respectively. 


3.  ANALYSIS  AND  DISCUSSION 

The  intriguing  feature  of  the  data  assembled  in  the 
previous  section  is  that  the  apparent  strain  sensed  by 
a  fluorescing  chromium  ion  is  difierent  depending  on 
whether  the  lattice  is  strained  by  the  application  of  a 
hydrostatic  pressure  or  internally  by  ^e  addition  of 
substitutional  solute.  Specifically,  for  the  same  lattice 
strain  measured  by  changes  in  the  lattice  parameter 
by  X-ray  diffraction  the  frequency  shift  of  the  chro¬ 
mium  line  produced  by  doping  is  about  one-third 
[equation  (6)]  that  produced  by  an  applied  pressure 

Av®~iAv'  (7) 

We  believe  that  the  different  pressure  and  solute 
concentration  effects  may  be  reconciled  by  postulat¬ 


ing  that  the  frequency  shift  of  the  fluorescence  lines 
is  sensitive  to  the  image  strain  whereas  X-ray  diffrac¬ 
tion  monitors  the  average  lattice  strain.  To  prove  this, 
we  proceed  to  perform  a  set  of  imaginary  experiments 
(Fig.  3),  assuming  for  simplicity  that  the  medium  is 
elastically  isotropic  and  that  the  concentration  of 
solute  is  sufficiently  small  that  any  solute  interaction 
effects,  such  as  dipolar  or  higher,  can  be  ignored. 

We  start  with  a  body  of  sapphire  of  volume  F, 
free  of  external  pressure  (/*  •  0),  containing  only  a 
few  Cr^^  ions  to  sense  the  strain,  so  that  the  solute 
concentration / -*0.  The  frequency  of  the  fluorescence 
line  is 

V  «  v(/»  =  0,/-0).  (8) 

This  frequency  is  expected  to  be,  in  a  complicated 
manner,  dependent  on  the  local  strain  in  the  vicinity 
of  each  Cr^^  ion  as  it  is  a  substitutional  solute,  larger 
than  the  Al’^  ion  it  replaces.  Fortunately,  by  choos¬ 
ing  this  state  as  the  reference  state,  a  knowledge  of  the 
effect  of  this  local  strain  is  unimportant  since  it  is  the 
same  for  both  the  state  of  interest  and  the  reference 
state  and  hence  cancels  out. 

In  the  first  experiment,  a  hydrostatic  tensile  stress 
— /*  is  applied  to  the  body,  so  that  its  volume 
increases  by  AF.  The  strain  dependence  (due  to  an 
applied  stress  —P)  of  the  frequency  shift  [9]  for  an 
isotropic  material  can  be  in  general  expressed  as 

Av»A£,  (9) 

where  the  convention  of  summation  of  repeated 
indices  is  assumed  and  where  X  is  the  piezo-spectro- 


v+av,-p.f-o  v+av.P-o.f-N/v 

Fig.  3.  Illustration  of  a  set  of  imaginary  experiments  to 
produce  a  volume  change  by  either  external  stress  or  solute 
doping.  The  region  W  is  small  but  finite  volume  of  arbitrary 
size  that  contains  a  single  probe  ion  used  to  monitor  the 
strain. 


QING  MA  and  CLARKE:  OPTICAL  FLUORESCENCE  IN  SAPPHIRE 


ISIS 


icopic  volumetric  strain  coefficient.  Therefore,  the 
frequency  shift  due  to  the  stress  -  P  is  simply 

Av'  -  v(  -  P,/-0)  -  v(P  =  0,/-*0) 


(10) 


Alternatively,  in  the  second  experiment,  we  can 
produce  such  a  volume  change  by  increasing  the 
number  of  Cr^*  ions  in  the  body  from  virtually  zero 
in  our  reference  state  to  N,  distributed  at  random. 
Following  Eshelby,  the  displacement  field  generated 
in  an  infinite  body  by  all  the  substitutional  ions  may 
be  written  as 


■*(r)  = 


C 


I 


Ir-M’ 


(11) 


where  C  is  the  dilational  strength  of  an  individual 
CH*  ion,  proportional  to  its  misfit. 

In  a  finite  crystal,  the  displacements  are  modified 
by  the  existence  of  the  image  term,  required  so  that 
the  tractions  at  the  external  surface  are  zero.  The 
additional  displacement  at  position  r  produced  by  the 
image  force  is  [3] 


u“(r)  = 


8«(l-2<r)C)V_ 
3(1 +<r)  y  ' 


(12) 


where  a  is  Poisson’s  ratio.  The  total  displacement  is 
then  the  sum  of  these  two  displacements 


o“^(r)  =  C 


8n(l  -2ff)CAr_ 
3(1+0)  V  ’’ 


(13) 


The  volume  change  produced  by  the  solute  is  then 
simply  the  integral  of  the  displacement  field  over  the 
surface  of  the  body 


AF=J  o'«*'(r)  ds 

=  4„C)vilLl^  =  4nCNy.  (14) 
1+0  ' 


From  this  volume  change,  the  mean  strain  in  the 
body  produced  by  solute  doping  is 

«/.  =  ^  =  4jivC  p  =  4jty  Q'  (1 5) 


image  term  due  to  all  substitutional  ions  (equation 
(12)],  and  the  displacement  field  due  to  the  probing 
ion  itself,  since  the  displacements  due  to  other  ions  in 
equation  (II)  not  included  in  If' do  not  contribute  to 
the  integration.  Because  the  displacements  generated 
by  the  probing  ion  are  already  in  the  reference  state, 
to  which  the  volume  change  is  being  compared,  its 
contribution  is  subtracted  by  the  second  term  in  the 
equation  (16).  So,  effectively,  only  the  image  displace¬ 
ment  field  due  to  the  presence  of  all  substitutional 
ions  in  tlie  body  contributes  to  AH'.  We  may  then 
write 


AIF* 


8ii(l  -2g) 
1  +ff 


iVCf. 


(17) 


Therefore,  the  strain  sensed  by  the  solute  ion  is  simply 
the  image  strain 

c“"  =  ^  =  4x(y-l)C/  (18) 


Comparing  with  equation  (15),  the  corresponding 
frequency  shift  of  the  fluorescence  from  the  sensing 
solute  ion  would  then  be 


Av“  =  v^P  =  0,/=  -  v(P  =  0,/- 

V  -  1  . V  -  1 AK 

=  =  - — . 


•0) 


(19) 


If  one  assumes  that  the  Poisson's  ratio  is  1/3,  then  a 
comparison  of  equations  (10)  and  (19)  indicates  that 
the  effect  of  solute  concentrabon  on  the  fluorescence 
shift  is  about  1/3  that  produced  by  an  external 
pressure,  as  found  experimentally  (equation  (6)}. 
Since  in  the  case  of  solute  concentration,  the  strain 
affecting  the  shift  is  entirely  due  to  the  image  term, 
we  therefore  conclude  that  the  fluorescence  probes 
the  image  strain  whereas  X-ray  diffraction  provides  a 
measurement  of  the  total  strain. 
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Abstract 

In  cooling  sintered  co^^x)sites  of  Al203-Zr02  from  their  fabrication  temperature  residual 
stresses  are  created  as  a  result  of  both  the  difference  in  thermal  expansion  between  the  two 
phases  and  the  aystallographically  anisotropic  thermal  expansion  of  the  AI2O3  phase.  In  this 
work  the  first  and  second  moments  of  the  residual  stress  distribution  have  been  measured  as  a 
function  of  volume  fraction  of  zirconia  from  0.01  to  0.90.  The  measurement  technique  used  is 
piezo-spectroscof^  based  on  the  optical  fluoresceiKe  from  Cr^  dopants  in  the  alumina  phase. 
For  zirconia  volume  fractions  up  to  035  the  averse  stress  accurately  fits  the  pn-edictions  given  by 
the  upper  Hashin  bound  and  this  fit  provides  a  value  of  the  average  thermal  strain  in  the 
composites.  Using  this  value,  the  effective  medium  approximation  produces  an  excellent 
description  of  the  average  stress  over  the  entire  volume  fraction.  It  is  also  shown  that  the 
fluorescence  broadening  due  to  stress  fluctuations  lies  close  to  the  predicted  upper  and  lower 
Hashin  bounds  modified  by  the  restrictior^s  imposed  by  the  principle  of  maximum  entropy.  The 
measured  moments  and  those  predicted  by  stochastic  stress  analysis  compare  well  suggesting 
that  the  stochastic  analysis  provides  a  reliable  method  of  calculating  residual  stresses  in 
composites. 


1  Introduction 


The  determination  of  the  average  elastic  properties  of  composite  and  poly  crystalline  materials 
is  one  of  the  classical  problems  in  the  mechanics  of  materials.  Over  the  years  a  variety  of  bounds 
on  the  elastic  modulus  have  been  established,  ranging  from  the  early  Voigt  and  Reuss  bounds  [1], 
to  the  Hashin-Strikman  bounds  [2]  and  more  recently  a  variety  of  generalized  effective  medium 
equations  [3].  Whilst  much  of  the  attention  has  been  focused  on  predicting  the  elastic  modulus  of 
two-  (or  multi-)  phase  materials  from  the  elastic  properties  of  the  constitutive  phases,  the  models 
also  allow,  in  principle,  the  calculation  of  the  stress  distribution  in  the  phases  [3,4,5].  However, 
since  the  calculation  of  residual  stresses  can  be  complex,  three  quite  distinct  approaches  have 
been  adopted  depending  on  the  type  of  information  about  the  residual  stresses  desired.  One  is 
the  approach  embodied  in  the  Eshelby  calculation  of  the  strun  field  due  to  a  misfitting  inclusion 
[4,6].  Once  a  stress  free  transformation  strain  and  the  shape  of  the  inclusion  is  specified,  the  stress 
field  throughout  the  material  can  be  calculated  in  detail.  The  analysis  can  then  be  extended  to 
a  composite  containing  other  inclusions  by  the  linear  superposition  of  stresses  using  assumptions 
about  the  distribution  of  the  inclusions  [4].  The  second  approach  is  to  presume  some  particular 
geometry  of  the  inclusions  and  their  arrangement  and  compute  the  detailed  stress  distribution  by 
application  of  large-scale  finite  element  methods  [7].  The  third  approach  is  to  describe  the  spatial 
distribution  of  stresses  by  statistical  methods  [3].  The  stress  field  is  described  by  an  appropriate 
probability  functional  from  which  the  average  value  of  the  stresses  (in  the  various  phases)  are 
expressed  in  terms  of  the  first  moments  of  the  stress  distribution.  Similarly,  the  spatially  varying 
stresses  are  expressed  in  terms  of  the  second  moment  of  the  distribution.  In  the  framework  of 
these  stochastic  models,  the  morphology  of  the  microstructure  and  the  spatial  variation  of  phases 
are  described  in  terms  of  appropriate  correlation  functions.  By  applying  the  maximum  entrop> 
principle,  further  restrictions  can  be  placed  on  the  properties  of  the  statistical  distribution  ami 
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hence  on  the  first  and  second  moments  of  the  stress  distribution. 

Although  it  might  .be  expected  that  such  a  statistically-based  approach  might  be  most  appro¬ 
priate  in  describing  the  average  stresses  in  a  polycrystalline,  two-  (or  multi-)  phase  material  the 
predictions  have  yet  to  be  critically  tested  by  experiment.  In  this  contribution  we  attempt  to 
provide  such  a  test  using  measurements  of  the  optical  fluorescence  shifts  from  the  alumina  phase 
(A)  in  a  series  of  Al203-Zr03  composites  in  which  the  volume  fraction  of  the  cubic  zirconia  phase 
(Z)  is  varied.  The  composites  were  prepared  by  sintering  of  alumina  and  cubic  zirconia  particles 
and  as  a  result  the  spatial  distribution  of  grains  is  essentially  disordered  (figure  1).  In  this  ma- 
teri2d  system  residuaJ  stresses  form  on  cooling,  due  to  both  the  thermal  mismatch  between  the 
alumina  and  zirconia  phases  and  the  crystallographic  thermal  contraction  anisotropy  of  the  alu¬ 
mina  material.  Since  both  phases  have  very  high  yield  stresses,  the  residual  stresses  are  believed 
to  be  accommodated  elastically  to  a  large  extent.  Furth''rmore,  the  }rttrium  doping  of  the  zirconia 
ensures  that  it  remains  cubic  and  no  phase  transformation  can  occur  on  cooling. 

The  technique  of  optical  fluorescence  is  used  to  determine  the  stress  distribution  in  the  alumina 
grains  in  the  composite.  The  particular  fluorescence  used  is  the  R2  fluorescence  line  from  electronic 
transitions  of  Cr^'*'  substitutional  dopants  in  the  alumina,  excited  by  an  optical  probe.  The 
theoretical  basis  of  the  fluorescence  technique  has  been  described  in  detail  in  an  earlier  publication 
[8].  Briefly,  the  frequency  shift  of  a  characteristic  fluorescence  line  is  stress  dependent  and  is  given 
by: 

Ai/  =  n.j  <7*  ( 1 ) 

where  11, ^  are  the  piezo-spectroscopic  coefficients  relating  the  frequency  shift  to  the  stress,  and 
<7'  is  the  stress  tensor  represented  in  the  crystallographic  coordinate  frame.  .As  shown  earlier, 
for  a  polycrystalline  material,  the  shift  in  fluorescence  frequency  <  Ai/  >  is  a  measure  of  the 
first  moment  of  the  stress  distribution  and  the  broadening  of  the  fluorescence  peak  \/<  Ai/^  >  i< 
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dependent  on  the  second  moment  of  the  stress  distribution.  By  using  a  probe  large  compared  to 
the  grain  size  and  by  ptiiizing  the  extensive  multiple  scattering  of  the  light  within  the  composite, 
the  fluorescence  signal  is  obtained  from  a  very  large  number  of  grains.  Thus,  by  measuring  the 
fluorescence  peak  the  first  and  second  moments  of  the  stress  distribution  from  the  material  can 
be  obtained  directly  and  compared  with  the  predictions  of  the  stochastic  stress  analysis. 

2  Stochastic  Analysis  of  Residual  Stresses  In  Two-Phase 
Materiab 

In  this  section  we  outline  the  principal  results  of  the  stochaistic  analysis  of  the  stresses  in  two- 
phase  materials  appropriate  to  a  composite  of  Al^Oa-ZrOj  [3].  The  underlying  assumptions  being 
that  the  stress  distribution  in  the  composite  can  be  described  by  a  probability  functional  and 
that  the  spatial  and  angular  distributions  of  the  zirconia  and  alumina  grains  are  random.  The 
emphasis  is  on  the  first  and  second  moments  of  the  stress  distribution  since  they  are  those  that 
are  measured  directly  by  the  fluorescence  technique. 

To  apply  the  stochastic  stress  analysis,  a  series  of  . restrictions  can  be  imposed  so  as  to  succes¬ 
sively  refine  the  predictions.  The  first  step  is  a  general  statement  of  the  first  and  second  moments 
of  the  stress  distribution  in  terms  of  the  elastic  strain  energy,  u,  stored  in  the  stress  field.  As  will 
be  described  below,  the  elastic  strun  energy,  u,  depends  on  the  effective  elastic  modulus  of  the 
composite,  which  in  turn  depends  on  the  topological  distribution  of  the  constituent  phases.  By 
considering  two  extreme  topologies  corresponding  to  the  upper  and  lower  bounds  on  the  elastic 
moduli,  limits  on  the  magnitude  of  the  residual  stress  distribution  can  be  established.  The  Hashin 
bounds  are  used  in  this  work  as  the  microstructure  appears  to  be  random.  Lastly,  a  further  re¬ 
striction  on  the  stress  distribution  is  imposed  by  invoking  the  principle  of  maximum  entropy.  The 
values  for  the  resulting  stress  distribution  depend  on  the  topological  nature  of  the  phase  distri- 
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bution,  but  for  the  Al203-Zr02  system  the  upper  and  lower  bounds  appear  to  be  close  to  one 
another. 

To  begin,  the  residual  stress  field  needs  to  be  defined  and  from  it  the  stored  elastic  strain 
energy  expressed.  In  the  AljOa-ZrOj  material,  the  residual  stresses  are  a  result  of  both  the 
anisotropic  thermal  expansion  of  the  alumina  grains  and  the  difference  in  thermal  expansion 
between  the  alumina  and  zirconia  phases.  The  residuad  stress  in  any  particular  alumina  grain, 
whose  crystallographic  axes  are  oriented  at  an  angle  fi  to  a  fixed  coordinate,  can  then  be  written 
as  the  linear  superposition  of  two  parts: 

=  ^ik(s)  -h  <r“fc(x|n)  (2) 

where  denotes  the  stress  in  the  alumina  averaged  over  all  possible  crystallographic  orienta¬ 
tions  of  the  grains  and  therefore  corresponds  to  the  isotropic  component  (»)  due  to  the  presence 
of  ZrOj  grains,  and  o'®*  gives  the  orientation  dependent  part,  corresponding  to  the  anisotropic 
component  (a)  of  the  stress  distribution  due  to  the  thermal  anisotropy.  These  two  components 
can  be  considered  separately  in  the  following. 

Given  the  residual  stress,  the  elastic  strain  energy  can  be  expressed  [3]  by  the  relationship: 

2u  =  -  <  >=  2u'  -I-  2u®  (3 ) 

with 

where  and  Sz  are  the  elastic  compliance  tensors,  cj  and  c|  are  the  average  stress  free  strain 
tensors  of  the  AI2O3  and  Zr02  ph2ises,  c*  is  the  macroscopic  stress  free  strain  tensor,  and  = 
t'.4€l+i’z«7  (I’A  and  I'z  are  the  volume  fractions  of  AI2O3  and  Zr02)  [3].  The  anisotropic  component 
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of  strain  energy  is  given  by 


2u“  =  3 


(5) 


where  e^  =  3(2e2'  +  ef)  is  the  hydrostatic  component  of  cj,  auid  cf  and  ef  are  the  stress  free 
strains  in  the  alumina  grains  in  the  directions  of  transversal  isotropy,  and  in  the  3-direction 
(corresponding  to  the  a-  and  c-  crystallographic  axes  of  the  trigonal  lattice  of  sapphire).  The 
terms  5a  =  Sukk^  Sc  =  Ssskk  are  sums  of  components  of  the  elastic  compliance  tensor  and  e*  is 
the  isotropic  macroscopic  stress  free  strain  (assuming  no  texture). 

According  to  the  work  of  Kreher  and  Molinari  [9],  the  first  and  second  moments  of  the  stress 
distribution  in  a  heterogeneous  solid  can  be  quite  generally  expressed  in  terms  of  the  stored  elastic 
str2un  energy  in  the  stress  field  by  the  following  relationships: 


U,  <  <7,  >  6{€jk)r  * 

(6) 

/V 

<  (TfSSrtTr  >—  ^ 5( 5” j” ")  ^( ^•**"* )’’ 

(7) 

where  ^(c,>)r  and  SSr  describe  virtual  variations  of  the  stress  free  strain  and  elastic  compliance 
tensor  of  the  particular  phase  r. 

For  clarity  of  exposition,  the  evaluation  of  the  first  and  second  moments  of  the  stress  distribu¬ 
tions  are  considered  separately  below. 

2.1  Average  Stress 

The  average  stress  in  the  alumina  grains  can  be  expressed  as  the  stresses  in  equation  2  averaged 
over  all  the  alumina  grains  in  the  material.  This  is  given  by  the  first  moment  of  the  stress 
distribution  over  the  ensemble  of  grains: 


<  <TikpA  =  <  <T\k  >A  +  < 


(  V 


5 


where  <  <7{*  >.<  is  the  first  moment  of  the  isotropic  stress  field  and  <  is  the  first  moment 

of  the  distribution  for. orientation-dependent  part  of  the  stress  field  due  to  anisotropy  of  AI2O3. 

The  isotropic  part  of  the  average  stress  in  AI2O3  is  simply  given  in  terms  of  the  effective  modulus 
by  the  expression: 

<<^ik>A=  ^A^ik  (9) 


3  KaKz 
VA  Ka  -  Kz 


(  Kz  Ka-K‘\  , 
ir-  Ka-Kz) 


(10) 


with 


oa 


20,  -f  Qc 
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(11) 


AT  is  the  temperature  drop  over  which  the  residual  stresses  are  created  on  cooling.  From  force 
balance  condition,  the  average  stress  in  the  zirconia  grains  is  given  by 


vz  <  (r'iii  >z  >A=  0  ( 12) 

Since  the  zirconia  phase  is  cubic,  <  >7  is  the  total  average  stress  in  zirconia  grains.  The 

effective  modulus  K’  is  a  function  of  the  elastic  constants  of  the  two  phases,  and  their  volume 
fractions  as  well  as  the  topological  arrangement. 

.Assuming  that  the  microstructure  is  statistically  isotropic,  the  effective  bulk  modulus  A'*  is 
restricted  to  lie  between  the  upper  and  lower  Hashin  limits,  >  K"  >  K~.  The  upper  limit 
corresponds  to  a  self-similar  composite  sphere  topology  in  which  Zr02  inclusions  are  embedded 
within  a  .AI2O3  matrix.  The  lower  limit  corresponds  to  AI2O3  inclusions  embedded  within  ZrO^. 
The  Hashin  limits  may  be  expressed  in  terms  of  the  volume  fraction  of  AI2O3  by  the  equations  [3): 


K*  =  vaKa  +  I’zA'z  -  vavz 


{Ka  -  Kz? 

+  VzKa  +  ‘^I^Ga 


(13 1 
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(14) 


K"  =  vaKa  +  vzKz  -  vavz 


{Ka  -  Kzy 
va^z  +  vzKa  +  “i/ZGz 


where  K  and  G  are  the  bulk  and  shear  moduli  of  the  constituent  phases.  Whereas  the  effective 
Young’s  moduli  for  the  different  topologies  differ  by  leas  then  5  percents,  the  predicted  average 
stresses  corresponding  to  the  alternative  topologies  differ  considerably  more.  Using  the  upper  and 
lower  Hashin  limits  for  the  effective  modulus  in  equation  10,  the  average  stress  can  be  rewritten 
as 

sa  =  -3vz/C  (ctA  —  otz)^T  (15) 


with  (a)  for  the  upper  Hashin  limit. 


a:-*-  =  MKa  +  vaIKz  +  Z/AGa)'^ 


(16) 


(b)  for  the  lower  Hashin  limit, 

AC-  =  {vz/Ka  +  VaIKz  +  3/4Gz)-»  (17) 

.As  will  be  shown  in  figure  5,  the  two  topologies  are  characterized  by  significantly  different  average 
internal  stresses. 

The  upper  and  lower  Hashin  limits  for  the  effective  modulus  are  related  to  the  so-called  effective 
field  approximation  (EFA)  for  the  internal  stress  field.  In  this  approximation,  one  of  the  compo¬ 
nents  is  the  matrix  phase  while  the  other  is  the  inclusion  phase.  For  a  microstructure  of  random 
distribution  of  the  two  phases,  the  effective  medium  approximation  (EMA)  often  describes  the 
topology  more  satisfactorily  [3}.  In  contrast  to  the  Hashin  topologies,  interconnecting  networks 
of  grains  are  a  feature  of  topologies  for  intermediate  volume  fractions  in  the  EMA  method.  The 
effective  bulk  modulus  of  a  composite  according  to  the  effective  medium  approximation  (EM.Xi 
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fulfills  the  following  condition: 


and  it  normally  has  to  be  evaluated  numerically,  as  is  done  later.  As  will  be  shown  figure  5,  the 
average  stress  calculated  by  EMA  is  shown  to  have  values  intermediate  between  the  two  Hashin 
bounds. 

The  first  moment  of  the  anisotropic  part  of  the  stress  field,  that  due  to  the  thermal  misfit 
between  the  differently  oriented  alumina  grains,  can  be  expressed  as: 


<  fffjn  >=<  >=  -Gb{qc  -  <ar«)Ar  =<  > 


<  <^33!^^  >=  -"^Gb{oic  -  aa)Ar  =<  <T*  > 


with 


Ga 

5(3/rM+4G4) 


(20) 


This  expression  is  valid  provided  the  alumina  grains  are  randomly  oriented  and  the  elastic  anisotropy 
of  alumina  is  small,  which  it  is  [10]. 

Thus,  combining  the  above  equations,  the  average  stress  in  the  alumina  grains  along  the  a-  and 
c-  crystallographic  directions  can  be  written  as: 


<  >;4=  -3u^AC  (q,i  —  Ctz)Sr  +  -Gsicie  “  Cla)^T  (21) 

<  >A~  -^VzK  (q^  —  Qz)Ar  -  -Ge(Qc  “  0‘a)^T  (22) 


where  K.  has  the  value  given  in  equations  16  and  17,  appropriate  to  the  topological  description  of 
the  material  microstructure. 
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2.2  Stress  Fluctuations 


The  second  moment  of  the  stress  distribution  describes  the  spatial  fluctuation  of  the  residual 
stresses,  and  as  described  in  reference  [8]  is  characterized  by  the  width  of  the  fluorescence  peak. 
Using  the  separation  of  stresses  described  by  equation  2,  the  second  moment  can  be  expressed  [3] 
as: 

<  (Tii^aik  >=<  a\k  crik  >  +2  <  >  +  <  <y“k<y“fc  >  (23) 

If  the  cross-correlation  of  the  two  stress  fields  is  small  relative  to  the  other  two  terms,  it  can 
be  neglected  and  the  second  moment  becomes: 


<  <Tik(Tik  >=<  <rik  ffik  >  +  <  >  (24) 

In  order  to  derive  an  expression  for  the  second  moment  of  a  particular  stress  component,  we 
assume  that  the  second  order  stress  invariant  equals  the  sum  of  the  second  order  moments  of  the 
three  principle  stresses  defined  in  the  lattice  coordinate  of  a  given  crystal.  This  ignores  any  possible 
correlation  between  the  second  order  moments  of  the  shear  stresses  with  the  surroundings.  By 
making  the  further  assumption  that  the  magnitude  of  the  fluctuations  in  each  of  the  three  isotropic 
principle  stresses  is  the  same,  then  the  second  order  moments  may  be  written  as: 

<  >=<  >=<  >=<  (25  . 

These  values  of  the  second  order  moments  for  the  isotropic  components  <t,‘*  can  now  be  obtained 
explicitly  by  imposing  limits  to  the  phase  topology  using  the  Haishin  bounds  on  the  modulus. 
Using  equation  4  and  7,  the  upper  limit  is  given  by  the  upper  Hashin  bound,  namely 

C  =  3(^)‘/'AC+{q4  -  az)M  v2t, 
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The  lower  bound  corresponds  to  AljOa  as  the  embedded  inclusion.  Since  in  this  geometry  the 
stress  in  the  inclusion. is  purely  hydrostatic,  the  stress  fluctuations  are  necessarily  zero,  i.e., 

=  0  (27) 


Both  limits  are  shown  in  flgure  6. 

The  second  order  moments  for  the  isotropic  stress  can  be  refined  by  application  of  the  max¬ 
imum  entropy  principle  of  the  information  theory  [3].  The  calculation  takes  into  account  of  all 
possible  topologies  which  are  compatible  with  a  given  set  of  macroscopic  thermal-elastic  proper¬ 
ties,  i.e.,  effective  modulus  and  stress  free  stain.  Under  the  assumption  of  the  effective  medium 
approximation,  i.e.,  only  one  phase  forms  a  continuous  network  (composite  sphere  topology),  two 
functions  corresponding  to  the  upper  and  lower  Hashin  bounds  are  derived 

-  az)Ar  (28) 

-  az)AT  (29) 

These  two  functions,  evaluated  for  the  AljOa-ZrOj  composite,  are  plotted  in  figure  6,  where  they 
are  compared  with  our  experimental  data. 

.\  further  restriction  cac  be  imposed  by  applying  the  maximum  entropy  principle  to  the  stress 
variation  and  using  the  macroscopic  thermal-elastic  properties  derived  from  the  effective  medium 
approximation,  rather  than  those  given  by  the  Hashin  bounds.  It  can  be  shown  [3]  that  the 
function  calculated  this  way  lies  between  the  Hashin  bounds,  i.e., 


„EMA 


>  (7 


I- 

w 


.\s  indicated  by  the  curves  in  figure  6,  the  differences  between  these  predictions  are  small. 
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The  second  order  moment  of  the  orientation  dependent  stress  component  can  be  calculated  ap¬ 
plying  the  derivative  o.f  the  stored  elastic  energy  as  well.  Under  some  mild  assumptions  concerning 
the  neglected  elastic  anisotropy  and  correlation  between  the  different  grains  for  polycrystalline  sin¬ 
gle  phase  materiaJs,  an  analytical  expression  for  the  density  of  the  elastic  stored  energy  u  can  be 
derived  [9].  The  second  order  stress  moments  follow  immediately  by  derivatives  with  respect  to 
the  bulk  modulus  Ka  and  the  shear  modulus  Ga  of  the  isotropic  polycrystal: 

<  >=  3  <  s**  >  +  <  >  (30) 


with 


<  =  2  il2KAGA?  , _ 

BKa  90  {3Ka  +  4G^  )2  ^  "  o« )  at 


15^^  (ZKa+^Ga)^  (ac-Q,)A7’ 


mA 

+  »iki  4k  =  0 


<  =  j(<  >  -|  <  < 


(31) 

(32) 

(33) 

(34) 


Thus,  the  first  and  second  moments  of  the  residual  stress  distribution  are  dependent  only  on 
the  elastic  constants  of  the  constitutive  phases,  their  volume  fractions,  their  thermal  expansion 
coefficients  and  the  stress-free  temperature  (AT).  In  practice,  the  last  is  the  parameter  that  is 
generally  unknown. 


3  Experimental  Detaib  and  Results 

3.1  Materials 

Nineteen  composite  samples  with  Zr02  volume  fraction  vz  ranging  from  0.01  to  0.9  were  mea¬ 
sured.  The  samples  were  all  sintered  with  the  exceptions  that  the  three  with  f2^=0.025.  O.U' 
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and  0.9  which  were  hot  pressed  in  order  to  attain  full  density.  All  samples  were  subsequently 
air-annealed  at  1650*. for  10  hours  after  desiflcation.  The  grain  size  of  both  the  AI3O3  phase  and 
the  ZrO]  phase  were  measured  for  all  the  samples  and  the  results  are  plotted  in  figure  2. 

3.2  Measurement 

The  techniques  used  in  fluorescence  measurement  and  analysis  have  been  described  in  detail 
in  [8].  Briefly,  an  unmodified  optical  microprobe*  was  used  to  both  excite  the  fluorescence  and 
to  collect  and  analyze  the  resulting  fluorescence  spectrum  using  an  attached  spectrometer.  The 
514  nm  line  of  an  argon  ion  laser  was  used  to  excite  the  fluorescence.  The  fluorescence  signals 
were  collected  from  a  region  of  about  ZOfxm  diameter  in  size  as  defined  by  the  laser  beam  so 
as  to  probe  many  grains  simultaneously.  The  intensity  of  the  Ri  amd  Rj  fluorescence  lines  were 
typically  scanned  by  integrating  over  0.5  second  intervals  at  a  spacing  of  0.2  wavenumbers,  with  the 
intensities  being  recorded  under  computer  control.  The  collected  data  was  subsequently  analyzed 
with  curve  fitting  algorithms,  included  in  the  LabCalc  software  package^,  to  identify  the  position  of 
the  fluorescence  peaks.  All  the  measurements  were  made  with  the  samples  at  room  temperature. 
The  peak  shift  due  to  temperature  fluctuation  was  corrected  using  the  temperature  dependence 
coefficient  of  ~  0.14  cm"*/*C  [11].  Local  heating  due  to  the  incident  leaser  beam  was  avoided  by 
using  a  small  laser  power  of  ~0.5  mW  at  the  sample.  The  instrumental  shift  due  to  variations  in 
room  temperature  was  also  corrected  by  simultaneously  monitoring  a  characteristic  Neon  line  at 
14564  cm"*,  close  to  the  R  lines. 

3.3  Analysis 

Several  spectra  were  taken  at  different  spots  on  each  sample.  The  resulting  peak  shift  and 
broadening  are  plotted  ais  a  function  of  Zr02  volume  fraction  vz  in  figures  3  and  4.  Both  the 

'  Instrumentt  SA,  R«nutn  Microprobe  Model  T64000. 

^Galictic  Industries  Corp. 


12 


shift  and  broadening  were  measured  using  a  strain  free  sapphire  as  the  zero  shift  and  broadening 
standard. 

From  equation  1,  the  frequency  shift  of  the  Rj  line  from  a  volume  element  with  constant  stress 
distribution  within  a  AI3O3  grain  is 


At/ =  2n«<r,  +  nc<Te  (35) 

where  Ila  and  lie  piezo-spectroscopy  coefficients  for  the  a  and  c  directions.  They  have  been 

measured  recently  at  room  temperature  to  have  values  of  2.75  and  2.10  cm“'GPa“^,  respectively 
[12].  A  positive  shift  implies  tension  and  a  negative  one  implies  compression.  Following  the 
assumption  of  section  2,  we  can  write  the  residual  stress  in  AI3O3  as  the  sum  of  the  contribution 
from  the  thermal  mismatch  between  the  two  phases  and  from  the  anisotropy  of  AI2O3  itself.  The 
fluorescence  shift  can  similarly  be  decomposed: 

Aj/  =  Aj/*  +  Ai/‘  (36) 

and  we  <issume  that  the  stress  due  to  thermal  expansion  anisotropy  is  independent  of  the  volume 
fraction  of  ZrOj.  Using  the  piezospectroscopic  relationship  in  equation  35,  we  have 

<  Al/“  >=  He  <  <7*  >  -|-2na  <  <7*  >=  -(Ila  “  He)  <  <7®  >  (37) 

Thus,  the  intercepts  of  figure  3,  <  Aj/®  >,  gives  directly  the  average  stress  along  the  c-direction 
due  to  the  thermal  expansion  anisotropy  above. 

The  shift  due  to  the  average  stress  caused  by  Zr02  inclusion  <  ct’  >  is 

<  Ai/'  >=  (lie  +  2na)  <  a'  >  (3S' 
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A  least  square  fitting  of  the  data  upto  vz  —  0.8  using  a  second  order  polynomial  gives  < 


Av"  >=  —0.037  cm“*.and  <  Ai/*  >*:  —  6.0u^  +  2.2u^  cm"  .  The  shift  is  independent  of  zirconia 
volume  fraction  from  0.8  to  0.9  and  have  an  average  value  of  —3.4  cm~^  Therefore,  using  equation 
37,  the  smisotropic  components  of  the  average  stress  along  the  a  and  c  directions  of  AI2O3  grains 
aue  —28  and  56  MPa,  respectively.  The  average  isotropic  stress  component  in  the  AI2O3  phase  as 
a  function  of  vz  is  also  derived  from  data  in  figure  3  using  equation  38 


<  a*  >A=  —790v2  +  290 ul  MPa  for  vz  <  0.8  ;  <  a*  >a=  —450  MPa  for  vz  >  0.8  (39) 


The  average  stress  in  the  Zr02phase,  <  a*  >z,  is  obtained  by  using  the  force  balance  require¬ 


ment  (equation  12),  thus 


<  ff'  >^=  790  —  290  Vz  MPa 


The  peak  broadening  is  directly  related  to  the  width  of  stress  distribution.  To  obtain  an  explicit 
relationship,  we  assume  again  that  <7,  and  (Tg  are  (a)  independent  of  each  other,  (b)  independent 
of  <7',  and  (c)  they  all  follow  Gaussian  distribution: 


/  ax  ,  ax  /  ix  >)’l  >)’l  >) 

pMpi<^c  )Pi<^  )  =  -/Vexp - ^ -  exp - - - 


where  <7“  and  are  the  breadths  of  the  distribution  in  stre-'es  <7“  (or  <7“)  and  er'  respectively. 
The  peak  broadening  can  then  be  shown  [8]  to  consist  of  two  parts 


<  A:/^  >  =  V  <  Ai/“*  >  <  At/‘*  > 


<  Ai^“*  >  =  \/4n*  +  n*a 


(44) 


v'<  Ai/'a  >  =  (2n,  +  neX 

where  n/<  >  is  the  broadening  due  to  thermal  anisotropy,  which  is  independent  of  17.  and 

>/<  >  is  due  to  the  thermal  mismatch  of  the  two  phases. 

Fitting  of  the  data  (for  sintered  materials  only)  in  figure  4  using  parabolic  dependence  of 
on  vz  gives  y/ <  Ai/“*  >  =  0.89  cm~'  and  <  Ai'**  >=  5.09  V7  —  5.16u|  cm"*.  Therefore,  from 
equation  43,  44,  the  widths  of  the  distributions  are 

<t“  =  153  MPa  ;  <  =  ^669  u7  -679v|  MPa  (45) 


4  Discussion 

The  data  of  figures  3  and  4  provides  a  basis  for  a  detailed  comparison  between  the  experiment 
and  the  predictions  of  the  stochastic  stress  analysis  outlined  in  section  2  for  both  the  average 
stress  in  the  alumina  phase  and  the  width  of  the  stress  distribution  as  a  function  of  the  volume 
fraction  of  the  zirconia  phase. 

The  observed  frequency  shift  of  the  fluorescence  peak  (figure  3),  and  the  corresponding  average 
stress  (figure  5)  calculated  using  equation  38,  exhibit  a  number  of  interesting  characteristic  fea¬ 
tures  as  a  function  of  volume  fraction  of  zirconia.  For  volume  fractions  up  to  0.35  the  data  is 
approximately  linear,  then  its  curvature  becomes  evident  until  vz  ~  0.8  it  becomes  independent  of 
volume  fraction  of  zirconia.  Although  the  values  of  the  average  stress  can  be  fitted  very  accuratcl\- 
to  a  linear  dependence  on  zirconia  volume  fraction  upto  ~  0.35,  they  can  equally  well  be  fitted  over 
the  same  range  to  equation  15,  the  prediction  of  the  effective  field  approximation  for  the  a\-erage 
stress  in  the  alumina  phase  in  terms  of  the  effective  modulus  and  the  thermal  e.xpansion  mismatcii 
strain.  This  is  shown  in  figure  5  by  the  lower  dashed  line  corresponding  to  the  upper  Hashin 
bound  on  the  effective  modulus  and  for  a  temperature  change  AT  of  —  1210°C.  The  excellent 
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correspondence  between  the  data  and  the  upper  Hashin  bound  upto  at  least  35%  zirconia  is  con¬ 
sistent  with  the  microstnictural  observations  of  zirconia  particles  embedded  in  an  alumina  matrix. 
This  is  despite  the  fact  that  the  percolation  linut  for  randomly  distributed  zirconia  particles  is 
—  17%  [13],  suggesting  that  alumina  may  actually  “wet”  the  zirconia  phue.  The  departure  of  the 
data  at  the  intermediate  volume  fractions  from  the  prediction  based  on  the  upper  Hashin  bound 
indicates  that  the  interconnecting  networks  of  grains  start  to  dominate  the  microstructure,  and 
the  effective  modulus  increases  with  increasing  ZrOj  volume  fraction.  At  high  volume  fractions, 
zirconia  becomes  the  matrix  and  alumina  grains  become  inclusions.  It  is  therefore  expected  that 
the  average  stress  dependence  should  follow  the  lower  Hashin  bound.  To  better  describe  the  stress 
dependence  over  the  entire  volume  fractions,  the  effective  medium  approximation  (EMA)  is  also 
included  since  it  takes  into  account  all  possible  topologies,  providing  intermediate,  values  between 
the  two  Hashin  bounds  in  the  intermediate  volume  fractions,  but  approaching  th.  ^  bounds  at 
the  extreme  volume  fractions.  This  function,  calculated  numerically,  is  plotted  in  figure  5,  and  pro¬ 
vides  a  close  fit  to  the  experimental  data,  except  at  the  highest  volume  fractions  where  the  stress 
becomes  independent  of  volume  fraction  instead.  This  plateau  region  cannot  readily  be  explained 
in  terms  of  elasticity.  One  possible  explanation  is  that  at  the  high  volume  fractions  of  zirconia 
the  stress-free  temperature  is  significantly  lower  than  that  for  the  adumina-rich  composites.  An 
approximately  linear  dependence  on  volume  fraction  would  nevertheless  be  expected.  The  other 
possible  explamation  is  that  the  residual  stress  exceeds  that  for  plastic  yielding  in  the  zirconia 
surrounding  the  individual  alumina  grains.  From  figure  5,  the  average  stress  in  the  alumina  grains 
in  this  plateau  region  is  ~  450  MPa,  which  would  correspond  to  a  resolved  shear  stress  for  plaistic 
flow  of  640  MPa  in  the  zirconia  phase.  It  is  difficult  to  know  whether  this  is  reasonable,  since 
the  complete  information  on  the  yield  strength  of  cubic  zirconia  in  the  entire  cooling  temperature 
range  is  not  available.  This  is  consistent,  however,  with  the  measurements  of  Cheong  et  al.  [14^ 
of  the  yield  stress  in  cubic  zirconia  single  crystals  at  temperatures  above  1200® C.  They  found  the 
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yield  stress  of  cubic  zirconia  drops  significantly  as  the  temperature  increases  due  to  activation  of 
secondary  slip  systems^ 

The  value  of  the  temperature  change  used  to  fit  the  average  stress  in  figure  5  is  of  some  interests. 
First,  since  the  materials  were  sintered  at  1650‘C,  the  calculated  temperature  change  over  which 
thermal  stresses  developed  indicates  that  diffusional  relaxation  accommodated  any  stresses  for  the 
first  400  degrees  of  cooling.  This  is  consistent  with  the  observation  that  zirconia  is  an  effective 
densification  aid  for  alumina,  lowering  the  temperatures  necessary  for  densification.  Secondly, 
using  this  temperature  change  the  anisotropic  stress  components  in  AI3O3  are  evaluated  using 
equations  18  and  19  to  be  88  and  —44  MPa  along  the  c  and  a  axis,  respectively.  These  values  are 
comparable  to  those  directly  obtained  from  the  intercept  in  figure  3  (section  3.3). 

The  thermal  strain  AaAT  calculated  from  the  linear  portion  of  the  data  in  figure  5  provides,  in 
principle,  the  information  necessary  to  evaluate  the  predicted  values  of  the  second  moment  of  the 
stress  distributions  (equations  26-29).  Comparison  of  these  stress  fluctuations  with  the  measured 
broadening  data  involves  the  addition  of  terms  quadratically  using  equation  42.  To  do  this,  the 
data  in  figure  4  is  replotted  in  figure  6  together  with  the  upper  and  lower  Hashin  bounds  on  the 
broadening  calculated  using  equations  28  and  29,  and  43  and  44.  The  lower  Hashin  bound  on 
the  broadening  is  independent  of  the  zirconia  volume  fraction  because,  as  mentioned  in  section 
2.2.  this  would  correspond  to  isolated  alumina  grains  embedded  in  ZrOj  matrix  in  which  case  the 
stress  resulting  from  the  phase  mismatch  would  be  purely  hydrostatic.  The  fluctuations  are  then 
the  same  as  in  a  polycrystalline  alumina  matrix.  The  upper  Hashin  bound  on  the  broadening, 
whilst  also  bounding  the  data,  evidently  does  not  fit  the  data  well.  Much  closer  to  the  data  are 
the  bounds  on  the  broadening  calculated  by  imposing  the  additional  restriction  of  the  maximum 
entropy  principle  (equations  28,  29).  This  is  particularly  so  if  the  data  from  the  hot-pressed 
materials  (the  open  data  points)  is  ignored.  (  The  hot- pressed  materials  have  exceptional  grai.n 
sizes  compared  to  the  sintered  materials.)  Although  the  fit  between  the  data  from  the  sintered 
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materials  and  the  bounds  is  good,  more  detailed  comparisons  are  not  possible  at  the  present  time. 
This  is  primarily  because  the  observed  variation  in  peak  broadening  from  place  to  place  within 
the  materials,  especially  at  the  lower  and  higher  volume  fractions  of  zirconia,  is  larger  than  the 
difference  between  the  bounds.  The  experimental  variation  is  real  and  is  believed  to  be  due  to  the 
difficulty  in  sampling  a  sufficiently  large  region  that  is  statistically  homogeneous  at  the  extreme 
volume  fractions. 

The  above  comparison  can  also  be  carried  out  in  terms  of  the  stress  fluctuation  caused  by  the 
phase  mismatch  component,  alone.  Plotted  in  figure  7  aire  the  bounds  on  the  second  moment 
of  the  stress  distribution  calculated  using  the  thermal  strain  derived  from  the  first  moment  of 
the  stress  distribution.  The  upper  and  lower  Hashin  bounds,  and  their  modifications  obtained  by 
applying  the  restriction  of  maximum  entropy  are  indicated.  Also  shown  is  the  curve  corresponding 
to  the  least  square  fit  to  the  experimental  broadening  data  from  figure  4  (equation  45),  which  has 
values  close  to  but  smaller  than  both  of  the  bounds  from  the  information  theory.  This  is  not 
entirely  surprising  because  the  difference  between  the  two  bounds  are  extremely  small,  hence  any 
systematic  error  in  the  experimental  data  and  of  the  theory  introduced  by  simplifications  of  model 
can  easily  cause  discrepancies  large  compared  to  the  spacing  between  the  two  maximum  entropy 
bounds. 

Having  analyzed  the  experimental  data  it  is  now  possible  to  offer  an  assessment  of  the  validity 
of  the  stochastic  stress  analysis.  The  agreement  between  the  average  stress  determined  from 
the  fluorescence  frequency  shift  and  that  calculated  using  the  stochaistic  analysis  based  on  the 
upper  Hashin  bound  is,  up  to  about  35%  zirconia.  excellent.  However,  this  by  itself  is  not  a 
particularly  stringent  test  of  the  analysis  since  the  same  prediction  can  be  made  by  a  variety  of 
other  approaches.  For  instance,  according  to  Esfaelby's  calculations  [15]  the  interaction  energy 
between  centers  of  dilatation  is  dependent  only  on  their  volume  fraction  and  the  effective  elastic 
modulus.  On  the  other  hand,  by  using  the  thermal  strain  obtained  from  the  upper  Hashin  bound 
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fitting  at  relatively  low  volume  fractions,  the  effective  medium  approximation  produces  a  close  fit 
to  the  data  through  out  the  entire  volume  fraction.  This  suggest  that  the  stochastic  analysis  can 
predict  the  first  moment  of  the  stress  distribution  accurately,  provided  that  a  realistic  model  for 
the  microstructural  topology  is  used.  A  more  exacting  test  of  the  stochastic  analysis  is  the  second 
moment  of  the  stress  distribution.  The  experimental  data,  although  showing  greater  variability 
from  place  to  place  in  the  material  than  for  the  first  moment,  is  consistent  with  the  predictions  of 
the  stochastic  analysis  that  they  closely  fit  the  maximum  entropy  limit  of  the  Hashin  bounds.  With 
the  exception  of  the  hot-pressed  composites,  the  majority  of  the  data  exhibits  the  characteristic 
convex  dei>endence  on  volume  fraction.  Since  the  only  data  used  in  predicting  the  second  moment 
of  the  stress  distribution  are  the  single  crystal  elastic  constants,  the  zirconia  volume  fraction, 
and  the  thermal  strain  obtained  from  the  data  for  the  average  stress,  the  agreement  represents  a 
further  validation  for  the  stochastic  analysis. 

5  Concluding  Remarks 

The  thrust  of  the  majority  of  studies  of  the  elastic  behavior  of  composite  materials  has  hitherto 
been  in  the  determination  of  the  elastic  moduli,  or  equivalently,  the  thermal  conductivity,  as  a 
function  of  volume  fraction.  From  such  investigations,  bounds  on  the  elastic  moduli  have  been 
established  and  compared  with  experiment.  By  contrast,  in  this  work  we  have  determined  the 
residual  stress  distribution  in  a  series  of  two-phase  Al203-Zr03  composites,  rather  than  the  moduli, 
and  compared  it  with  the  predictions  of  a  stochastic  stress  analysis.  This  has  proven  possible  by 
using  the  technique  of  optical  fluorescence  spectroscopy,  which  provides  a  measure  of  both  the 
first  and  second  moments  of  the  stress  distribution  from  the  frequency  shift  and  broadening  of 
the  characteristic  R-line  fluorescence  from  Cr^'*’  dopants  in  aluminum  oxide.  The  predictions  of 
the  stochastic  stress  analysis,  outlined  in  section  2,  are  in  good  agreement  with  the  measurement  - 
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made  by  fluorescence  spectroscopy  indicating  that  confidence  can  be  placed  on  the  predictions. 
These  are  expressed  in  terms  of  the  single  crystal  elastic  and  thermal  expansion  constants,  the 
second  phase  volume  fraction  and  the  thermal  strain.  The  last  is  obtainable  from  the  volume 
fraction  dependence  of  the  average  stress,  determined  from  the  shift  in  the  fluorescence  frequency 
as  a  function  of  volume  fraction.  More  stringent  tests  of  the  stochastic  analysis  than  that  presented 
here  can  be  envisaged  but  they  will  require  both  a  better  method  of  experimentally  averaging  over 
large  volumes  of  material  in  order  to  establish  statistically  significant  data  as  well  as  methods 
of  malting  random  composites  in  which  scale-dependent  stress  relaxation  processes  do  not  occur. 
In  the  absence  of  such  developments,  the  experimental  measurements  of  the  first  and  second 
moments  of  the  residual  stress  distributions  compare  exceptionally  favorably  with  the  predictions 
of  the  stochastic  stress  analysis. 
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Figure  Captions 

1.  Scanning  electron  micrograph  of  three  of  the  alumina-zirconia  materials  studied  to  illustrate 
the  distribution  of  the  two  phases.  The  top  is  10%  zirconia,  the  middle  is  50%  and  the  bottom 
85%.  The  zirconia  phase  appears  as  the  brighter  of  the  two  phases  due  to  its  higher  atomic 
number. 

2.  The  average  grain  size  of  the  AI3O3  and  ZrOj  phases  as  a  function  of  zirconia  volume  fraction. 

3.  Observed  shift  in  R3  fluorescence  frequency  as  a  function  of  the  volume  fraction  of  the  cubic 
zirconia  phase.  The  full  data  points  are  from  the  sintered  materials  whereas  the  open  data  points 
are  from  the  hot-pressed  materials.  The  line  through  the  data  corresponds  to  a  least  squares  fit 
to  the  fluorescence  data  (section  3.3). 

4.  Observed  broadening  of  the  R3  fluorescence  peak  as  a  function  of  the  volume  fraction  of 
zirconia  phase. 

5.  The  data  for  the  average  stress  from  figure  3  compared  with  the  predictions  of  the  upper  and 
lower  Hashin  bounds,  and  the  EMA  on  the  first  moment  of  the  stress  distribution.  The  upper 
Hashin  bound  is  given  by  equation  15  with  the  modulus  function  K  of  equation  16  and  an 
effective  freezing  temperature  of  —  1210"C.  The  data  for  the  entire  volume  fraction  range  is 
plotted  in  the  top,  the  lower  volume  fraction  part,  from  which  the  thermal  strain  was  obtained, 
is  replotted  in  the  bottom  for  more  detailed  comparison. 

6.  Comparison  of  the  observed  fluorescence  peak  broadening  and  the  predictions  of  the 
broadening  based  on  bounds  derived  from  the  stochastic  stress  analysis.  The  broadening  is 
related  to  the  stress  distributions  using  the  piezo-spectroscopic  equations  42-44.  H-f  and  H-  are 
the  upper  and  lower  Hashin  bounds.  1+  and  I—  are  the  upper  and  lower  bounds  obtained  by 
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Applying  the  principle  of  maximum  entropy. 

7.  Bounds  on  the  second  moment  of  the  stress  distribution  calculated  using  the  thermal  strain 
derived  from  fitting  the  first  moment  of  the  stress  distribution  in  figure  5.  Also  shown  as  the 
solid  curve  is  the  second  moment  derived  from  the  least  squares  fit  to  the  experimental  data  for 
the  observed  broadening  of  the  fluorescence  peak  (equation  45). 


23 


0.2  0.4  0.6 

Zirconia  Volume  Fraction 


Average 


-400 

-600 

-800 


.  .  .  ."-1 

0.2  0.4  0.6  0.8  1 

Zirconia  Volume  Fraction  v 


TIm  American  Society  of 
Mechanical  Engineers 


ft»print»d  From 
AMD>Vol.  181.  Experiments  in  Smart 
Materiais  and  Structures 
Editor:  Kytmg^uk  Kim 
Book  No.  H00888  - 1993 


STRESS  MEASUREMENT  USING  OPTICAL  FLUORESCENCE 

Qing  Ma  and  David  R.  Clarke 
Department  of  Materiais 
University  of  California 
Santa  Barbara,  California 


Abstnet 

A  y«s«rs/  m€tktdoUn  itttrmmmt  la  mafenaJi  tad  tiraciant  w  dtaehped  kttti  «•  Ikt 

ftat’^ttctratetpic  froftrUu  tf  erftlaJi.  GtaertI  rtkiitatkipt  Utwtea  tkt  mtttartd  flatttteeact 
Hat  aii/b  tad  At  ttrtM  tltU  art  ptattaitd.  Aa  tapernatal  mtatahaf  At  ruidatl  tfrttttt  la  ta 
tmMdtd  atpfkirt  /Utr  it  attd  U  Hhttratt  ttau  tf  At  aAmtayc*  a/  At  tfUcal  auAtd  tad  Ht  ate 
ftr  fTatia§  auerotetpie  ftaimrtt.  Tie  ttptrmtat  aita  tktwt  Att  ttntt  itattr  eta  be  dtttrmiaed 
detpift  At  feet  Ati  At  Jherttetact  tkifi  it  a  aca/ar  taeatitg. 


1  Introductran 

As  technology  enables  the  development  of  more  advanced,  and  often  more  complicated  mate¬ 
rials  and  stmctures,  it  also  demands  new  stress  measnrement  technk)ues.  Very  ^ten,  it  is  also 
required  to  measure  the  stress  state  with  high  spatial  resolution.  For  example,  it  is  desirable  to 
know  the  localized  stresses  associated  with  the  domain  structure  and  their  variations  due  to  do¬ 
main  switching  in  ferrodectric  materials,  because  these  stresses  can  affect  the  device  reliability 
by  inducing  crack  formation  ann  growth  under  electrical/medranical  loads.  Cmiventional  stress 
measurement  tedmiques  based  o^  diffraction  methods,  such  as  X-ray  and  neutron  di&action, 
although  accurate,  are  incapab^<:  of  such  selective  measurements,  wUch  usually  requires  both 
an  imaging  csq>ability  and  Ugh  spatial  resolution. 

The  development  of  the  optical  fluorescence  method  was  partially  motivated  by  such  de¬ 
mands.  Many  materials,  especially  ceramics,  are  known  to  possess  fluorescence  if  containing 
i^ppropriate  ^pant  ions.  The  duiacteristic  flutvescence  lines  shift  systematically  with  stresses 
(piezo-spectroscopic  effect),  and  therdore  can  be  employed  for  stress  measurement.  In  fact,  the 
frequency  shifts  d  the  flwxescence  lines  of  ruby,  generated  by  the  Cr^  impurities,  have  been 
used  to  nxmitor  the  pressure  in  diamond  cells  since  the  early  1970s  (Forman,  1972). 

Recently,  we  have  significantly  extended  this  technique  and  applied  it  for  studying  the  me¬ 
chanical  pr^>erties  of  ceramics  and  composites  (Ma  and  Clarke,  1993a,  b,  Ma  et  al,  1993a,  b). 
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By  using  sn  optics!  microscope  to  inspect  the  sample,  and  to  focus  a  laser  beam  on  a  small 
i^on  of  interest  to  excite  the  fluorescence  signals,  selective  nteasurement  of  stresses  can  be 
easily  achieved  with  spatial  resolution  determined  essentially  by  the  diffraction  limit  of  the 
laser  beam.  In  this  paper,  we  will  first  review  the  basic  aspects  of  the  fluorescence  method  and 
then  illustrate  its  utility  with  an  experiment  in  which  the  residua!  stress  in  a  sapphire  fiber 
embedded  in  a  cubic  zirconia  matrix  was  measured. 


2  Fluorescence  Method 

2.1  General  Formulations 

Fluorescence  lines  of  materials  are  usually  generated  horn  the  electronic  transitions  of  impurity 
i<ms.  When  a  stress  is  ^pUed  to  a  material,  the  resulting  strain  changes  the  crystal  field 
around  a  impurity  ion  and  therefore  shifts  the  electronic  levels  of  the  ion.  A  fluorescence  line 
shifts  if  the  cmresponding  initial  and  final  electraiic  leveb  of  the  radiative  transition  shift  with 
different  amount  of  energies.  The  relationship  between  the  applied  stress  and  the  resulting 
frequency  shift  depends  on  the  detailed  electronic  structure  of  the  ion  in  the  particular  lattice, 
and  ther^ore  is  complicated.  However,  for  the  purpose  of  stress  measurement,  it  is  sufficient 
to  establish  phenomenological  relationships  and  calibrate  them  experimentally. 

Conrider  a  crystal  subject  to  a  homogeneous  stress  (in  crystallogr^hic  frame  of  refer¬ 
ence),  the  frequency  shift  of  a  fluorescence  line  of  a  impurity  ion.  At',  as  a  function  of  the 
applied  stress,  can  be  expanded: 

Ai/  **  UijO-j  +  ■  •  •  (1 ) 

Since  the  second  and  higher  order  terms  are  usually  much  smaller  than  the  linear  terms  for 
most  ^plicatioitt  when  the  stress  levels  are  not  unusually  high,  they  are  often  neglected  from 
the  analysis,  therefme, 

AV  »  HijOy  (2) 

Here  By  are  the  i»ezo.spectrasoopic  coefficients  relating  shift  to  the  stress  a-j  in  the  crystallo- 
gr^thic  frame  of  reference.  As  previously  pointed  out  by  Grabner  (1978),  Hy  fmm  a  second  rank 
tensor  because  the  above  relationship  must  be  covariant,  and  due  to  the  symmetry  the  stress 
tensor,  oy  »  erji,  the  11  tensor  is  also  S3rmmetric.  A  major  limitation  of  piezo-spectroscopy  is 
evident  from  this  equatirm.  Unless  the  fiuotesdng  species  emits  at  several  different  frequencies, 
so  that  a  set  of  simultaneous  equations  (one  for  each  transition)  can  be  generated,  there  may  be 
an  insuflkient  number  of  equatirms  to  solve  for  each  of  the  maximum  of  six  independent  stress 
compmients  needed  to  completely  describe  the  stress  state.  Further,  the  rite  symmetry  of  the 
flwnesdng  ion,  whidi  governs  the  symmetry  properties  of  the  11  tensor,  may  also  preclude  the 
detemunarion  of  aU  the  independent  stress  components. 

Vinth  the  exceptions  that  the  rite  point  symmetry  group  belcmgs  to  either  monoclinic  or 
triclink  systems,  the  11  matrix  is  diagonal  in  the  orthogonal  crystallographic  frame  of  reference 
(Nye,  1967).  Consequently,  the  frequency  shift  is  insensitive  to  pure  shear  stress  as  far  as  the 
linear  terms  are  concerned.  Further,  if  the  rite  tymmetiy  belongs  to  one  of  the  uniaxial  systems 
(tetragmial,  hexagonal  or  trigonal),  we  have  Hu  «  IIss  (axes  1  and  2  are  perpendicular  to  the 
principle  synometry  axis),  as  in  the  case  <ff  Cr^  in  sapphire,  only  the  summation  oji  -f  can 
be  determined,  not  oj,  and  vjg  individually,  even  with  multiple  fluorescence  lines  available.  If 
the  symnoetry  is  cubic,  such  as  for  substitutional  ions  in  MgO,  the  n  matrix  reduces  further, 
Hn  w  H]]  w  Hss,  hence  only  the  trace  of  the  stress  tensor  can  be  obtained. 

While  the  pieso-spectroso^ic  coefficients,  Ily,  are  usually  calibrated  in  the  crystallographic 
coordinates,  it  is  sometimes  mme  convenient  to  represent  the  stress  in  other  coo^nates  (e.g., 
sample  coordinates).  The  stress  on  the  crystal  structure  will  then  be  given  by  the  ^plied  stress 
components  resolv^  onto  the  crystallographic  axes,  i.e.,  through  an  orthogonal  coordinate 
transformation: 

v'j-OitafiVu  (3) 


28 


*T, 

*T2 

*T2 

*T, 


*A2 


I 

I 

•  Phonon 

f 

Excitation 
19435  cm'' 

Ri  R2 

cm’' 

nm 

eV 

25000 

400 

3.1 

21000 

480 

2.6 

18000 

550 

2.2 

15160 

659.6 

1.88 

14448 

692.1 

1.793 

14419 

693.5 

1.790 

0 

0 

Pifiif*  1:  Lawr  dactMoic  mtufy  kwdt  of  •  Cr*^  too  in  tapphin  (Ndioii  and  Storge,  IMS),  deaignated  by 
th«r  wpiaaiMatioai  in  tht  Oa  (Kup. 

where  Oij  is  the  tnosformntioii  mstrix  relating  the  frame  of  reference  in  whidi  the  applied 
stresses  are  defined  with  respect  to  the  ciystallogra{duc  axes  of  the  lattice.  Thus,  the  fre¬ 
quent  shift  ci  a  fluonsoenoe  line  in  a  fluorescing  crystal  oriented  at  an  arbitrary  angle  to  a 
superimposed  stress  is  given  by  the  following  tensoriaJ  rdation: 

^umtliiOn^iiOia  (4) 

For  oertain  site  symmetriei  described  above,  the  11  matrix  is  diagonalized,  and  the  above 
equatioo  can  then  be  expressed  explicitly: 

Av  *  nii(oi]  +  ojs  +  ff») 

+  (Iln  — nn)(«j,on-bo|,02j  +  0jj0sa) 

+  (npi  —  nii)(a3]  Oil  +  ^») 

+  2(nss  —  Dll  )(asia32  ou  +  eaittn  ow  +  osjoss  o») 

+  2(n»  —  nii)(asia32  oij  +  032033  033  +  031033  031)  (5) 

2.2  FhioreaoeBcc  of  Chnmaiiim  lone  in  S«]^hire 

Chroatam  in  sapphire  is  a  widl  known  dirainophore,  9ving  rise  to  the  characteristic  red  flu¬ 
orescence  of  rult-  The  fluoresoenoe  oripnates  from  dectnmic  transitions  in  the  chromium 
ion  produced  by  ^ropriate  exdtatioo  as  illustrated  in  figure  1.  Electrons  excited  from  their 
grouiMl  state,  ^As,  can  jump  to  the  ’E  level  by  emitting  phonon.  The  chrwnium  ion  can  substi¬ 
tute  for  any  of  the  aluminum  ions  in  the  cmundum  crystal  stnictuxe  of  sa^hite  with  the  result 
that  they  aw  octahedrally  coordinated  to  ndghboring  oxygen  and  aluminum  ions.  The  octahe¬ 
dron  is  not  regular,  but  elongated  along  the  c-axis  of  the  lattice,  resulting  in  splitting  (rf  the 
energy  levds  by  0.004  eV,  in  turn  giving  rise  to  the  two  distinct  radiative  transition  lines  in  ruby, 
Ri  and  Rs  at  1.790  and  1.794  eV  respectively.  Since  the  Ci**  ion  site  has  3-fold  rotational  axis 
almt  the  o«ds,  Dn  mHu  ^  Has.  For  simplicity,  we  will  use  the  notatirms  Iln  s  IIss  w  lU 
and  lls3  w  He  for  both  the  Ri  and  R3  lines.  The  pieso-spectroecopic  coefficients  TI,  and  II,  were 
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fraqiMBcr  shifts  aa  a  function  of  stress. 

pievkmaly  measured  (Fdier  aad  Sturge,  1968,  K^lyanskii  and  Pizbevuskii,  1962)  for  Ri  to  be 
3.2  and  1.4  an~’GPa'',  and  for  Rj  to  be  2.7  and  2.15  cm~’GPa~*,  respectively*.  A  positive 
shift  inqrUei  tension  and  a  negative  shift  im{dies  compresskm.  For  illustration,  tbe  shifts  of  the 
Ri  and  Rs  peaks  when  a  single  crystal  sapphire  was  compressed  along  its  c>axis  are  shown  in 
figure  2.  The  sharp  neon  line  it  used  for  fr^uency  calibration. 


2.3  Procedures  for  Optical  Fluorescence  Measurement 

We  use  aa  unmodified  <g>tical  microprobe’  to  both  excite  the  fluorescence  aad  to  collect  and 
analyse  the  resulting  ^ectrum  using  an  attached  spectnaneter.  The  optics  of  the  microprobe 
is  illustrated  in  figure  3.  An  argon  ion  laser  with  a  number  of  possible  line  frequencies  is  used 
to  excite  the  fluorescence.  To  start  a  measurement,  a  i^km  of  interest  on  the  sample  is  first 

*  A  n«m  wmmmmmm  pi*  aaS  CMm,  ItSS)  vwmlidb'  Ib***  »«bm,  hot  ah*  i««*«l«d  nMll  dHhnaeit  wb**  nm. 

wm  uwfcdla  p«*p*udad»  Simti—  Im  Uw  b«*ri  plia*  oaniBig  tht  d»  S-W  lywtry  mty  b«  htolwn.  Heiwcvtr.  ler 
•h*  (f  eaebcii]Mh*i*  lOhuMM*  Madid  bw*. 

*bMmmaH*  SA.  OpUal  Menpiab*  M*de  TSeMO. 
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DOUBLE  CONCAVE  HOLOGRAPHIC 

VIEWING  GRATING  MONOCHROMATOR  FILTER 


FigBM  i:  S^HMtic  diacnm  of  the  optical  mictoprobe. 


adcctod  whh  the  atudtad  optical  microacope,  then  the  laaer  beam  i»  fociued  to  a  spot  on  the 
adectad  faatuie  and  the  flnorcaoenoe  ptodn^  ia  odlected  and  analyzed.  The  intensity  of  the 
fltaonaeenoe  lines  aie  typcally  scanned  by  integrating  over  0.5  second  intervals  at  a  spacing  of 
0.2-0.d  wavenumbers,  with  the  intensiticB  being  recorded  under  computer  control.  The  collected 
data  are  subeequently  analysed  with  curve  fitting  algorithms,  included  in  the  LabCalc  software 
P«biCe*.  to  idntify  the  position  of  the  fiuorescence  peak.  By  using  objective  lenses  of  SOx 
and  lOOx  magnifyiag  powers  miniinum  spot  nses  of  ~  5pm  and  ~  1  -  2pm  can  be  produced. 
Larger  probe  saes  can  be  farmed  by  using  lower  magnification  lens  or  ^  partial  d^ocusing. 
Abo,  by  appropriate  dioioe  of  the  coUection  ^>ertnre  size,  the  attainable  axial  resolution  could 
be  varM  by  t^  ooUection  optics  6om  ~  10pm  to  ~  200pm. 

It  b  known  that  fluorcMenoe  lines  shift  systematically  with  the  change  of  temperature  ( Wun- 
der  and  Schoen,  1961).  For  the  R|  and  Ri  lines  of  ruby,  their  frequencies  shift  to  smaller 
wavenumbers  with  incieasiag  temperature  as  shown  in  figure  4*.  The  linear  fitting  coefficients 
are  -0.138  and  -0.132  cm~V*C  for  the  Ri  and  R}  lines,  respectively.  Therefore,  local  beating 
due  to  the  laser  exdtatkn  bew  has  to  be  avoided.  Thu  b  achieved  1^  using  low  laser  intensity 
sudi  that  no  line  drift  can  be  measured  when  the  spot  site  u  systematically  decreased  (thereby 
increasing  the  power  density  incident  on  the  sample),  indicating  that  no  significant  heating  of 
the  probed  vohnne  occurred.  Variatkms  were,  however,  noted  with  variations  in  room  temper¬ 
ature  and  so  oorroctiont  were  made  to  the  peak  shifts.  Thu  involved  correctioiu  for  the  change 
in  room  temperatnie  using  the  above  temperature  dependence  d  line  frequency  and  changes  of 
the  spectreineter  dimennons  at  a  result  <rf  thermal  expannon.  The  latter  can  be  corrected  by 
riamhaneoesbr  monitoring  a  diaracteristic  neon  calibration  line  dose  to  the  fluorescence  lines 
of  interest,  as  iOnstrated  in  figure  2. 

«Pn  lUyarf  la  wSsa— Uw  wfcS  J.  Il«.  OCSB. 
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Figuc  4:  FtaqiMScy  ahilte  of  Rt  an4  Ri  tiae*  a*  •  fvaetioa  of  tempctatiue. 

3  Determination  of  Residual  Stresses  in  a  Embedded  Sapphire 
Fiber 

To  illustrate  the  ^iplication  of  the  fluoiesceoce  method  we  have  measured  the  residual  stresses 
in  Si^hikaa  sapi^hire  fibers  embedded  in  a  cubic  sircooia  matrix.  The  fibers  have  sufficient 
levds  of  chronuum  impurity  to  produce  suffidently  intense  fluorescence  lines  Rj  and  Rs  of  ruby 
when  excited  by  as  atfoo  ioa  laser.  For  the  problems  of  a  c-axis  fiber  embedded  perpendicular 
to  the  OQoqxmte  surface,  the  axiaymmetric  stress  field  in  the  fiber  is  diagonalized  with  the 
radial  stress  equal  to  the  tangential  stress.  Therefiwe,  the  two  independent  stress  components, 
the  radial  stress  o’r  and  the  axial  stress  o,  can  be  determined  by  solving  the  following  two 
nmultaneous  equatkms  from  the  Ri  and  Ra  lines: 

Ani  s  2n«i0v  *1*  HciOi  (6) 

Asi  *  2n«39r  +  11*10,  (7) 

sdiere  the  subscripts  1  and  2  denote  R|  and  Ri  Hnes,  respectively. 


8.1  The  Through*!^  w  Depth  Profiling  Method 

The  method  of  measuring  subsurface  stresses  in  a  transparent  fiber  was  developed  in  a  earlier 
work  (hfa  and  Clarke,  1991b)  and  is  described  briefly  in  the  following.  As  illustrated  in  figure  5, 
the  fluoicsoenoe  from  depth  z  bebw  the  top  surface  is  collected  by  focusing  a  laser  beam  through 
the  fiber  top  surface.  By  focusing  to  succesnvely  greater  depths  and  measuring  the  frequency 
shift,  a  ]wofile  along  the  length  ^  the  fiber  it  obtained.  Bewuae  of  the  small  depth  of  field  of 
the  lens,  the  fluorescence  excited  will  be  collected  from  only  a  small  volume  on  either  side  of 
the  focal  plane.  The  measured  frequent  shift  £iv{z)  is  hence  a  weighted  average  of  that  within 
the  effsctive  excitation  volume.  Therefore,  the  actual  depth  profile  of  the  frequency  shift  has  to 
be  decouvohited  from  that  measured  by  using  the  depth  of  field  function  of  the  microprobe.  As 
described  previously  (Ma  and  Clarice,  1993b),  the  d^th  of  field  function  near  the  top  surface 
can  be  measured  systematkally  moving  the  plane  of  focus  from  a  position  above  the  top 
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Figan  S:  Sthwnatic  diapam  of  tha  optical  arraactmaat  «aad  to  eoUact  flnorttccBcc  from  a  region  in  an  em- 
baddad,  traaaparaat  ftbar  bdo*  Ha  r  tlaea.  Whan  tha  axeitatioa  laaar  ia  focuaad  at  a  ^atanea  a  balow  tha  top 
antihea,  only  tha  aigaal  ganarated  in  a  amatt  vehima  above  and  baiow  tha  focal  piaae  ia  affKtivaly  eoUactad. 


mifaoe  to  «  position  boloir  tbe  toiface  whHe  recordiag  the  Awnrescence  inteasity  I{z).  The 
depth  ci  fidd  furction  near  the  t^  surface  is  then: 


P(r)» 


ds 


(8) 


For  illustratkm,  the  depth  of  field  function  is  plotted  in  figure  6  for  our  optical  microscope  using 
a  40/ M  water  immecnon  lens  and  a  50pm  collection  ^>erture. 

Having  estahliihed  the  depth  fidd  function,  tbe  measured  shifts  can  then  be  related  to 
the  true  shift  profile  Av(r)  by: 


m-.C;Ay(»-l-t.)p(«)du 

jr,3(u)iu 

la  order  to  evaluate  the  deconvolution  it  was  assumed  that  the  true  shift  Ai'(r)  is  a  polynomial 
function  of  x  but  with  unknown  coefficients.  By  fitting  the  polynomid  function  convoluted 
with  the  depth  of  fidd  fimetioo  udng  equation  9  to  the  measured  distribution,  the  values  of 
the  coefficients  were  obtained.  The  initially  assiuned  function,  with  .the  coefficients  obtained 
by  fitting,  is  then  regarded  as  the  deconvoluted  shift  ^v{g).  The  accuracy  d  this  method  has 
bm  preWoody  demonstrated  (Ma  and  Claihe,  1993b)  and  can  be  judged  later  in  this  paper  by 
viewing  figure  7,  where  the  dashed  curves  are  the  true  shift  Av  (the  solid  curves)  convoluted 
fay  tbe  dq>th  of  fidd  function.  They  have  exceOent  fit  with  the  measured  data  in  generd. 

3.3  Shift  smd  StiwM  ProSlna  Along  End>cddcd  Sapphire  Fibers 

Appiying  the  dq>th  profiling  method,  the  residnd  stresses  along  a  c-axis  sapphire  fiber  embed- 
in  a  cubic  sitcania  matrix  were  measured.  The  diameter  of  the  fiber  is  about  130pm.  The 
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Fifuic  6:  The  depth  of  field  function  of  the  microprohe  need.  It  detchbc*  the  telntive  collection  efficiency  nc  a 
function  of  dietanee  fiom  the  focal  plane.  Detaik  ace  (iven  in  the  text. 

sapphire  fiber  was  perpendicular  to  the  composite  top  surface,  which  was  polished  to  obtain  an 
i^tical  finish  at  the  fiber  end. 

The  measured  and  the  decoovoluted  firequency  shift  as  a  function  of  depth  into  a  sapphire 
fiber,  as  well  as  the  convoluted  function  ate  plotted  in  figure  7.  Since  cubic  sirconia  has  thermal 
expansion  coefficient  larger  than  that  of  sapphire  (TaUe  I),  the  residual  stress  devdoped  by 
tbcnnal  mismatch  in  the  fiber  should  be  compressive  in  both  radial  and  axial  directions.  This  is 
consistent  with  the  large  negative  shifts  observed  when  the  probe  is  focused  deep  in  the  interior 
of  the  specimen. 


Table  I.  Properties  of  Fiber  and  Matrix. 


Property 

Sapphire  Fiber 

ZrOr  AI2O3  Matrix 

E{GP\) 

434 

220 

a  (10-*/*C) 

8.3  (a),  9.0  (c) 

11.0 

The  stress  distributions  calculated  from  the  equations  6  and  7  are  plotted  in  figure  8.  The 
axial  and  the  radial  stresses  ^proach  cemstant  values  of  1180  and  440  MPa,  respectively, 
toward  the  interior  the  cemyrosite.  The  axial  stress  is  close  to  zero  near  the  top  surface  as 
force  balance  condition  requires.  We  recognize  that  equations  6  and  7  are  not  strictly  valid  near 
the  fiber  end,  so  that  the  stress  values  derived  are  not  reliable  within  distances  to  the  surface 
comparable  to  the  fiber  diameter.  This  apparently  gives  rise  to  the  small  positive  values  for 
the  axial  stress  near  the  surface.  The  radial  stress  also  decreases  its  magnitude  approaching 
the  fiber  end,  due  to  the  bending  effect  associated  with  the  end. 

It  is  interesting  to  compare  the  interior  residual  stress  values  obtained  above  with  those 
calculated  using  dastic  solutions.  Since  the  s^phire  fiber  has  relatively  large  thermal  expansion 
anisotropy,  exact  expressiwrs  frv  the  stresses  are  lengthy  and  not  particularly  informative. 
However,  since  the  matrix  is  considerably  more  compliant  than  the  fiber,  Pmsson’s  effect  can 
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Fifun  7:  Tbe  Bee  ihift  depth  prcAlci  lor  the  (•)  Ri  and  (b)  R*  line*  (z  =  0  at  the  foiface)  in  a  e^axit  sapphire 
Sbcr  embedded  in  a  cubic  snconia  matzix. 
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Hfuie  8:  The  radial  and  axial  ttreases  at  a  functioa  of  depth  along  the  sapphire  fiber  determined  from  the 
fluotaseence  data. 
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b«  neglected  and  the  following  relations  are  obtained: 

o,  =  (10) 

a,  =  -£/(a„-Oe)AT  (11) 

Using  the  materials  parameters  listed  in  Table  I,  the  effective  temperature  difference  calculated 
using  the  radial  and  axial  stresses  are  AT,  =  1120‘’C  and  AT,  =  1360*’C.  The  difference  between 
AT,  and  AT,  indicates  that  the  system  did  not  behave  elastically  during  the  cooling  period  from 
the  processing  temperature.  It  is  conceivable  that  the  stress  relaxation  in  the  axial  direction 
is  less  effective  because  the  fiber  can  not  slide  over  long  distances,  so  that  the  temperature 
difference  calculated  from  the  axial  stress  is  closer  to  the  processing  temperature  of  1500°C. 


4  Concluding  Remarks 

The  preceding  analyses  and  experiments  illustrate  how  the  technique  of  fluorescence  spec- 
trosct^y  can  be  used  to  make  measurements  of  commonly  arising  types  of  stresses  in  materials. 
This  is  possible,  despite  the  fact  that  the  fluorescence  shift  is  a  scalar  quantity  whereas  the  com¬ 
plete  stress  state  is  tensorial,  because  considerable  simplification  of  the  equations  (equation  ■*) 
can  be  made  in  many  practical  situations.  For  instance,  it  was  demonstrated  (Ma  and  Clarke, 
1993a,  c,  Ma  et  al.,  1993c)  that  by  taking  advantage  of  the  random  orientation  of  polycrys- 
talline  ceramics,  simpliflcaiions  can  often  be  made  concerning  both  applied  stresses  and  residual 
stresses.  In  other  cases,  stresses  can  be  applied  at  prescribed  orientations  to  the  crystallographic 
axes  of  single  crystals,  where  the  symmetry  of  the  stress  held  provides  additional  relationships 
between  stress  components  and  therefore  enables  determination  of  the  stress  tensor.  This  was 
illustrated  in  the  case  of  residual  stress  in  embedded  fibers,  where  two  independent  stress  com¬ 
ponents,  the  axial  and  radial  stress  were  determined  by  two  piezo-spectroscopic  relations  from 
Rt  and  R}  lines.  This  experiment  also  provides  an  excellent  example  of  the  advantage  of  an 
optical  method.  A  single  fiber  can  be  easily  located  and  inspected  by  optical  microscope,  and 
probed  using  piezo-spectroscopy.  In  addition,  taking  advantage  of  the  fact  that  sapphire  is 
transparent  optically,  it  was  able  to  measure  stresses  from  the  the  deep  interior  of  the  sample 
and  therefore  obtaining  the  true  values  of  the  residual  stress  in  the  bulk. 
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